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A mixture of N“ and N* gas has been bombarded with 
deuterons of 1 Mev and the disintegration particles ob- 
served at right angles to the direction of the deuteron 
beam. The ranges of the disintegration particles were 
measured by means of a shallow ionization chamber and 

amplifier. The shortest measurable range was 2.1 
cm (15°C, 760 mm Hg). The reaction N'4(d, a)C™ gave 
rise to one group of alpha-particles at 5.25 cm range with a 
Q value of 7.54 Mev. Three groups of alpha-particles 
(11.97 cm, 6.54 cm and 3.47 cm range) were observed 
arising from the reaction N“(d,a)C*%. The Q values 
corresponding to these three groups are 13.39 Mev, 9.02 
Mev and 5.77 Mev, respectively. From the reaction 
N¥(d, p)N!* two valid groups of protons were observed 
having ranges of 90.5 cm and 21.02 cm with Q values of 


8.51 Mev and 3.15 Mev. The ranges given above are from a 
single measurement and the Q values are the average of the 
results of several measurements. A group at 66 cm pre- 
viously reported due to this reaction was found to be 
due to boron contamination. No proton groups were found 
which could be attributed to the reaction N**(d, p)N", 
although a radioactivity of half-life 9.5+1.0 sec. was 
found to result from deuteron bombardment of the mixture 
of N“ and N". This radioactivity was not found when N** 
alone was bombarded. A method of range measurement 
is described which allows accurate corrections to be made 
for the effect of chamber depth and of the bias on the grid 
of the counting thyratron on the output of the pulse 
amplifier. 


INTRODUCTION 


INCE the careful work of Cockcroft and 

Lewis (1936)! very little work has been done 
on the measurement of the energy of the heavy 
particles emitted when nitrogen, carbon and 
oxygen are bombarded with high energy deu- 
terons. Such heavy particle measurements are of 
particular interest in that they can give very 
accurate information about energy levels in the 
residual nuclei, although selection rules quite 
often prevent this method from giving informa- 
tion about all of the energy levels that one could 
on energy considerations expect to be involved. 
The present investigation is limited to a study 
of the (d, p) and (d, a) reactions with N™ and 


J. D. Cockcroft and W. B. Lewis, Proc. Roy. Soc. 154, 
261 (1936). 


N™ as the target nuclei. The results given in this 
report are based on more measurements with 
better corrections than those reported pre- 
viously** and are believed to be more accurate. 


EXPERIMENTAL ARRANGEMENTS 


The source of the high energy ions used in this 
experiment was the Cornell cyclotron. The beam 
of ions could be directed outside the cyclotron 
chamber to a position 14 to 30 cm from the 
chamber wall. As is shown in Fig. 1 the target 
cell is bolted to the exit port of the cyclotron. 


' The deuterons enter the target cavity through 


an aluminum foil at the end of the foil tube. 


1939) G. Holloway and B. L. Moore, Phys. Rev. 56, 705 
? M. G. Holloway, Phys. Rev. 57, 347 (1940). 
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The foil is attached to the supporting gridwork 
with Glyptal. The target gas is admitted to the 
target cell through the tube at the left and the 
pressure of the gas is measured by means of a 
Bourdon type vacuum gauge which is not shown 
in the drawing. Such a gauge can be used in this 
experiment since the gas pressure in the target 
cell is usually less than 7 cm of Hg and errors in 
this measurement amount to only _ several 
hundredths of a millimeter range. 

The disintegration particles to be observed 
pass out through the mica foil at the side of the 
target cell. The direction of these particles is 
generally at 90° to the direction of the deuteron 
beam, but no collimation is introduced at this 
point. After passing through the foil the dis- 
integration particles traverse the range cell in 
which the pressure of dry air can be varied by 
means of an external system. This pressure is 
measured by a mercury manometer. The tem- 
peratures of the target cell, range cell, and 
ionization chamber are measured with three 
thermometers in copper tubes soldered to the 
respective parts. The particles enter the ioniza- 
tion chamber through a collimating head which 
supports an aluminum foil, which is one elec- 
trode of the parallel plane ionization chamber. 


Foil tube 


Range Cell 


lonzation ¥ 
Chamber 


Fic. 1. Sketch of the assembly of target cell, range cell 
and ionization chamber used in measuring the ranges of the 
particles emitted in a direction 90° to the direction of the 
deuteron beam. The gaskets and screws which hold the 
components together are not shown. 
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The collecting electrode and its guard ring cap 
be moved parallel to its axis on threads of 1-mm 
pitch turned on the guard ring cyclinder. The 
physical depth of the ionization chamber jg 
measured by moving the collecting electrode 
forward until it makes electrical contact with the 
collimating head and foil, then the guard ring 
cylinder is screwed back. The number of revoly. 
tions of this cylinder determines the physical 
depth of the ionization chamber. In some cases 
the ionization chamber was modified, but these 
modifications will be mentioned in connection 
with the particular experiment. 

The collimating head on the ionization cham. 
ber allows a maximum angle of 5.4° with the 
direction of the axis of the range cell and ioniza- 
tion chamber. It is necessary to have such a 
collimator because for some of the short range 
measurements the physical separation of the 
two foils bounding the range cell was only 1.4 
cm and the obliquity correction in this case was 
0.5 percent and less in all other cases. All three 
gridworks were drilled from the same pattern 
and have a ratio of open to total area of 0.54. 

The collecting electrode of the ionization 
chamber is connected directly to the grid of the 
first amplifier tube of the 6-stage linear pulse 
amplifier described previously.‘ The output of 
this amplifier is connected to the grid of a 
discriminating thyratron which operated a 
mechanical counter. This grid is biased so that 
the output voltage pulse must be above a 
certain value before the thyratron will flash. 
With no signal on the grid, the thyratron will 
flash at —15 volts grid bias, so for any negative 
bias, V, greater than 15 volts, the pulse height 
required to flash the tube is V—15 volts. It will 
be shown in a later section that the method of 
range measurement used in these experiments 
makes it necessary that the amplifier be linear 
in its response to voltage pulses. For this reason 
it was necessary to make frequent measurements 
of the response of the amplifier. This was done 
with the thyratron oscillator described in 
reference 4, p. 24. 

It was necessary to use a current integrator 
because the deuteron beam was somewhat un- 
steady. This integrator is similar to the type 


M. G. Holloway and M. S. Livingston, Phys. Rev. 
18 (1938). 


| 
des 
has 
in 
tak 
late 
coll 
ing 
— 7} sup 
arg nece 
= | tu 
the 
S the 
| perc 
| 
Fo 
rd rubber 
thick 
quent 
| foils 
! The ; 
lot w 
polon 
| 
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described by Herb, Kerst and McKibben® and 
has been found to be very satisfactory, although 
in the summer months precautions must be 
taken to prevent surface leakage across insu- 
lators. In Fig. 1 it can be seen that the grounded 
collimating shield prevents the beam from strik- 
ing any part of the foil tube but the foil and 
supporting gridwork. Such a precaution is 
necessary, for if the beam shifts in position then 
a part of it will strike the solid edges of the 
supporting gridwork or the walls of the foil 
tube and no fraction .of this current can reach 
the target gas although it will be measured by 
the integrator. However, with the beam re- 
stricted to the gridwork it is certain that 54 
percent of the measured beam does get to the 
target gas. 

For some of the longer range groups it would 
have been impractical because of intensity 
considerations to obtain sufficient absorber 
thickness by removing the ionization chamber a 
large distance from the target cell and conse- 
quently mica absorbers were used. A set of mica 
foils was made from one large piece of mica. 
The air equivalents of two thin pieces of this 
lot were measured with a collimated beam of 
polonium alphas; these foils were placed in the 


5R. G. Herb, D. W. Kerst and J. L. McKibben, Phys. 
Rev. 51, 691 (1937). 


alpha-particle beam so that the average energy 
of the alphas in the mica corresponded to a resid- 
ual range of 3 cm of air. The areas and masses 
of the two pieces were measured and from these 
measurements the mass/cm?/cm air was found 
to be 1.44 mg/cm?/cm air and 1.41 mg/cm?/cm 
air; the average of these is 1.425 mg/cm?/cm air 
and this value was taken as the stopping power 
for that sample of mica for alphas of 3-cm residual 
range. From the mass/cm? of the thicker pieces 
their thickness in air centimeters could be ob- 
tained. However, some correction had to be 
made to the thicknesses so obtained because the 
factor 1.425 mg/cm?/cm air was obtained at a 
residual range of 3 cm and all of the thicker 
absorbers were used at much higher ranges. 
Livingston and Bethe® have given a table of the 
stopping power relative to air of aluminum for 
various velocities of alpha-particles and protons. 
The correction for mica was taken as half the 
correction for aluminum. For velocities above 
1.5X10*® cm/sec. alpha-particles and protons 
have almost the same range. While these cor- 
rections are somewhat uncertain they are of 
second order in the range measurement and 
errors in the corrections are not too serious. 
The corrections for the one platinum foil used 


* M. S. Livingston and H. A. Bethe, Rev.,Mod. Phys. 9, 
272 (1937). 
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4 Thyratron 
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Fic. 3. Diagram showing the 
range for ionization c 


in these experiments were made with the curve 
for gold. 


METHOD OF RANGE MEASUREMENT 


The determination of the “range” in an 
absorber of a group of heavy particles is in 
most cases the problem of measuring the range 
of the group in such a manner that the average 
initial energy of the particles in the group can be 
determined ; that is, a knowledge of the actual 
length of the path of a particle is useful only in 
that it is a means of obtaining the kinetic energy 
of the particle at some particular point. It is 
the purpose of this section to describe the method 
used in these experiments in obtaining the energy 
of the particle groups by means of range meas- 
urements. All ranges given in this paper are for 
dry air at 15°C and 760 mm of Hg. 

Each alpha-particle or proton as it moves 
through a gas produces ion pairs along its path. 
The number of ion pairs per cm path length 
dI(x)/dx is not constant for any one particle, 
but is a decided function of the velocity of the 
particle. In air the dependence of dI(x)/dx upon 
velocity is such that dJ(x)/dx has a maximum 
several millimeters from the end of its path. 


Residual Range-cm 


uantities involved in the correction to the 
mber depth and thyratron bias. 


This end of path is defined as that point along 
the path at which the particle has so little 
energy (<33 ev) that it cannot produce even 
one ion pair. In Fig. 2 is shown the ionization 
collected in parallel plate ionization chambers of 
different depths as a function of the distance of 
the face of the ionization chamber from the end 
of the path of an average alpha-particle. The 
differential specific ionization dI(x)/dx is given 
by the curve (0 mm) for zero chamber depth.’ 
The chamber depth for each curve is indicated 
by the number to the right of the curve. These 
curves for chambers of finite depth are obtained 
from the differential specific ionization curve by 
taking the integral 


dI(c) 
f 
z-8 dc 
where 6 is the chamber depth, and then normal- 
izing the ordinate so that the maximum of the 
curve for each chamber depth is 1.0.° 


7 Reference 4, p. 29. 
’From the work of I. Curie and F. Béhounek (J. de 


eer et rad. 7, 125 (1926) and Ann. de physique 3, 299 
1925)) it can be deduced that the ion pairs correspondi 
to the maximum of the differential specific ionization curve 
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In Fig. 3 is shown a curve giving the ioniza- 
tion collected in a S-mm chamber as a function 
of the distance of the face of the chamber from 
the end of the range of the alphas in the group 
having the mean range of the group. Since under 
rather wide operating conditions the voltage 
pulse on the grid of the first amplifier tube is 

portional to the ionization produced by the 
individual alpha-particle in the chamber, the 
ordinate of the curve in Fig. 3 will also represent 
the voltage pulse on the discriminating thyratron 
as a function of the ion chamber position. This 
jast is true only if the amplifier is linear in its 
response to voltage pulses. The drawing indi- 
cates a bias on the discriminating thyratron 
which would correspond to a pulse height 0.8 of 
the maximum average obtainable with the 5-mm 
chamber. With such a bias the chamber face 
could be at any position within 8 and the 
alpha-particle would still be counted by the 
discriminating thyratron. If the ion chamber 
face is anywhere outside of 8 the alpha would not 
be counted. If the range straggling is of the form 
exp [—x*/a?] it can be shown easily that the 
peak of the curve connecting number of counts 
with the position of the chamber face will be 
at the center of 8 or displaced by y+48 from 
the mean range. 8 and y are both functions of 4, 
the chamber depth, and of P, the pulse height 
below which the discriminating thyratron will 
not flash. 

For a very high bias corresponding to ioniza- 
tion pulses in the near neighborhood of 1.0, 
(6, P) is very small and (6, P) is very nearly 
equal to the distance of the position of maximum 
ionization from the end of range. In this case 
the peak of the counting (high bias number- 
distance) curve will occur at the same position 
as that of the maximum of the ionization curve 
for that particular chamber depth. 

It is very easy to realize experimentally these 
last conditions. The procedure is to take a 
number-distance curve and use a bias known 
from experience to be slightly higher than the 
mean pulse height at the maximum of the 


is 6760 ion pairs/mm. With this value, the maxima of the 
curves for finite chamber depth are: 2 mm, 13,480; 4 mm, 
26,400; 5 mm, 32,450; 6 mm, 38,400 ion pairs/mm. The 
numbers give the number of ion pairs produced in the 
chamber when the face of the chamber is set at the posi- 
tion of the maximum of the curve. 
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ionization curve. The peak of this number- 
distance curve will be very near the position 
of the peak of the ionization curve; the accuracy 
of this being limited by the rather poor statistics 
coming from the few counts at high bias. The ion 
chamber is then placed at the position of maxi- 
mum counting rate and a number bias curve 
taken (see Fig. 8, reference 4). At this position 
of the chamber the bias for half-maximum count- 
ing rate is then the bias corresponding to the 
mean of the ionization pulses arising from parti- 
cles having their maximum ionization in the 
ion chamber. The bias of the discriminating 
thyratron is set at this bias which corresponds 
to ionization pulses of height 1.0 in Fig. 3. 
One must remember that the curve in this 
figure represents the ionization pulse such that 
half the particles in the group produce larger 
and half smaller pulses than that given by the 
ordinate of the curve. The number distance 
curve is then repeated with the bias correspond- 
ing to ionization pulse height 1.0. The true 
mean range of the group of particles is the sum 
of the distance between source and the position 
of the peak in the counting rate curve and of 
y(6, 1.0) for the particular chamber depth used. 
If the bias is for a pulse height less than 1.0 
then the correction is not y(é, 1.0) but (6, P) 
+ 38(6, P) appropriate for the particular bias. 
The reason for using a bias for pulse height as 
near 1.0 as possible is that this bias gives a 
larger ratio of maximum to half-width of the 


ALPHA-PARTICLES 


5 -mm 


wit 


Fic. 4. The correction to the range of alpha-particles for 
thyratron bias and ionization chamber depth as a function 
of the chamber depth for different pulse heights (different 
thyratron biases). 
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Residual Range ~cm 


Fic. 5. The average differential specific ionization 
quot by a proton. 


counting rate vs. position curve than any other 
bias. 

From the set of curves in Fig. 2 the correction 
(y+48) for different depths 6 and pulse height P 
can be obtained. Figure 4 gives this correction 
as a function of 6 for different P’s for alpha- 
particles. Again it should be mentioned that 
for most experiments P21.0 and the range 
correction for bias and chamber depth is 38(6, 1.0) 
+~/(6, 1.0); for P<1.0, the correction is $8(6, P) 
+7(6, P). 

In order to check this method of range meas- 
urement the mean range of polonium alphas was 
measured with a 2-mm chamber and with a 
4-mm chamber. The 2-mm chamber measure- 
ment gave 3.84; cm and the measurement with 
the 4-mm chamber gave 3.85; and 3.85» in two 
measurements. Since we have taken 3.842 +0.006 
cm as the standard range of the polonium 
alpha-group it would appear that this method of 
measurement gives the correct range to an 
accuracy limited by errors in measurement of 
the dimensions of the apparatus. The range and 
energy of the polonium alpha-group gives a 
point on the range-energy relation and thus, if 
the mean range of some new group of alpha- 
particles were measured by the method described 
above and found to be 3.84 cm one would be 
certain that the mean energy of the group is 5.29 
Mev. Similar points betweeen 2 Mev® and 10.5 
Mev"® based on measurements of the range and 
aS: Mano, Ann. de physique 1, 407 (1934); G. H. 

s, Proc. Roy. Soc. 114, 341 (1927). 
. B. Lewis and B. V. Bowden, Proc. Roy. Soc. 145, 


235 (1934) ; S. Rosenblum and G. Dupouy, J. = phys. et 
rad. 4, 262 ‘(1933). 


energy of alphas from naturally radioactiye 
sources make the method reliable in this region, 
The theoretical extension of the range ene 
relation to energies higher than 10 Mey j, 
straightforward and is quite accurate. The curye 
in the region below 2 Mev is obtained by inte. 
grating the differential specific ionization curve 
and fitting the integrated curve to the experi. 
mental curve at about 4 Mev; this may be in 
error because of a possible variation in energy 
loss per ion pair in the low energy region. 

The measurement of the range of protons is a 
problem similar to that for alpha-particles, 
However, the differential specific ionization for 
an average proton of a group is not as well 
known as that of the alpha-particle. Parkinson, 
Herb, Bellamy and Hudson" in their work on 
the range of protons in air used a very shallow 
(0.04 cm) ionization chamber and obtained 
curves of ion current vs. energy (for a fixed air 
path) showing the familiar specific ionization 
curve for heavy particles. These curves of 
specific ionization vs. energy can be transformed 
into curves of specific ionization vs. range. In 
particular, curve VI, Fig. 5, reference 11, was so 
transformed. The transformed curve is a ‘‘Bragg” 
curve for a proton group having in this case a 
very small straggling. The effect of the finite 
chamber depth and range straggling is to put a 
small tail at the foot of the curve and to decrease 
the slope of the face of the specific ionization 
curve (this is discussed in reference 4). Initially 
it was assumed that the true differential specific 
ionization curve was the same as curve VI, 
Fig. 5, reference 11, except that the linear face 
of the specific ionization curves was extended 
directly to the axis, i.e., no tail. To this curve 
was then applied a straggling of 0.5 mm (more 
than would be expected from the Wisconsin 
high voltage generator) and a new curve com- 
puted. There was no appreciable change in the 
curve except that a tail was added to the foot 
and the curve displaced by about 1 mm. The 
effect of the finite chamber depth of 0.4 mm 
was also to add a tail and decrease the slope 
by about 15 percent. The differential specific 
ionization curve of the average proton of the 
group was then taken as curve VI, Fig. 5, 


"1 —D. B. Parkinson, R. G. Herb, J. C. Bellamy and C. M. 
Hudson, Phys. Rev. 52, 75 (1937). 
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reference 11, with the linear face extrapolated 
to the axis and a 15 percent increase in slope of 
the face. This curve is shown in Fig. 5. In 
measuring ranges the same procedure as for 
alpha-particles is followed and the correction 
48(5, P)+~7(6, P) for protons is given in Fig. 6. 
Because the curves for P=0.9, 0.8 and 0.7 are so 
near the curve for P=1.0, only the correction 
for P=1.0 is given. In as much as the proton 

uces only about a third the maximum ioniza- 
tion per millimeter as an alpha-particle it is 
very easy to distinguish between the two 
types of particles by the height of the voltage 
pulse which they produce in the output of the 
amplifier.” 

The range-energy relation for protons was 
taken from the work of Parkinson, Herb, 
Bellamy and Hudson with some corrections to 
their data to get mean ranges and mean energies. 
Above 1.95 Mev the relation was obtained by 
transforming the- alpha-particle range-energy 
relation by means of the relations 


Ru(v) (Mu/ Ma) (Za/Zu)*?Ra(v) 


and 
En(v) = (Mu/M,)E,(v). 


The constant c takes care of the difference in 
the capture and loss effect between proton and 
alpha-particles, and is equal to 0.30.06 cm. 


MEASUREMENT OF THE DEUTERON ENERGY 


It was necessary to determine as accurately as 
possible the energy of the deuterons from the 
cyclotron, since the energy of the bombarding 
particle enters into the calculation of the dis- 
integration energy of a nuclear reaction. Theo- 
retically the energy can be calculated from the 
frequency of the oscillating “dee” voltage and 
the geometry of the cyclotron chamber. However, 
uncertainties in the determination of the radius 


“For alpha-particles and protons of energies greater 
than 1 Mev and 0.1 Mev, respectively, the loss and gain 
of charge is negligible, and for alphas and protons of the 
same velocity the ratio of the differential specific ioniza- 
tion of the alpha-particle to that of the proton should be in 
the ratio of the square of their charges, namely 4/1. Com- 
parison of the differential specific ionization curves of the 

-particle and the proton in the region of 3-cm residual 
range gives an absolute ordinate for the proton curve such 
that the maximum corresponds to 2320 ion pairs/mm. If a 
very deep chamber is used it is necessary to extend the differ- 
ential specific ionization curve to — residual ranges. 
This can be done by using Eq. (749), reference 6. 
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of curvature of the last half-circular path make 
this method unsatisfactory for accurate work. 
A more direct method is to measure the range of 
the particles, but because of the large beam 
intensity a pulse amplifier could not be used to 
count the individual particles and as a result 
the correction for chamber depth could not be 
made accurately. This inaccuracy in the determi- 
nation of the range of the very short range 
deuterons was too large to allow the use of this 
method. 

A third method, which was the one adopted in 
the present work, is to use some nuclear reaction. 
The one chosen here was the reaction H?(d, p)H*. 
If the Q value of this reaction and the energy of 
the emitted proton (say at 90°) is known the 
energy of the bombarding particle can be com- 
puted. However, by measuring the energy of the 
protons emitted in the forward direction and at 
90° to the deuteron beam it is possible to de- 
termine both the Q value and the energy of the 
incident deuterons. Once the Q value is known 
it is necessary to measure the proton range only 
in the 90° direction in order to obtain the 
deuteron energy. 

The experimental arrangement for measuring 
the protons in the 90° direction was the same as 
that previously described with a 19.2-cm range 
cell. For the measurement at 0° a 33-cm range 
cell and a target cell with a foil on the forward 
end was used. The depth of the target cell was 
1.1 cm. Since the energy of the deuterons 
incident on the target gas was about 1 Mev and 


30 
20 
= 
= 
+10 — 
WA PROTONS 
10 20 30 
6-mm 


Fic. 6. The correction to the range of protons for ioniza- 
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the air equivalent of the gas was about 0.1 cm 
there was left about 1.4 cm of the deuteron 
range to be stopped in the exit foil. At this 
bombarding energy there is a group of protons 
from aluminum of about 35-cm range which 
might have interfered with the measurement of 
the group from the deuterium. In order to 
eliminate this the last foil was made up of a 
platinum foil of 1.1 cm air equivalent and an 
aluminum foil of 0.46 cm air equivalent. This 
last foil was not entirely made of platinum 
because the correction of the stopping power of 
platinum for different velocities is not very 
well known, and it was desirable to keep this 
correction as small as possible. The excitation 
function for the aluminum reaction becomes very 
small for 0.4-Mev deuterons and no protons from 
the small thickness of aluminum were observed 
in the region of 33-cm range. 

Equations (1) and (2) relate the Q value and 
the deuteron energy £;, with the energies of 
the protons emitted in the forward and 90° 


directions." 
(M2+Ms3)? 
Me EY 


M; (M3— M,)(M3+ M2) 
M; 


—E,*° | (2) 
Ey J 


(1) 


The following notation is used: The subscript 
0, 1, 2 and 3 refer to the initial nucleus, incident 
particle, emitted particle and residual nucleus, 
respectively. The superscript 0 will refer to the 
measurement in the 0° direction and 90 to the 


18 From conservation of mass and energy Q= E:+£E;— Ei 
and from the conservation of momentum 


M;E;= MiE:+ cos 6, 
where @ is the angle between the direction of the incident 
and emitted particles. Measurements of the energy of the 
protons emitted in the forward direction and in the 90° 
direction allow the use of two independent equations 
involving the unknown quantities Q and E,. 


M:+M; 
(a) E, at 90°, 


(M2E,°)' Ei at 0°. 
M; 


(b) Q=£,°+ 


From these two equations can be obtained the Eqs. (1) 
and (2) above. 


90° direction. The cyclotron was not shut down 
during the change of target cells so that condi. 
tions in the cyclotron would be the same for the 
two measurements. The mean range of the 
protons in the forward direction was 30.63 cm 
corresponding to an energy of 4.71; Mev; in the 
90° direction the range and energy was 15.94 om 
and 3.22 Mev, respectively. Using these values 
of E,° and E,* in Eqs. (1) and (2) one is led to 
the value 0.95 Mev for the energy of the dey. 
terons in the center of the gas in the target cel] 
and 3.98 Mev for the Q value for the reaction 
H*(d, p)H*.% Since the air equivalent of the 
aluminum foil through which the deuterons 
passed into the target cell was 0.80 cm the 
energy of the deuterons incident on the foil was 
1.27 Mev. From this energy value and the 
frequency of the cyclotron oscillator the radiys 
of the last circle in the cyclotron can be computed 
and used for future determinations of the 
deuteron energy. However, this must be checked 
frequently because adjustments of the cyclotron 
can very easily change the geometry enough to 
require a new calibration of the beam energy. 
In order to see how errors in E,° and FE,” 
affect E, and Q we take the differentials of 


Eqs. (1) and (2). 
4MiM, Ls 


2 


6(E2° — 


1 


M:+Ms; (Ms—M,) 
50 =————6E. —_— (4) 
M; M; 


The first item in (3) involves inaccuracies in the 
difference (E,°—E,®) and the second term 
inaccuracies in E,° alone, and since the errors in 
the two terms are approximately independent of 
each other the two terms must be added as the 
square root of the sum of their squares. 

In the case of the measurement of the beam 
energy by the method of observing the emitted 
protons in two directions from the disintegration 


“4M. L. E. Oliphant, A. R. Kempton and Lord Ruther- 
ford, Proc. Roy. Soc. 149, 406 (1935), give 3.98+0.02 Mev 
and F. E. Myers, R. D. Huntoon, C. G. Shull and C. M. 
Crenshaw, Phys. Rev. 56, 1104 (1939), give 3.93+0.10 


Mev. 
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of deuterium, the peak of the counting curve 
could be determined to about 2 mm of Hg 

ure and the temperature to about 0.5°C. 
The estimated error in the measurement of the 
length of the 19.2 cm long range cell was about 
0.01 cm. These errors give together a total error 
in the air equivalent of the range cell of 0.06 cm 
of air. There are two foils whose air equivalents 
are known to about 0.01 cm. The combined 
error in the chamber depth measurement and 
the chamber depth corrections to the range is 
about 0.02 cm. These constitute the experimental 
errors in the range of the 15.94-cm protons (90°) 
and total about 0.065 cm error for the 19.2-cm 


range cell. The longer range cell of length 33.0 


cm has a somewhat larger error amounting in all 
to 0.08 cm in the 30.63-cm group observed at 0°. 
The difference (E2°—E,*°) does not involve the 
inaccuracy in one foil or the error in the chamber 
depth and its corrections, and consequently the 
error in this difference is 0.094 cm rather than 
0.103 cm as would be expected from adding the 
errors in the two ranges. The above errors are all 
the result of inaccuracies in measuring physical 
quantities in the course of the experiment. 

The proton range-energy relation of Parkinson, 
Herb, Bellamy and Hudson includes the region 
up to about 2 Mev and is accurate to about 
1 percent in range and energy. As was mentioned 
previously the range-energy relation for protons 
above 2.0 Mev is obtained by transforming the 
relation for alpha-particles. This includes the 
additive constant c which comes out to be 0.3 cm. 
By taking the probable errors in the alpha- 
particle curve and that of Parkinson, Herb, 
Bellamy and Hudson, the error in c is found to 
be 0.06 cm. Other errors involved in the absolute 
range scale for the protons are small compared 
to this error of 0.06 cm in c. However, for 
differences on the range energy relation for 
protons the accuracy is much better being about 
0.02 cm. 

The combination of the experimental errors 
and the errors in the range-energy relation gives 
the following errors for the ranges of the protons 
from the reaction H?(d, p)H*: In the forward 
direction, 0.10 cm and in the 90° direction, 
0.09 cm. These correspond to inaccuracies in the 
energy of E,° of 11 kev and in E,* of 10.8 kev, 
the slope of the range energy relation in these 
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regions being 115 kev/cm. The difference 
(E,°—E,®) does not involve the absolute error 
in the range-energy relation, but rather a smaller 
error of about 0.02 cm. The error in the energy 
difference is then about 5 kev. 

Experimentally it was not possible to deter- 
mine the direction of the deuteron beam to 
better than 2°, and as a consequence the energy 
values 4.71; Mev and 3.22 Mev, while being 
accurately the energy of the protons observed, 
may not be £,° and £,* since these two quanti- 
ties are defined as the proton energies in exactly 
the 0° and 90° directions. The inaccuracy in @ 
introduces an error in Q of 72 kev and in £, of 
52 kev. These errors are large compared to the 
errors arising from the range measurements, and 
are connected only with the determination of the 
beam energy and the Q value, and not with the 
method used in determining the energy of the 
emitted particle. The errors combine to give an 
error in Q of 74 kev and in £, of 54 kev. 


RESULTS 


The target gases used in these experiments 
contained a high percentage of the heavy 
isotope N®. Some of the earlier work was done 
with a sample containing 20 percent of the N™ 
and the later work with a 70-percent sample. 
The appreciable amount of the N™ in all of the 
samples made necessary the study of the normal 
nitrogen under deuteron bombardment so that 
the contribution of the rarer isotope could be 
determined. This work on normal nitrogen was 
essentially a repetition of the work of Cockcroft 
and Lewis, but since the measurements on N™ 
had to be made it was thought worth while to 
make them carefully. 

In Fig. 7 are shown two groups of alpha- 
particles resulting from the bombardment of 
tank nitrogen which contains only about 0.4 
percent of the N®. The reaction involved is 
N*(d, a)C”. In view of some work to be men- 
tioned later the small hump occurring at about 
5 cm range is probably due to the presence of 
the small amount of N® in the tank nitrogen. 
The long range group at about 11.67 cm was 
measured in two ways.” The first measurement 

The range scale for the curves of counting rate vs. 


distance from the target cell is not absolute, but may 
differ by as much as 1 mm from the absolute value. The 
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Fic. 7, Two groups of alpha-particles from the reaction N“(d, a)C. The target was tank nitrogen. 


was with a foil holder between the range cell and 
the ionization chamber and with a copper screen 
as the face of the ionization chamber. The end 
of the range cell was covered in this case with a 
mica foil of about 0.88 cm air equivalent. A 
thick mica foil of 4.50 cm air equivalent from 
the calibrated set mentioned before was in the 
foil holder. With corrections for the stopping 
power of the mica and for slight obliquity the 
range of the group was 11.67 cm air, the bom- 
barding energy being 1.01 Mev. The second 
measurement was with the 19.2-cm range cell 
used previously in measuring the range of the 
protons from the reaction H?(d, p)H*. There 
was no thick absorbing foil used in this measure- 
ment and the ionization chamber face was the 
Al foil over the collimating gridwork. The range 
found by this measurement was 11.97 cm at 
accurate ranges are obtained from the counting rate vs. 

essure curves. It would be much too laborious to trans- 
orm each experimental point to an absolute scale, for the 
corrections required to get an accurate range scale are the 
following: (@) The contribution of the target cell must be 
ometal for pressure and temperature, (b) the range cell 
must have a similar correction, (c) the foils on the ends of 
the range cell as well as any thick absorbers must be 
corrected for the ym 4 in stopping power with velocity, 


(d) the chamber depth must yay and the ap- 
propriate correction made, and (e) the obliquity factor 


computed and applied to the range. The ordinates of these 
curves do not represent absolute yields. 


1.02-Mev deuteron energy. The Q values for 
these two measurements are 13.30 Mev and 
13.48 Mev, respectively ; the average of these is 
13.39 Mev as compared with the mass Q value 
of 13.36 taken from the mass tables of Barkas.'* 

The group at 6.54 cm range (deuteron energy 
1.02 Mev) is also due to the reaction N"(d, a)C® 
in which the C® is left in an excited state. The 
average Q value as taken from 5 separate 
measurements was 9.02 Mev. The two different 
ionization chamber heads and range cells used 
for 11.7-cm group were used in getting the 5 
measurements of this group, although no thick 
mica absorber was involved. 

The 6.54-cm group is shown again in Fig. 8 
along with the short range group of alpha- 
particles at 3.47 cm, the deuteron energy being 
1.02 Mev. Since these curves were taken with 
only tank nitrogen in the target cell it is probable 
that this group is from the N™. Checks with 
the target cell evacuated and with CO; in the 
cell indicate that it is not from a casual contami- 
nant. The fact that its intensity was not en- 
hanced by the presence of N® in large quantities 
makes it likely that it does not arise from N®. 


This group at 3.47 cm appeared with the 


* W. H. Barkas, Phys. Rev. 55, 691 (1939). 
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N“—N® mixture, and with tank nitrogen in 5 
distinct observations which makes it very prob- 
able that it is due to the N™ reaction. The 
average Q value from two range measurements 
is 5.77 Mev. Cockcroft and Lewis did not find 
this group of alphas, but an inspection of their 
published curves shows no experimental points 
in the region of this small peak. 

In Fig. 9 are shown the alpha-particle groups 
obtained in the region 2- to 8-cm range when the 

t cell contained the gas mixture of 70 
percent N® and 30 percent N™. In addition to 
the two groups at 6.54 cm and 3.47 cm there is 
a strong group at 5.25 cm; these ranges are for a 
bombarding energy of 1.02 Mev. The average 
of three range measurements on this group with 
slightly different deuteron energies gives a Q 
value of 7.54 Mev. This must correspond to the 
transition to the ground state of C™ in the 
reaction N(d,a)C"%, for no alpha-groups of 
longer range were found which could be at- 
tributed to the presence of N® in the target. 
The Q value as computed from Barkas’ mass 
values is 7.51 Mev. 

The very intense distribution of alphas below 
about 3-cm range was observed by Cockcroft 
and Lewis. The shortest range which could be 
studied by the present apparatus was about 
2.1 cm, consequently any definite information 
concerning the short range distribution must be 
left for future measurements. 

In the first report of this work? three groups 
of protons were attributed to the reaction 
N¥(d, p)N®. The ranges of these groups were 


Ronge ~ cm 


Fic. 8. Alpha-particle grou from the reaction 
pal a)C", The group at about 6.5 cm is the same group 
shown in Fig. 7. The target was tank nitrogen. 
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given as 90.8 cm, 66.1 cm and 21.0 cm. The 
statistics in the early work were very poor in 
the case of the 90.8-cm group because of its 
small intensity compared to the neutron recoil 
background. Recently it was thought advisable 
to repeat the measurement of the 90-cm group 
with an ionization chamber of 1-cm physical 
depth filled with krypton at two atmospheres. 
The ionization produced in this chamber by a 
neutron recoil was small compared to the ioniza- 
tion from a proton which went completely 
through the chamber. With the high thyratron 
bias method very few of the neutron recoils were 
counted. The 90-cm group is shown in Fig. 10. 
The small dotted curve is the one taken with the 
5-mm ionization chamber filled with air. The 
thick mica absorbers used in this measurement 
amounted to 71.6 cm of air, including the 
correction for the change in stopping power with 
velocity. Rather than make corrections for the 
large chamber depth which would of necessity 


N“(d.«)C* 


Alphas min 
8 


Range ~ cm 


Fic. 9. Alpha-particle groups from a target of 70 percent 
5.2-cm range is from the 
a e other two groups were seen previously 

to be from the reaction N“(d, ace, 


require a knowledge of the stopping power of 
krypton the ionization chamber was calibrated 
by measuring the ranges of the protons for the 
reactions H?(d, p)H* and C"(d, p)C™. These two 
groups have a range of about 16 cm and have 
been carefully measured in this laboratory with 
the thin ionization chamber. The range of the 
group from the reaction N“(d, »)N® as measured 
with the thick ionization chamber is 90.5 cm. 
The older measurements with the thin ionization 
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Fic. 10. The two from the reaction 
»P)N”. 


chamber gave 92.5 cm when corrected for change 
in stopping power of mica with velocity although 
with a different deuteron energy. Because of the 
much better statistics it is thought that the 
90.5-cm range is the better value, and the Q 
value obtained with this range is 8.51 Mev as 
compared with the mass Q value of 8.57 Mev. 
In the early work with the thin ionization 
chamber the proton group at 66 cm was more 
intense than the 90-cm group. However, with 
the thick krypton chamber the 66-cm group 
could not be found at all! Various common 
possible contaminants were tried but none gave 
the 66-cm group. Finally because the range of 
the group was about that to be expected from 
boron a target of borax was tried, and the 66-cm 
group appeared with a very large intensity. The 
problem was then to find the source of the 
boron in the earlier experiments and why the 
group did not show up when the target had been 
evacuated or filled with carbon dioxide but was 
present when tank nitrogen or the N*—N*® 
mixture was used. It was remembered that a 
Bourdon type vacuum gauge had been used on 
the target cell for testing the cell for leaks and 
that this gauge had been discarded because of a 
leak in its flattened tube. Further investigation 
revealed that this gauge had been used several 
years previously in the filling of boron trifluoride 
ionization chambers, and had been laid aside 
because of the leak. This explained the source of 
the boron in the early experiments, for pre- 
sumably the cavity into which the gauge had 
been screwed had become saturated with the gas. 
- The time required for taking the points for one 
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of the long range proton groups was of the order 
of several hours and the N*—N® mixture or 
the tank nitrogen might remain in the target cel} 
for as much as a day. This would allow ample 
time for some of the boron trifluoride to 

into the bombarding cavity. However, when the 
tests were made with the cell evacuated or with 
carbon dioxide very little time was required, 
since the method of checking the groups was to 
set the ionization chamber at the position of the 
peak of the curve and measure the counting rate, 
If the counting rate was the same as the back. 
ground it was assumed that the group did not 
come from the cell walls or the carbon dioxide. 
In the time required for the check it is probable 
that very little of the boron trifluoride got into 
the target cell. 

The strong group of protons at 21.02 om 
(deuteron energy of 1.01 Mev) is from the 
reaction N'(d, p)N®. The Q value is 3.15 Mey 
which indicates an excitation level in N® of 
5.42 Mev. 

A group. of protons might be expected to 
result from the reaction N"(d, p)N**. The resid. 
ual nucleus N’* is known to be radioactive with 
a half-life of 8 sec. Naidu and Siday"’ measured 
about 240 tracks in the beta-ray spectrum of N' 
and obtained a visually extrapolated end point 
at little over 6 Mev. Fowler, Delsasso and 
Lauritsen'* found a visual end point of about 
6 Mev and a K-U limit of 6.5 Mev. If 6.5 Mey 
is taken as the energy of the beta-particles given 
off from in the disintegration N'*-O'%+¢- 
the mass of N"* has the lower limit of 16.007 MU 
and may be heavier by the energy of any 
gamma-rays involved in the disintegration. In 
the present experiment a thin-walled Geiger 
counter tube was set up very near to the foil 
window of the target cell. A chronograph was 
used to record every sixteenth count of the 
Geiger counter. With the target cell filled with 
the N'*— N*® mixture two half-life measurements 
were taken of the activity produced in the gas. 


The first measurement showed an 8.1-sec. half- 


life superimposed on a 5.6-min. half-life. The 
second measurement gave a 10-sec. half-life on 


17R. Naidu and R. E. Siday, Proc. Phys. Soc. London 


48, 332 (1936). 
18 W. A. Fowler, L. A. Delsasso and C. C. Lauritsen, 


Phys. Rev. 49, 561 (1936). 
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a’ 7-min. half-life. It was determined that the 
longer periods came from the cyclotron chamber 
and not from the target cell and presumably 
this longer half-life radioactivity came from the 
N® formed in the cyclotron chamber from 
carbon contamination. The somewhat small half- 
life value for the N® is due probably to the 
diffusion of the radioactive nitrogen from the 
walls into the pumping system. With the target 
cell filled with tank nitrogen there was no short 
life component found at all, but rather one of 
6.7 min. half-life, which again was found to 
come from the chamber. Since N’* has a half-life 
of 8 sec.!* the short-lived activity detected as 
coming from the bombardment of the N® with 
deuterons is probably from N'* formed in the 
reaction N*(d, If this is assumed to be 
the case then the fact that it is possible to 
produce the above reaction with 1-Mev deu- 
terons sets an upper limit of 16.0126 MU for the 
mass of N'*. This value is for zero energy of the 
proton which is given off in the deuteron reaction. 
If the value 16.007 is used for the mass of N"*, 
the Q value for the reaction N"*(d, p)N"* is 4.19 
Mev and with 1-Mev deuterons the energy of 
the emitted proton should be 4.75 Mev or a 
range of 31 cm of air. The entire region from 2- 
cm to 100-cm range was searched for protons 
with the N“*—N*® gas mixture in the target 
cell, and there were found 3 groups of protons at 
21 cm, 66 cm and 90 cm; these groups were also 
found with the tank nitrogen and their source is 
known. No new proton group was found. Below 
12-cm range it was very difficult to detect 
protons because of the large numbers of alpha- 
particles in that region. Figure 11 shows a 
typical curve taken with a low bias for detecting 
protons. Since the weak group of alphas at 3.5 
cm does not show on this curve, it is evident 
that a proton group would have to be rather 
strong to be detected over the strong background 
of alpha-particles. With only the results of the 
present experiment the matter of the protons 
from the reaction N(d, p)N** must be left in a 
rather inconclusive state. 

In Table I the results of all these measurements 
are summarized. £; is the energy of the deuterons 


in the center of the gas target; E, is the energy 


J. J. Livingood and G. T. Sea , Rev. Mod. Phys. 
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of the emitted particle whose range in air (at 
15°C, 760 mm Hg) is given in the third column. 
Qexp is the Q value as computed from the 
present measurements and in the next two 
columns are given the Q values from Cockcroft 
and Lewis and from the mass tables of Barkas. 
The absolute values of the cross sections are 
probably accurate to a factor of 5, but the 
relative values are good to 50 percent. 

The agreement between the Q values found 
in the present work and those of Cockcroft and 
Lewis is good except in the case of the alpha- 
group corresponding to the transition to the 
ground state of C”; the difference is 0.17 Mev, 
which is within the experimental errors given. 
The agreement is excellent in every case with the 
corrected values deduced from Cockcroft and 
Lewis’ experiments by Livingston and Bethe.”° 
The experimental Q values reported here and 
the Q values obtained from mass values agree to 
at least 60 kev. 

Consideration of the energy levels shown in 
the table would lead one to expect that the 
bombardment of nitrogen with deuterons would 
produce gamma-rays having energies 4.37, 7.62, 
5.42 and 3.25 Mev; the last gamma would come 
from a transition between the 7.62- and 4.37- 
Mev levels in Gaerttner and Pardue” have 
observed five gamma-ray lines resulting from 
the bombardment of nitrogen with deuterons; 
the energies of these lines were given as 2, 4, 
5.340.4, 7.2+0.4 and 11 Mev. The 4- and 
7.2-Mev lines could correspond to the 4.37- and 
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Fic. 11. Curve taken in the search for protons of range 
than 8 cm. 


2° Reference 6, pp. 279, 321 and 327. 


1 E. R. Gaerttner and Louis A. Pardue, Phys. Rev. 57, 
386 (1940). 
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TABLE I. Data on the bombardment of nitrogen by deuterons. 


RANGE E: 


REACTION AIR CM MEV MEV 


Cross 
SECTION 


MEV cm? 


N¥(d, 11.67+0.11 10.60+0.07 


4.98 +0.06 


3.50+0.07 
4.96+0.06 


3.47+0.07 


11.97+0.07 10.76+0.06 13.394%0.08 13.22+0.1 


9.02+0.07 


5.77+0.07 


13.36+0.16 1x10" 


8.93 +0.03 3X 10-7 


6.37 +0.06 


5.10+0.07 
6.48 +0.06 


5.25 +0.07 


7.54+0.07 


7.51+0.28 1x10" 


8.75+0.1 
(8.86) 


3.78 +0.02 


N¥(d,p)N® 0.95 90.5 +1.6 
(1.01) (92.5) 


1.01 21.02+0.11 


8.51+0.1 
(8.57) 


3.15+0.05 


8.53+0.1 8.57+0.23, 10-8 


3.25+0.07 3x 1077 


7.62-Mev levels in C” and the 5.3-Mev line with 
the 5.42 level in N”. The relative intensities of 
the particle transitions to the 4.37-Mev level in 
C® and the 5.42-Mev level in N"“ agree approxi- 
mately with the relative intensities of the 4- 
and 5.3-Mev gamma-ray lines. However, it 
should be pointed out that the ratio of the 
numbers of the particle transitions to the 4.37- 
and 7.2-Mev levels in C” is 10 whereas the ratio 
of the numbers of the 4- and 7.2-Mev gamma-ray 
quanta as given by Gaerttner and Pardue is 0.45. 
In order to explain this one might assume the 
existence of another strong group of alphas 
corresponding to a still higher excited state of C”. 
Such a group would of course fall into the very 
short range region and could not be detected in 
the present measurements. The corresponding 
excited state of C” would be unstable against 
alpha-emission, but it is still easily conceivable 
that such a state could not actually emit alpha- 
particles because of selection rules. If one as- 
sumes that such a state exists at about 9 Mev 
it is easy to account for the intensity of the 
7.2-Mev line as well as the existence of the 2-Mev 


line which could originate in the transition from 
the 9 Mev to the 7.62-Mev level. This picture 
is given some support by the fact that the 
intensities of the 2-Mev and the 7.62-Mey 
gamma-ray lines are about the same. Crane 
Delsasso, Fowler and Lauritsen,” in addition to 
the first five lines reported by Gaerttner and 
Pardue, observed a gamma-ray line of 3.1-Mev 
energy. This could result from a transition 
between the 7.62- and the 4.37-Mev levels in C®. 
Crane, Halpern and Oleson* report gamma-rays 
of energies 8.2, 6.6, 5.1, 4.1, and 2.5 Mev as 
coming from deuteron bombardment of nitrogen. 
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The Specific Charge of the Positron 


A. H. Spees* anp C. T. Zaunt 
Depariment of Physics, University of Michigan, Ann Arbor, Michigan 
(Received September 20, 1940) 


Previous measurements of the specific charge of beta-particles have been supplemented to 
include data on positrons by making a direct comparison between electrons and positrons. The 
positrons and electrons were obtained from artificially radioactive copper. Comparisons were 
made for beta-particles having approximately the velocity 8=0.72. The results indicate that 
the ratio ¢/mp for positrons is the same as for electrons within an error of approximately two 


percent. 


INTRODUCTION 


COMPARISON of the specific charges of 

negative and positive beta-rays by a de- 
flection method was made by J. Thibaud! 
shortly after the discovery of the positron. 
Employing the trochoid method of collecting 
charged particles in a magnetic field, he observed 
the deflections of positrons and electrons in a 
combination of magnetic and electric fields and 
concluded that the specific charges did not differ 
by more than fifteen percent. Another experi- 
ment reported by E. Rupp? presented results 
indicating equality to within five percent. This 
work together with other articles on positrons 
was, however, retracted.*? From cloud-chamber 
data Chadwick and his co-workers‘ obtained a 
comparison not of the specific charges but of the 
masses of positron and electron and these were 
found to be the same to within ten percent. 


EXPERIMENTAL PROCEDURE 


In view of the rather low accuracy of these 
previous determinations it was decided worth 
while to attempt measurements yielding more 
precise results. The method and apparatus used 
in the previous work on the specific charge of 
beta-rays® proved sufficient for this purpose. 
For convenience the method is here briefly 
reviewed. When the Lorentz mass formula is 
introduced into the formula for the curvature k 


* Now at the Detroit Institute of Technology. 

t Now with the American Viscose Corporation. 

1 {: Thibaud, Phys. Rev. 45, 781 (1934). 

*E. Rupp, Zeits. f. Physik 92, 485 (1934). 

*E. Rupp, Zeits. f. Physik 95, 801 (1935). 

‘Chadwick, Blackett and Occhialini, Proc. Roy. Soc. 
144, 235 (1934). 

°C. T. Zahn and A. H. Spees, Phys. Rev. 53, 357 (1938) ; 
53, 366 (1938). In the present article these references are 


referred to as papers I and II. 


of the circular orbit that a charged particle e 
takes in a uniform magnetic field H, one obtains: 


e/mo= 


Upon the determination of k and 8 for a fixed 
magnetic field H this formula yields a value of 
the specific charge of the particle. In the experi- 
mental arrangement employed, a diagram of 
which is given in Fig. 1, the curvature & is 
determined by geometric conditions in the mag- 
netic field before the rays enter the velocity 
selector represented by the electric condenser. 
These geometric conditions include the relative 
position of the beta-ray source and the condenser, 
the size of the source, and the condenser spacing. 
Because of the latter two features the rays that 
are permitted to enter the condenser will have, 
not just a single curvature koo, but curvatures 
in a spread dk about a curvature ko. This 
curvature is that of a central ray which leaves 
the source at the center and enters the condenser 
in the central plane, in a direction parallel to 
the condenser plates. 
According to the above formula for e/mpo, a 
definite value of k=koo as thus determined for a 
fixed field intensity H and the assumption that 
e/my is a constant, determine a value of B= Boo. 
Therefore since at all times only particles having 
approximately k=koo enter the condenser space 
(the region of crossed electric and magnetic 


Fic. 1. Diagram of experimental arrangement. 
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Positrons 
Execrrons 


BC OUNTS BER 


Vours — 
Fic. 2. Transmission peaks for positrons and electrons from 
a Cu source. 


fields) it follows that the particles transmitted 
by the condenser will have approximately a 
velocity 8 =8o0. But this transmission, it is to be 
remembered, may occur only for a compensating 
setting of the electric field EZ, namely when 
Boo=E/H. As a consequence of such action a 
plot of the counts of particles emerging from the 
condenser space against various settings of the 
electric field (magnetic field constant) should 
show a peak indicating maximum transmission 
at a certain value E=Eoo. An intensity peak 
obtained in this manner must, due to the finite 
resolution of an actual apparatus, have a definite 
spread. The quantity E determined experi- 
mentally in the fashion just outlined gives a 
measure of the velocity Bo which is to be 
inserted in the formula for e/mo. Particles having 
different specific charges will obviously yield 
peaks that have different values of Eoo. 

The procedure for the problem of comparing 
the specific charges of electron and positron was 
based on the use of a beta-ray source which 
emits both kinds of particles. Copper when 
bombarded with deuterons has just this property. 
The unstable isotope Cu“ that is produced 


decays with the emission of both positrons and 
electrons, the half-life being 12.8 hours. For the 
present experiment this branching reaction wag 
convenient because measurements on both kinds 
of particles could be made without disturbing 
the geometry of the apparatus. In order to 
change from an electron to a positron source all 
that was required was a reversal of the electric 
and magnetic fields. The procedure, therefore, 
was to obtain a transmission peak for the 
electrons and then after reversing the fields to 
obtain another for the positrons. The relative 
positions of two such peaks should give a direct 
comparison of the specific charges that is practi- 
cally independent of any geometrical features. 

The particular beta-ray sources used in the 
experiment were obtained from copper soldered 
to the end of a water-cooled probe that was 
placed directly between the duants of the 
Michigan cyclotron. With this arrangement a 
bombardment of 8 hours with a beam of 9-My 
deuterons gave a source of several millicuries 
strength. In the apparatus the source had to be 
shielded with lead in order to absorb some of 
the 0.5-Mv positron annihilation radiation that 
would have produced a very high background 
count in the particle counter. 


OBSERVATIONS 


The experimental peaks for two sets of 
measurements are given in Figs. 2 and 3. The 
sets differ in that the positions of the sources 
were slightly different and hence correspond to 
slightly different values of koo and Boo. The 
measured points of the present peaks are slightly 
more scattered than those of paper II obtained 
from a RaE source. This resulted from the 
necessity of using shorter counting periods due 
to the time limitation set by the shorter half-life 
of 12.8 hours and the additional requirement of 
obtaining two peaks instead of just one. As for 
the spreads, the present peaks are seen to have 
a ten percent half-width as compared with the 
four percent of the Ra E peaks. This change is 
to be expected because of an increase in source 
size. It was necessary to use a 1-cm square 
copper source because of intensity considerations; 
whereas the Ra E source was a 1-cm length of 
0.5-mm wire. In view of this relatively large 
increase in source size it is indeed surprising at 
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first sight that the peaks turned out to be as 
narrow as they did. According to resolution and 
intensity calculations for a finite source of 1 cm 
the expected half-width for over-all resolution 
js about 6.5 percent which is to be compared 
with the previous result of 2 percent for a line 
source. As previously, the greater width of the 
observed peaks may be attributed to scattering. 
It might have been possible to increase the 
resolving power in the present case by using the 
defining slits mentioned in paper II. These, 
however, had to be removed to make way for 
the lead shielding and so were not available for 
the measurements with positrons. 

In the transmission curves of paper II there 
occurred some pronounced subsidiary peaks to 
each side of the main peaks. An examination of 
the present peaks gives a slight indication of the 


game features. The principal reason they are 


not more pronounced is that the time factor did 
not permit drawing them out. In some prelimi- 
nary work with a small condenser it was, 
however, definitely established that the side 

also occur for an extended source. Thus 
since the side peaks arise from sources differing 
in size and shape, as do the Cu and RaE 
sources, they can hardly be ascribed to source 
conditions. Furthermore an extended and more 
thorough resolution calculation precludes the 
explanation that the side peaks may arise from 
some second-order effect inherent in the char- 
acteristics of the instrument as a whole. Their 
origin, therefore, is probably the one mentioned 
in previous discussion, namely, scattering from 
the condenser plates. 

Since the present investigation is concerned 
with relative values of the specific charge it was 
considered unnecessary to attempt an accurate 
determination of absolute values. It seems worth 
while, however, to discuss briefly the results in 
this connection. For the extended copper source 
used it is rather difficult to choose a proper 
central curvature ko in that one has to allow 
for a non-uniform distribution of radioactivity. 
Hence one cannot expect the accuracy that was 
attained with the small RaE source. The 
measured values of ko for the two source 
positions used for Figs. 2 and 3 were corrected 
for stray electric field and non-uniform magnetic 
field, in the manner outlined in paper II. The 
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corrected values were then used with values of 
the velocity as given by the graphs, namely, 
Boo=0.71 for Fig. 2 and Bo.=0.73 for Fig. 3, 
to obtain values of e/mo. In the former case the 
obtained value is lower than the accepted 
precision value by about 7 percent and in the 
latter case by about 10 percent. In the previous 
investigation with the Ra E source the observed 
peak position differed from the expected position 
by about 2.5 percent which means a discrepancy 
with the accepted value of e/my of about five 
percent. 


CONCLUSION 


In conclusion, the answer to the main problem 
of how the specific charges of the electron and 
the positron compare may be read from Figs. 2 
and 3. The peaks for electrons and positrons 
have positions on the voltage scale that corre- 
spond to voltages differing by no more than 
one percent. To find the corresponding difference 
in terms of specific charge one uses the following 
formula derived from the above formula for e/mo: 


A(e/mo) /(e/mo) =[1/(1 —Boo*) V/ Vo]. 
Thus the percentage difference in e/mp is larger 


Positrons 
ELectrons 


BCouwrs Yin. 
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Fic. 3. Transmission peaks for positrons and electrons from 


a source position slightly different from that of Fig. 2. 
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than that in the voltages by a factor containing 
the compensated central velocity Boo. For the 
actual values of Boo used in the experiment this 
factor is approximately 2.0. From the percentage 
difference in voltage of about one percent one 
would therefore conclude that the specific charge 
of the positron is the same as that of the electron 
to within approximately two percent. 
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If it be assumed that fission of heavy nuclei takes place in competition with the escape of a 
neutron from the highly excited compound system, we should expect that, for sufficiently high 


excitation of the system, fission of the residual nucleus left after neutron escape may still occur. 
Since, in this second stage of the process, the conditions for the competition with neutron 
escape are in several cases more favorable than in the first stage, such effects may give rise to 


AS has been shown in earlier papers,' it is 
possible to explain the principal features 
of the fission of heavy nuclei on the basis of the 
assumption that the process involves a com- 
paratively long-lived intermediate state of the 
compound system in which the excitation energy 
is distributed over all degrees of freedom, as in 
temperature equilibrium. In fact, the excessive 
deformations of the compound nucleus leading 
to rupture are to be attributed to fluctuations 
in this energy distribution occasionally resulting 
in the concentration of a considerable part of 
the excitation energy in particular modes of 
oscillation of the closely coupled system of 
nuclear particles. The probability that a fission 
of the compound system takes place is therefore 
determined by a competition with other disinte- 
grative or radiative processes leading to a 


1N. Bohr, Nature 143, 330 (1939); Phys. Rev. 55, 418 
(1939), and especially N. Bohr and J. A. Wheeler, Phys. 
Rev. 56, 426, 1065 (1939), (refe to in the following 


as BW). 


much increased cross sections for the fission process. 


decrease in the excitation energy of the residual 
system to such an extent that fission is no 
longer possible. 

In the ordinary cases, where the compound 
system has an energy not greatly exceeding that 
required for fission, the occurrence of one or the 
other competing process will reduce the energy 
available below the critical value. If, however, 
the excitation of the compound system is very 
high, the residual system may have a sufficient 
excitation to permit fission. In the second stage, 
the probability of fission will, of course, again 
depend on a competition with other disinte- 
grative or radiative processes. Such progressive 
transformations have already been briefly dis- 
cussed (BW, p. 449) especially in connection 
with deuteron incited fission, but at that time 
no experimental evidence for their occurrence 
was available. Recent experiments on fission with 
high speed neutrons as well as with deuterons 
seem, however, to afford definite evidence of 
successive transformation and, as they at the 
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same time elucidate the competitive character 
of the fission process, it may be of interest here 
to consider them in some detail. 

In the first place, it has been found by Ageno, 
Amaldi, Bocciarelli and Trabacchi? that the 
cross section for fission in uranium by neutron 
impact remains practically constant for neutron 
energies from about one million volts to about 
ten million volts, but that it increases consider- 
ably for neutrons with still higher energies, 
obtained by the bombardment of lithium with 
deuterons. This result can easily be understood 
by considering that, for the lower energies, we 
have simply to do with a competition between 
neutron escape from the compound nucleus 92U™ 
and its fission, the ratio between the probability 
of these competing processes being nearly con- 
stant for neutron energies above 1 Mev. If, 
however, the neutron energy is above 10 Mev, 
there is a considerable probability that the 
residual nucleus 92U** left after neutron escape 
will have a sufficient excitation itself to undergo 
fission. Moreover, the condition for fission is, in 
this case, especially favorable since in 92:U™* we 
have to do with a nucleus of even charge and 
mass number for which, just as in the compound 
nucleus 92U*** formed by neutron impact on the 
rare uranium isotope 92U™*, the critical energy 
for fission is somewhat lower than the binding 
energy of a neutron. Thus, for excitation energies 
of the nucleus 92U™* just above the critical energy 
for fission, no neutron escape can take place; 
and, even at higher excitation, the probability 
for fission in this process is much greater than 
for neutron escape. 

An estimate of the difference between the 
critical fission energy E; and the neutron binding 
energy E, can, for all the nuclei concerned, be 
arrived at by simple considerations (BW, pp. 430, 
433), the result of which may be summarized in 
the following approximate formula: 


SE=E,—E,=0.27(A — 238) —1.32(Z —92) 
—0.6| (A —Z) even 
(A—Z) odd, 
giving AE in Mev for a compound nucleus with 
*M, E. Amaldi, D. Bocciarelli and G. C. 


Trabacchi, Atti Acc. d'Italia, in press; kindly communi- 
cated to the author by Professor Amaldi. 


mass and charge number A and Z. The difference 
in the two constants represents just the difference 
of about 1 Mev between the binding energy of a 
neutron in a heavy nucleus for even and odd 
numbers (A —Z) of nuclear neutrons. 

For and formula (1) gives 
AE=+0.7 Mev, and AE=—0.6 Mev, from 
which one deduces that, for sufficiently high 
excitation the ratio p between the probabilities 
for neutron escape and fission, which for 92U?* 
is about 4 : 1, will for 92U%* be less than 1 : 3. 
For the average cross sections ¢,’ and a,” for 
uranium fission with fast neutrons in one and 
two stages, respectively, we get 

1 1 4 3 3 


4+1 5 


where go is the cross section for formation of 
the compound system 92U”* and the factor a 
roughly represents the fraction of the residual 
nuclei 92U%* having excitation greater than the 
critical fission energy. Since the experiments 
show that the average cross section for D+Li 
neutrons is 40 percent higher than that for 
D+Be and D+B neutrons, a should in this 
case have an average value somewhat larger 
than 5, which appears reasonable from the 
known spectrum of the D+Li neutrons.* 

A comparatively larger effect from successive 
transformations should be expected in the case 
of the fission of thorium with fast neutrons. In 
fact, from formula (1) we get for 9oTh™ and 
90h approximately AE=1.7 and +0.4, corre- 
sponding to values for p of about 24:1 and 
2 : 1. In this case, we therefore get 


=00/25, of’ =8a-00/25. 


With the same value of a as above, we shall 
thus expect that the average cross section for 
D+Li neutrons will be almost twice that for 
D+Be neutrons. Still more conspicuous effects 
should, of course, be expected in the case of 
neutrons with well-defined high velocities for 
which @ can be almost 1, giving a total fission 
cross section nearly 10 times greater than that 
obtained for neutrons with only a few million 
volts energy, for which a is still zero. 


a 935) W. Bonner and W. Brubaker, Phys. Rev. 48, 748 
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In this connection it may be of interest to 
point out that similar effects can also be expected 
in the fission of protoactinium with fast neutrons. 
Here, for the compound nuclei »,;Pa”* and 
91Pa™!, we get from formula (1), 


AE=+0.1 and AE=—1.2 


giving values for p of about 1:1 and 1 : 10, 
leading to 

=400, of’ =Sa-oo/11. 
Thus, we shall expect an increase in fission yields 
of about a factor of two when comparing the 
effects of neutrons of a few million volts with 
neutrons above ten million volts. 

As regards fission induced by deuterons, it 
follows from simple theoretical considerations 
(BW, p. 448) that a sufficient excitation of the 
compound system can only be obtained by impacts 
leading to a complete fusion of the deuteron and 
the original nucleus. The yield of the reaction 
will, therefore, in the first place depend on the 
ease with which the deuteron penetrates the 
electrostatic field around the nucleus, and cross 
sections of the order of magnitude of nuclear 
dimensions can only be expected for deuteron 
energies approaching 10 Mev. This is in agree- 
ment with the experiments of Gant‘ who reported 
a threshold for the process in uranium of about 
8 Mev, without attempting measurements of 
cross sections. Such measurements for uranium 
as well as for thorium have recently been made 
at this Institute by Jacobsen and Lassen® who 
found for both elements a rapid rise in the yield 
of the process between 8 and 9.5 Mev, corre- 
sponding to cross sections of the order of 10-*° 
cm? at the latter energy. 

In deuteron induced fission of uranium and 
thorium, we shall expect progressive transforma- 
tions involving the compound nuclei 9;EkaRe™, 
os3EkaRe*” and »,Pa™, respectively. From 
(1), we get 


AE=—0.4 Mev and AE=-—1.7 Mev 


*D. H. T. Gant, Nature 144, 707 (1939). ; 
5 J. C. Jacobsen and N. O. Lassen, Phys. Rev. 58, 867 


(1940), following paper. 
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for the successive uranium transformations, ang 


AE=+0.6 Mev and AE=-—0.7 Mey 


for the thorium transformations, correspond 
to values of p of. about 1:2 and 1:24 fo 
uranium and about 3 : 1 and 1 : 4 for thorium, 
Thus, we shall expect for uranium 


=309; =8a00/25 


and for thorium 


=3a00/5. 


While o,’ for thorium is only about 3} the valye 
for uranium, we see that the total cross sections 
a;'+0;"" become nearly equal for the two ele. 
ments if a approaches unity, as is to be expected 
from the high excitation (about 15 Mev) of the 
compound system. The experiments of Jacobsep 
and Lassen, showing that for 9-Mev deuterons 
the fission yield for thorium amounts to about 
0.7 of the uranium value, seem therefore to 
confirm that we have here to do with successive 
transformations. 

In the fission effects in uranium with fast 
neutrons and deuterons, the presence of the rare 
uranium isotope 92U*** will, in contrast to the 
case with slow neutrons where it is responsible 
for the whole effect, be of negligible importance. 
In experiments with separated uranium isotopes, 
we shall, of course, also for 92U*** expect to observe 
successive transformations of the same kind as 
those discussed here, but because of the fact that 
for the compound nuclei and »;EkaRe® 
directly formed by neutron and deuteron inm- 
pacts, we have from (1) the large negative values 

=-1.2 and AE=-—2.3, corresponding to 
very small values for p, such effects will be far 
less pronounced than in the cases considered 
above. Similar considerations apply to proto 
actinium fission by deuteron impact. 

From this brief discussion of successive trans 
formations in nuclear fission, it will be seen that 
the study of these phenomena offers a means d 
enlarging greatly the number of different nuclei 
in which the fission process may be investigated. 
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Deuteron Induced Fission in Uranium and Thorium 


I. C. JACOBSEN AND N. O. Lassen 
Institute of Theoretical Physics, University of Copenhagen, Copenhagen, Denmark 
(Received August 12, 1940) 


By means of deuteron beams from the cyclotron, the variation of fission yield with deuteron 
energy is determined for uranium and thorium. In both cases, a rapid rise is found in the 
region from 8 to 9.5 Mev. For an energy of 9 Mev, the cross section in uranium is estimated to 
be 5107?" cm*, and for this energy the ratio between the cross sections for uranium and 
thorium is found to be 0.7, supporting the theoretical expectation that we have in these fission 


processes to do with successive transformations. 


CCORDING to theoretical considerations,! 
A fission of heavy nuclei by deuteron impact 
cannot be expected to occur unless the deuteron 
has sufficient energy to penetrate the Coulomb 
field around the nucleus without breaking up. 
Thus, appreciable effects are not to be expected 
at energies much below ten million volts. This 
conclusion is confirmed by the experiments of 
Gant? who found a threshold for deuteron 
induced fission in uranium of about 8 Mev. 
Gant did not, however, attempt a quantitative 
determination of the cross section for the process. 
In view of the desirability of a closer study of 
the problem* and especially of a comparison 
between the effects induced by deuteron bom- 
bardment in thorium and uranium, experiments 
have recently been made in this Institute with 
the object of obtaining measurements of cross 
sections for fission at various deuteron energies 
in both elements. 

These experiments were made by means of 
the cyclotron built in the Institute in the course 
of the last two years, the construction of which 
will be described in a forthcoming publication in 
the Proceedings of the Copenhagen Academy. 
After various recent improvements, the cyclotron 
is now able to give fairly steady deuteron beams 
of a few microamperes up to about ten million 
volts. The beam has an energy spread of about 
0.5 Mev and has nearly constant intensity over 
this region. 

In the measurements with uranium, four 
aluminum targets upon which had been evapo- 
rated thin layers of uranium (about 1.5 mg/cm?) 


'N. Bohr and J. A. Wheeler, nO C Rev. 56, 426 (1939). 


*D. H. T. Gant, Nature 144, 707 (1939). 
*N. Bohr, Phys. Rev. 58, 864 (1940), preceding paper. 


were mounted in a rotating holder provided with 
absorbing foils. This holder could be inserted 
by means of a vacuum lock into the cyclotron 
chamber and each target in turn exposed to the 
deuteron beam, the foils serving to reduce the 
energy of the deuterons by a known amount for 
each individual target. The fission products were 
collected on aluminum foils fastened in the 
holder in such a way as not to intercept the 
beam. During the irradiation, the holder was 
rotated constantly in order to insure that all 
targets received the same deuteron current. 
After an irradiation period of ten to twenty 
minutes, the target holder was removed from 
the cyclotron chamber. Measurements of the 
beta-ray activity of the foils with the fission 
products were begun after the lapse of about 
ten minutes and continued for some ten hours. 
Separate experiments made to determine the 
background effect due to stray neutrons in the 
cyclotron showed that, for energies above 8 Mev, 
the number of fissions so produced was negligible 
compared with the deuteron effect. 

The variation of fission cross section in uranium 
with deuteron energy is shown in Fig. 1, where 
the open circles indicate differences between the 
measured yield values for energies differing from 
one another by an amount which is less than 
the energy spread of the beam. Thus, the curve 
through these points is an approximation to a 
“differential yield” curve, corresponding to the 
yield to be expected from an infinitely thin 
target bombarded by a monochromatic beam. 
For a determination of the absolute value of the 
cross section, a single point was taken with a 
fixed target to which the deuteron beam current 
was kept as constant as possible. From the 
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- Relative Fission Yield 
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Fic. 1. Relative variation of fission yield with deuteron 
enexgy full circles refer to uranium, open circles to 
t um 


geometry of the collecting and counting arrange- 
ment, we have a means of estimating the number 
of fission processes occurring in the target. 
Assuming the average number of beta-particles 
given off per fission fragment during the counting 
period to be about three, we calculate the cross 
section for fission in uranium of 9-Mev deuterons 
to be about 5X10-*? cm?. We observe from the 
curve that the rise with deuteron energy is very 
rapid in the interval investigated. Thus, the 
experiments indicate that the cross section will 
continue to increase for still higher energies, 
presumably reaching, for an energy somewhat 
above 10 Mev, a value comparable with nuclear 
dimensions, as is to be expected from the theory. 

In order to make a comparison between the 
effects in thorium and those in uranium, two 
targets coated with thick layers of U;0s and 
ThO, were mounted on the rotating holder in 
the vacuum chamber and bombarded alternately 
by a deuteron beam of 9- to 9.5-Mev energy. At 
this energy, it was found that the yield for 
thorium was 0.75 that of uranium. Because of 


the thickness of the layers of uranium and 


thorium oxide on the target, we should, jn 
estimating the ratio of the cross section for 
fission in the two elements for a given deuteron 
energy, take into account that the stopping 
power in the uranium layer for the fission 


fragments escaping from it will be somewhat 


larger in proportion to the amount of heavy 
element present. In fact, because of the greater 
number of oxygen atoms in the U;Og than in 
ThO:, the relative stopping power for fission 
fragments can be expected to be about 10 percent 
greater in the uranium than in the thorium 
target. We may conclude, therefore, that the 
average fission cross section of thorium for 9. 
to 9.5-Mev deuterons is about 0.7 that of 
uranium. As explained in the paper by Bohr 
cited above,’ this result confirms the view that, 
in the fission of uranium and thorium by deu- 
terons, we have to do with successive transfor- 
mations of the compound nucleus. 

The fission yields at various energies for 
thorium were determined in the same way as 
those for uranium, using, however, in this case 
thick targets of ThO:. The results are, in relative 
measure, represented by the full circles in Fig. 1, 
To facilitate comparison of the higher yields from 
the thick targets, the scale is chosen so that the 
values at a deuteron energy of 9.3 Mev for 
uranium and thorium are put equal. We see 
that the curve drawn fits equally well both sets 
of points, indicating that the cross section for 
uranium and thorium rises in very much the 
same way with deuteron energy in the two 
elements, as was to be expected from the theory 
on account of the small difference in charge 
number of the two elements. 

In conclusion, the authors wish to thank the 
director of the Institute, Professor N. Bohr, 
who suggested the problem, for his continued 
interest in the work and for many valuable 
discussions of its theoretical implications. Our 
thanks are due also to Mr. H¢ffer Jensen for 
preparing the uranium targets. 
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Capture Cross Sections for Thermal Energy Neutrons. II 


Franco RASETTI 
Laval University, Quebec, Canada 


(Received September 17, 1940) 


Cross sections for capture of slow neutrons by a number of elements were measured by 
determining the intensity of the 8-activity induced under standard neutron irradiation. The 
number of neutrons was measured by using as standard the activity induced in Mn, the cross 
section for capture by Mn being taken from former experiments. Data are reported for neutron 
capture by AP’, P®, K*, V5, Cu®, Cu®, As®, Br79, Sb!!, 1127, 


INCE the early observations on the activity 
induced by slow neutrons, the activities have 
usually been reported as “strong”’ or “‘weak,”’ but 
more quantitative data are very meager. Rela- 
tively accurate measurements are available only 
for those elements that, like Rh or Ag, possess a 
large cross section for activation. 

It was considered worth while to perform a 
series of measurements on other elements, in 
order to determine the cross section for the 
neutron capture reaction. The use of a radium- 
beryllium source limited the investigation to the 
stronger activities. The work could be greatly 
extended if a stronger source, as afforded by a 
cyclotron, were available. 

The source was placed inside a paraffin block, 
3 cm below the surface. The detectors, 55 cm? 
in size, were placed above the source, on the 
surface of the paraffin, and covered with another 
paraffin block, to increase the activation by 
thermal neutrons. 

The 8-activity was measured by a small, thin- 
walled glass counter, the flat detectors being 
placed along the side of the counter, in an 
accurately reproducible position. The background 
effect was about 8 counts per minute, the 
measured activity seldom giving less than one 
hundred pulses per minute, and often a great 
deal more. Each activity was followed through at 
least one decay period in order to be certain of its 
identity. 

It is believed that the main sources of error in 
the present work are due to the difficulty of 
securing thin uniform layers of the substances 
used as detectors, and to back scattering of 
electrons by the detector itself and by its sup- 
port. Metal foils were used when possible. Com- 
pounds were generally melted on a thin Al or Ni 


foil, and fairly uniform thin layers were obtained 
by tilting the support before the substance had 
solidified. The activity of the Al or Ni support 
was negligible under the circumstances. Na and 
K were rolled into thin sheets between Al foils to 
prevent their oxidation. 

No attempt was made to correct for back 
scattering by the support. Absorption and scat- 
tering by the detector itself and by the counter 
wall were accounted for by assuming an ex- 
ponential ‘“‘absorption,’”’ with an empirically de- 
termined coefficient. This coefficient is given, for 
about the same type of geometry, by Amaldi 
et al.,' except for a few cases where it was 
experimentally determined. The value used is 
given in Table I. No dependence of the mass 
“absorption” coefficient upon atomic number 
was considered, the experimental values generally 
referring to aluminum. 

Since the measurements were effected inside 
the paraffin, it was necessary to ascertain that the 
presence of the detector or its support did not 
appreciably alter the density of the’ thermal 
neutrons. In all doubtful cases, this was checked 
by observing whether the activity of a rhodium 
detector was appreciably decreased by placing 
under it the substance to be investigated. No 
decrease of more than a few percent was recorded. 

As part of the activity is generally due to 
resonance neutrons, a measurement was also 
taken with the detector covered by Cd on both 
sides, and the number of counts thus recorded 
was subtracted from the total. In some cases the 
activity with Cd was too weak to be measured. 

To express the activity in terms of a cross 
section, one must know the number of thermal 


1E. Amaldi, O. D'Agostino, B. Pontecorvo, F. Rasetti 
and E. Segré, Proc. Roy. Soc. A149, 522 (1935). 
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TABLE I. Capture cross sections for thermal neutrons. 


SATUR. 
AcTIV. 
DUETO 


PERIOD 
oF ACTIVE PER 
IsOTOPE 


G/CM? OF 
Cc 


MINUTE NEUTRONS 


PERCENT- 
AGE OF 
AcTIv. 


ACTIVITY 
DUE To 
THERMAL 
RESONANCE NEU- 
TRONS 


x 


2.6 hr. 
15 hr. 
2.4 min. 
14 day 
15 hr. 
3.8 min. 
12.8 hr. 
5 min. 
18 min. 
4.4 hr. 
2.8 day 
25 min. 


om 
~ 


Je 


neutrons traversing the detector per second. This 
can be determined by employing a detector 
whose capture cross section is already known 
from absorption measurements. The capture 
cross section of Mn® =9.4X10~-* cm?) was 
taken as standard since its value is one of the 
best known from the earlier measurements.” The 
cross section o for any isotope is then given by 


Nun A 


’ 

N Amn 
A being the activity (corrected for absorption), 
N the number of atoms of the isotope per cm’, 
Amn and Nya, respectively, the same quantities 
for Mn. 

Table I summarizes the results. The headings 
of the columns are self-explanatory. The cross 
section is referred to the isotope responsible for 
the capture process and not to the whole element. 
The last column gives the values found by 
Lapointe and Rasetti? in those cases where 
comparison is possible. 

The agreement for V* and I'?’ is as good as can 
be expected. For As’*, the present value is about 
half as large as that previously reported. A 
discrepancy by a factor of two cannot be alto- 
gether excluded as due to experimental errors; 


2 C. Lapointe and F. Rasetti, Phys. Rev. 58, 554 (1940). 


however, this discrepancy is more likely to be 
accounted for by the complex structure® of the 
B-spectrum of As’*. The absorption coefficient 
used, 4.5 cm*/g, would roughly correspond to a 
simple distribution with maximum energy of 2.7 
Mev. In addition to this distribution, there 
appears to be a large number of low energy 
electrons, which are scarcely effective in the 
present experiment because of the thickness of 
the detector and of the counter wall. The higher 
value of the cross section from absorption meas- 
urements is, therefore, believed to be more 
accurate. For the other radioactive isotopes 
investigated, it is unlikely that a possible com- 
plexity of the 8-spectrum will introduce a con- 
siderable error in the measurement of the 
activity. For the case of phosphorus, the measure- 
ment of the capture cross section from absorption 
was rather inaccurate, the effect being only about 
five times the statistical error. The present 
measurements indicate a considerably lower 
value, which is probably more accurate. 

The cross sections given for Br” are partial 
cross sections for the production of the two 
radioactive isomers. The difference in the fraction 
of the activity which is due to resonance neutrons 
for each of the two isomers is not significant, 
being within the experimental error. 


3G. L. Weil and W. H. Barkas, Phys. Rev. 56, 485 


(1939). 
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ASSUMED ACTIVITY 
ABSORP- Cor- . 
TION RECTED 
CoEFFI- For AB- 
CIENT SORPTION 
ELEMENT Com- G/cM? 
if AND POUND OF TRONS AND 
IsOTOPE UsED ISOTOPE POUND = cmM?/G COUNTER 
Mn* Mn 0.032 3700 0.04 3500 O04 
i Na* Na 0.07 350 < 
Al?? Al 0.18 750 
Pp P 0.028 ~70 1 
a K 0.003 ~45 
va V.0; 0.012 850 <0 6.8 
Cu® Cu 0.27 100 
Cu® Cu 0.13 1000 0 
As As 0.068 2000 0 6.5 
ta Br?9 NaBr 0.0035 280 0 
Br79 NaBr 0.0035 220 0 
Sb”! Sb 0.073 1060 0 
[27 PbI2 0.009 450 0 
m 
| 
i 
< 
193 
‘} 
i 
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Elastic Scattering in Helium at 15°C 


Ropert LEE Mooney 
Georgetown University, Washington, D. C. 
(Received September 24, 1940) 


An angular distribution function is obtained for the scattered intensity of head-on collisions 
of parallel beams of helium atoms with relative velocities appropriate to a temperature of 15°C. 
Margenau’s interaction energy is used in a form modified for small values of the separation 
parameter. The scattering curve does not exhibit the striking regularity obtained in calculations 
based on the smooth sphere model. The expected maximum at 90° is present. The total collision 
cross section is calculated. The “diffraction cross section’’ is compared with the viscosity 


cross section. 


T has been theoretically demonstrated that 
the angular distribution function, J(6), for 
helium atoms scattered by helium atoms ex- 
hibits maxima and minima similar to those 
observed in ordinary diffraction patterns formed 
by plane waves passing over spherically sym- 
metrical obstacles." Massey and Mohr,? using 
an interaction energy given by Slater and Kirk- 
wood,? have caleulated an angular distribution 
function for elastically scattered helium atoms 
with relative velocities corresponding to 20°C 
and to —185°C. Since the publication of their 
results Margenau has presented a more exact 
expression for the van der Waals forces in 
helium.‘ In this paper is given an approximate 
distribution function obtained by numerical 
integration of the appropriate Schroedinger 
equations with Margenau’s interaction energy— 
modified for small values of the separation 
parameter for the colliding atoms. The “‘diffrac- 
tion cross section’”’ for elastic collisions is com- 
pared with the experimentally obtained viscosity 
cross section. No attempt is made to take into 
account the Maxwellian distribution effect. 
The colliding atoms are moving with equal and 
opposite parallel relative speeds appropriate to a 
temperature of 15°C. 

In setting up the Schroedinger equation it is 
convenient to consider two steady beams of 
infinite width moving with equal and opposite 
gas kinetic speeds. By considering one atom from 
each beam a two-particle system is obtained 


1938) S. W. Massey and C. B. O. Mohr, Nature 130, 276 
?H. S. W. Massey and C. B. O. Mohr, Proc. Roy. Soc. 

141, 434 (1933); 144, 188 (1934). 

as fi Slater and J. G. Kirkwood, Phys. Rev. 37, 682 
‘H. Margenau, Phys. Rev. 55, 1137A (1939). 


with the center of gravity at rest. For such a 
system the appropriate equation is 


(2m*/h)[3m*?— V(r) (1) 


where m*=m/2; v is the relative velocity of 
approach; and V(r), the mutual interaction 
energy with r equal to the separation of the 
colliding particles. The solution of (1) must be 
of the form 


¥=exp (im*vz/h)+(0)r— exp (im*vr/h) 
+exp (—im*vz/h)+f(x—6) exp (im*vr/h)r-, (2) 
in which the first term represents the incident 


wave, and the second the scattered wave as 
observed from the center of gravity when the 


* colliding particles are distinguishable. However, 


since helium atoms are identically similar, it is 
necessary to add the third and fourth terms in 
(2) to secure a symmetrical wave function. It is 
then easily shown that the effective cross section 
for which either particle is deflected into the 
solid angle dw is given by 


\f(0) — 8) |*dw = I(0)dw. (3) 


It has been shown®® that an axially sym- 
metric asymptotic solution of (1) for unlike 


particles is 
(2n+1)i" exp (4) 


n=0 
where L,(r) are the solutions of the ‘‘radial part”’ 
of (1), and », are constants depending upon the 
gas kinetic and mutual interaction energies of the 
colliding particles. For similar particles the 


5 Vide N. F. Mott and H.S. W. Massey, The Theory of 
Atomic Collisions (Oxford, 1933). 
asa Faxen and J. Holtsmark, Zeits. f. Physik 45, 307 
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given by Jeffreys:’ 


L 


Fic. 1., Interaction potential reduced to reciprocal 
Bohr radii squared as a function of the separation of 
helium atoms. 


scattered intensity is given by 
1(6) = (1/48) (4n-+1) 
X (cos 
+E (4n-+1) sin 22nP2x(cos 6) (5) 


In order to evaluate this expression, 72, must 
be known. L,(r) must satisfy 


r2dr\ dr r? (6) 
where k=m*v/h in (5) and (6). Margenau gives‘ 


U(r) = (2m*/h?)(770 exp (—4.60r) 
— X 10-" erg, 


in which rf is in angstroms. In these calculations 
U(r) was modified as shown in Fig. 1 by making 
U(r) constant for 3.50 ao where radius of 
first Bohr orbit for hydrogen. Such modification 
has little effect upon the asymptotic solution of 
(6). The asymptotic forms of L,(r)r— for r large 
were found for different values of m and com- 
pared with sin (kr—4nz). For all n £50, n, was 
determined for a given m from the fact that the 
asymptotic form of L,(r)r~ is® 


L,(r)r~ sin (kr np). (7) 


For n>50, n, was calculated from the expression 


7H. Jeffreys, Proc. Math. Soc. London [2] 23, Part 6 
(1923). 
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+ f [(k?— U(r) —n(n+1)/r*)'—k Yar, 


where 1, is the first (and only) zero of 
[k? — U(r) —n(n+1)/r*]. 


Figure 2 shows 7, plotted as a function of y. 
(It can be shown that 7, is, in general, negative 
for a repulsive field and positive for an attractive 
field.) The phase shift, ., becomes asymp. 
totically less than 90° only for greater than 78. 
The contribution to (@) by terms involving 
for which »,<7/2 and is approaching zero 
asymptotically with reasonable rapidity is rela- 
tively negligible. 

The scattered intensity, J(@) is plotted in 
Fig. 3. In evaluating the many terms necessary 
for the determination of this curve, P,(cos @) for 
n>6 was computed at intervals of 2.5 degrees 
from the approximate formula‘ 


P,,(cos =( ) 


rn sin 


1 1 
x (:-—) sin g—— cot cos } 
4n &n 


where g=(n+4)0+7/4. 
The total cross section for elastic collisions 


was calculated from 


f "1(6) sin 648. 


On account of the orthogonality of the P,, this 
integral is equivalent to the sum 


xz 


Q,=4n/k? & (4n+1) sin? no,. 


Fic. 2. Phase shift as a function of n. 


8E. Jahnke and F. Emde, Funktionentafeln (Teubner, 
1938), p. 117. 
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X-RAY ABSORPTION 


Planimetric integration of the integral expres- 
sion for Q,, performed as a check on the accu- 
racy of the drawing of Fig. 3, agreed with the 
result of the summation expression to within 2 
percent. 

An inspection of Fig. 3 shows that the outer 
boundary of the large central maximum of the 
“diffraction pattern’’ occurs at 6,214°. This 
would correspond, in the optical analogy, to the 
diffraction of a plane wave of wave-length \ by 
a sphere of diameter oz=/6:. Here, \ is given 
by the de Broglie relation, \=A/mv. Using oa as 
defined, one might introduce the “cross section 
for diffraction,” Qa=70.’. 

Q, was found to have a magnitude of 
1228X10-'* cm*. The de Broglie wave-length 
of one of the colliding atoms for a relative ve- 
locity of approach as given above és about 
0.806A. The cross section for diffraction then 
would be given by Qa=d?/0,2=34.21X10-'* 
cm’. Trautz and. Zink® found for helium at 15°C 
a viscosity cross section Q,=14.93X10-'® cm?. 
The ratio Q2/Q,=2.29. It is of interest to note 
that for \=0, the diffraction cross section turns 


*M. Trautz and R. Zink, Ann. d. Physik 10, 81 (1931). 


of 
Fic. 3. Distribution in angle of scattered intensity. 


out to be just twice the classical value, and for 
A= 0, the ratio is 4. 

The most salient feature of the scattering 
curve here presented is that the striking regu- 
larity exhibited by the distribution given by 
Massey and Mohr’ is not duplicated. Naturally, 
the introduction of more complicated potential 
fields into the calculations leads to unanticipated 
results. The maximum of J(6) at 90° found by 
Massey and Mohr for diffraction by smooth 
spheres is present here, and offers an opportunity 
for experimental verification, since it appears at 
all impact velocities. Experimental detection of 
this maximum would lend support to the validity 
of the Bose-Einstein statistics as applied to this 
type of collision. 
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Comparison of X-Ray Absorption Fine Structure for Two Elements 
in the Same Crystal 


S. T. STEPHENSON 
State College of Washington, Pullman, Washington 
(Received August 26, 1940) 


The fine structures on the short wave-length side of the K-x-ray absorption edges of the 
elements Cu, Se, Br, Rb, Sr, and As in the compounds RbBr, SrBr2z, Cu:Se, CuBr, ZnSe, AgBr, 
CuSeO0,-5H,0, CuHAsO;, and Cu(AsO,)2-4H;0 have been investigated with a focusing trans- 
mission spectrograph. The structures of two elements in the same ionic crystal lattice are not 
identical, and some are different even at energy distances of 150 volts from the main edge. 
Some factors other than crystal structure alone are shown to influence the fine structure. 


INTRODUCTION 


HEN atoms are grouped together in a 
crystal lattice, energy zones, characteristic 

of the crystal as a whole and extending sometimes 
several hundred volts out beyond the ordinary 
atomic optical levels, are formed whose existence 


is demonstrated by extended fine structure of 
the Kronig! type to be found on the short 
wave-length side of x-ray absorption edges. 
Secondary absorption maxima occur when elec- 
trons are ejected into outer permitted zones. 


'R. de L. Kronig, Zeits. f. Physik 75, 191 (1932). 
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Fic. 1. Mo Kaia2 doublet (A decreasing to the right) 
Wri foe in the photographic emulsion (A decreasing to 
t t). 


Since these zones are characteristic of the crystal 
lattice, two different elements in the same 
crystal should have identical fine structures. 
This deduction has been verified for a few alloys 
(viz. Zn and Cu in brass?). However, it is not 
true in the case of some ionic crystal lattices. 
The fine structures of: Ca and S in CaS,? K 
and Cl in KCl‘ and KCIO;, Ca and Fe in 
(Ca, Mg, Fe)COs;,5 Ca and Ti in (Ca, Si, TiO;) 
and Fe and Cr in (Fe, Cr)[(Cr, Fe)Oz¢ ]. all differ 
significantly and sometimes completely. In no 
case, however, could the structures for both 
edges be followed out more than 30 volts from 
the main edge and it has been supposed that 
further from the edge the structures would be 
similar. Kronig® has attempted to explain the 


2 D. Coster and R. Smoluchowski, Physica 2, 1 (1935). 
3G. P. Brewington, Phys. Rev. 46, 861 (1934). 

4D. Coster and G. H. Klamer, Physica 1, 145 (1934). 
’V. P. Barton and G. A. Lindsay, Phys. Rev. 46, 362 


(1934). 
*R. de L. Kronig, Physik. Zeits. 36, 729 (1935). 


S. T. STEPHENSON 


TABLE I. Absorption fine structure. Absorption maxima A, B, etc., minima a, 8, etc., 
distances measured in volts from the center of the K edge. 


differences observed for ionic crystals as being 
caused by the presence of dissimilar electric 


fields around anion and cation which give rise 


to different probabilities of transition from the 
K state for the two ions. This explanation js 
only a qualitative one and will remain so untij 
theory is supplemented by more experimental 


work. 
PURPOSE 


It was thought worth while to extend the 
meager experimental data so far available on 
ionic crystals to include a comparison of the 
K-absorption fine structures of two elements 
in the same compound for those crystals con- 
taining elements of higher atomic number than 
the ones cited above. Available crystals con- 
taining elements from Cu to Sr were chosen for 
the test$. It was suggested to the author that 
it would be important in this connection to 
include some homopolar (non-ionic) crystals in 
this study. 


EXPERIMENTAL 


The curved crystal focusing spectrograph 
already described’? was used. Crystals of mica 
0.15-0.30 mm thick, oriented with the direction 
containing the two optical axes parallel to the 
generatrix of the focal cylinder, have given good 
results. More recently a quartz crystal 0.20 mm 
thick has proved even better. It is a laborious 


7S. T. Stephenson, Rev. Sci. Inst. 10, 45 (1939). The 
slit and ionization chamber there described are, in the 
_— work, replaced by the conventional photographic 
m. 


A 


LATTICE 


EDGE 


a B 


Ionic Crystals 
4 13 


Br in AgBr Cubic 20 29 

Rb in RbBr Cubic 11 18 26 46 58(?) 

Br in RbBr Cubic 9 19 32 49 62 104 137(?) 

Sr in SrBr: Cubic 17 29 49 75 101 138(?) 

Br in SrBr2 Cubic 9.5 17 31 47 

Cu and Se in Cu2Se Cubic No repeatable fine structure 

Cu in CuSeO,-5H:0 Triclinic 13 34 68 97 128 160(?) 

Se in CuSeO,-5H,0 Triclinic 12 50 100 150(?) 

Cu in CuHAsO; (?) 12 32 62 91 

As in CuHAsO; (?) 23 52 90 130 163 

Cu in Cu(AsO,)2-4H:0 (?) 14 31 60 

As in Cu(AsO,)2-4H:O (?) 21 53 93 133 157 
Homopolar Crystals 

Cu and Br in CuBr Cubic No repeatable fine structure 

Se in ZnSe Cubic No repeatable fine structure 
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ure to prepare a thin quartz plate of 
proper orientation from a large quartz crystal. 
A suggestion which saved the author much time 
and may be of general interest was that quartz 
crystals already cut for use in radiofrequency 
control might be quite suitable for this work. 
Accordingly, a Y cut commercial radio crystal 
about 0.60 mm thick was ground down to 0.20 
mm thickness with FF carborundum. This 
crystal gave doubled lines.* After etching for 16 
hours in dilute (about 30 percent) hydrofluoric 
acid nice sharp single lines were produced. The 
first microphotometer curve in Fig. 1 illustrates 
the Mo Kaya2 doublet obtained with this crystal. 
Mo Ka has a half-width at half-maximum 
(measured at the same height as the peak of 
Ka;) of 0.20 x.u. showing adequate resolution. 
A metal x-ray tube with either a tungsten or a 
molybdenum target was used at 28,000 volts 
peak (full wave rectification without condensers) 
and 10-50 milliamperes. An aluminum foil about 
0.001 inch thick served as a window. Exposures 
were from ten to one hundred hours on regular 
x-ray films or plates. The absorbers were pre- 
pared by sifting the powdered compound onto 
cigarette papers coated with castor oil and were 
placed close to the curved crystal and between 
this crystal and the film. 


RESULTS AND DISCUSSION 


Representative microphotometer traces are 
shown in the figures. At least two different 


11938) R. Hirsh and J. W. M. DuMond, Phys. Rev. 54, 789 
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RbBr, Cu and Se in CuSeO,-5H:;0, Sr H 
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exposures were made in every case and each 
film or plate was microphotometered at least 
three times. Also blanks of the background were 
taken to eliminate any variations caused by the 
instrument. It was found in this way that the 
peak a in the arsenic compounds is partly 
affected by an unidentified weak line while all 
structure beyond C for these compounds is 
masked by WL7.. 

Since the Br edges of the different absorbers 
might be affected by the AgBr of the photo- 
graphic emulsion, it was necessary to examine 
the Br K-absorption edge recorded in the 
emulsion itself when no absorbers were present. 
A study of the structure of this edge shows that 
it does not affect the results for the Br edges of 
the other absorbers. A typical curve is shown 
(inverted) in Fig. 1. The resolution is sufficient 
to give more detailed extended structure than 
has been observed before.* The close-in structure 
A, a, and B agrees with that previously observed 
with a double crystal spectrometer."° The fact 
that the fine structure is so pronounced relative 
to the main edge may find its explanation in the 
energy zone theory of photographic action." 
Whether or not a quantum is absorbed depends 
on whether an electron, after absorption of the 
quantum and ejection, would find itself in an 
allowed energy zone. If the electron can enter 
an allowed energy zone, it absorbs the quantum 


*J. D. Hanawalt, Phys. Rev. 37, 715 (1931); A. Sand- 
strom and E. Carlsson, Zeits. f. Physik 80, 597 (1933). 

%”S. T. Stephenson, Phys. Rev. 50, 790 (1936). 
"J. H. Webb, J. Opt. Soc. Am. 26, 367 (1936). 
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and is free to travel through the lattice to be 
trapped by an impurity and thus to contribute 
to the formation of a development center aiding 
photographic action. If the electron cannot enter 
a permitted zone, the quantum is not absorbed. 
Therefore, not only are there differences in the 
actual absorption of quanta by the AgBr in the 
emulsion for allowed and forbidden zones, giving 
rise to fine structure in the same fashion as for 
all other crystals, but those electrons going to 
permitted zones in the emulsion are very effective 
in producing development centers aiding photo- 
graphic action and enhancing the absorption 
maxima and giving rise to more pronounced 
fine structure. 

The results for crystals with ionic binding are 
illustrated in Fig. 2. It is apparent from an 
examination of these curves and the data in 
Table I that two elements in the same ionic 
crystal do not have identical structures in 
some crystals (SrBrz, CuSeO,-5H,O, CuHAsO;, 
and Cu(AsO,)2-4H2O). These results differ from 
previous work*~ in two significant ways: First, 
the crystals have elements of higher atomic 
number than those heretofore examined ; second, 
the structure differences are observable much 
farther from the main edge—for example, the 
structures of Cu and Se in CuSeQ,-:5H,0 are 
dissimilar at least as far out as 150 volts from 
the main edge. 

Coster and Klamer* emphasized that differ- 
ences in structure for the two elements K and 
Cl in KCIO; would be expected not only because 
of the difference in charge of anion and cation 
but also because of the fact that the K and 
Cl ions have totally different surroundings in 
the lattice. The same thing can be said of the 
two different ions in these crystals, SrBrs, 
CuSeO,-5H.0, CuHAsOs, and Cu(AsO,)2-4H.0; 
their surroundings are different in every case and 
therefore differences in structure should be 
expected. The surprising feature about the results 
is that the structure differs even so far from the 
edge as 150 volts, for it has been supposed by 
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most workers that the structures would be 
similar in position beyond about 50 volts from 
the edge. 

The case of RbBr is interesting because it jg 
the one ionic crystal for which significant 
differences in structure positions do not appear 
and it is the one crystal of those examined jp 
which each ion would have similar surroundings 
in the lattice; also the outer electron shells fo, 
RbBr should be complete. Therefore, one would 
expect the least differences in structure in this 
crystal for its two edges, in agreement with 
observation. 

In order to compare structures of two elements 
in the same compound unaffected by ionic 
charges, the simple homopolar” cubic compounds 
CuBr and ZnSe were studied. Unfortunately no 
repeatable structure could be observed although 
much time was spent photographing CuBr with 
varying absorber thicknesses under carefully 
controlled conditions. One can only come to 
the very unprofitable conclusion that Cu and 
Br in CuBr are similar, at least, in not showing 
any structure. 

A comparison of the structure positions for 
the Cu edges in CuSeO,-5H,O, CuHAsO,, and 
Cu(AsO,)2:4H,2O shows them to be quite similar. 
The same result is noted if one compares the Br 
edges of RbBr and SrBr2 or the As edges of 
CuHAsO; and Cu(AsO,)2-4H;0. 

In conclusion, these results serve to emphasize 
the growing list of experimental data pointing 
to the fact that, although the crystal lattice and 
its energy zones are a factor affecting fine 
structure in ionic crystals, there are other 
factors of importance, for example, the ionic 
charge and the immediate surroundings of the 
ion from which the electron was ejected. These 
factors seem, from the present work, to be 
effective even further than 100 volts from the 


main edge. 


%W. Hume-Rothery, The Metallic State (Oxford Uni- 
versity Press, 1931), p. 334. 
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Note on Possible Kronig Fine Structure in X-Ray Emission Spectra* 


S. T. STEPHENSON 
State College of Washington, Pullman, Washington 
(Received August 26, 1940) 


Weak copper emission bands are observed extending 34 x.u. out on the short wave-length 
side of Cu Kf2. A comparison with Kronig K-absorption fine structure for copper shows that 
these emission bands are actually caused by absorption of the x-rays while emerging from 


the copper target. 


INTRODUCTION 


URING the course of the experiments de- 
scribed in the previous paper! the K-ab- 
sorption edge and structure of a thin copper 
absorbing screen was obtained as a check on the 
ormance of the apparatus. The extended 
structure exhibited by copper is very pronounced 
and has been the subject of much study. (The 
excellent agreement between the results obtained 
here and those of Swedenborg and Claesson? 
shows how accurately this structure may now be 
measured. Compare columns 1 and 2 of Table I.) 
It occurred to the author that it might be worth 
while to search for emission structure from a 
copper target, of shorter wave-length than the 
absorption edge and in the same region as that 
in which this Kronig type absorption structure is 
found, because of the following considerations. 
In absorption an absorption maximum corre- 
sponds to a permitted energy zone for the 
ejected electron. In emission one could expect 
that electrons already in these permitted zones, 
normally unoccupied, might fall back to a vacant 
K level giving rise to an emission maximum. 
For this to be possible one would have to have: 
first, vacant K levels caused by electron bom- 
bardment of the copper target; second, electrons 
in the higher permitted energy zones of the 
copper crystal lattice caused again by the con- 
tinual electron bombardment of the copper 
target; and third, transition probabilities which 
allow this type of transition. Cauchois* and 
Hayasi* have reported emission bands for Cu 


*S. T. Stephenson, Phys. Rev. 57, 353A (1940). 

1S. T. Stephenson, Phys. Rev. 57, 873 (1940). 

*H. Swedenborg and M. Claesson, Ann. d. Physik 28, 
159 (1937). 

*Y. Cauchois, Comptes rendus 201, 1359 (1935). 

‘T. Hayasi, Sci. Rep. Tohoku 25, 785 (1936-37). 


and Ni which might be explained by these con- 
siderations. 


PROCEDURE AND RESULTS 


The apparatus has been described in the pre- 
ceding paper. For the copper absorption experi- 
ments a tungsten target was used with a copper 
foil absorber 12, thick.** For the copper emission 
experiments a copper target was used with the 
x-rays coming off at an angle of about 20° with 
the plane of the target face. The results obtained 
are illustrated by the microphotometer traces of 
Fig. 1 (curve 1 is absorption and curve 2 is 
emission) and the data in columns 2 and 3 of 


TaBie I. Co absorption and emission fine structure. 
hosted in x.u. from Cu 


PEAKS OR 


VALLEYS PRESENT 


6 
11 


WR 


** Some of the plates showed more marked absorption 


structure than that illustrated in Fig. 1, curve 1. Curve 1 
was chosen for illustration, however, because it was taken 
with the same positioning and dispersion as curve 2, thus 
facilitating comparison between the two. (The tungsten 
et consisted of a tungsten button sunk in the end of a 
hollow copper rod to facilitate cooling. Some electrons 
strike the sides of the rod and give rise to weak copper 
emission which is not strong enough to affect results in the 
region on the short wave-length side of the ee ab- 
sorption edge. Curve 1 is a trace of a twenty-five-hour 
exposure and this copper emission was sufficient, in this 
time, to give the Cu X@;, line whose base is observed.) 
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Fic. 1. Curve 1, copper K absorption with 12, Cu foil 
and W et (25-hr. exposure). Curve 2, copper K 
emission with Cu target and no absorber (70-hr. exposure). 


Table I. The emission structure on several of the 
films was quite visible to the unaided eye. This 
structure did not appear when other targets 
were substituted in the tube and therefore must 
have its origin in the copper target. 


DIscussION 


If the structure from the copper target is 
caused by excited lattice electrons falling into a 
vacant K level as postulated in the Introduction 
above, then the emission maxima (peaks) of 
curve 2 should have the same wave-length as the 
absorption maxima (valleys) of curve 1 and 
similarly the valleys of curve 2 should correspond 
to the peaks of curve 1. A comparison of the 
curves and of the data (columns 2 and 3 of 
Table I) indicates that this is not true. Rather 


the peaks and valleys for the absorption curve 
and for the emission curve come at the same 
places within experimental error (except perhaps 
for A and a, the closest structure’). This may be 
taken as evidence that the structure observeg 
from the copper target is self-absorption of the 
Kronig type in the target itself, brought about 
as the x-rays produced beneath the surface of 
the target pass through a thin layer of copper 
before emerging. 

It must be concluded that the observed struc. 
ture is caused by an absorption process and not 
an emission process and that the considerations 
discussed in the Introduction do not lead to 
emission bands strong enough to be observed jp 
this experiment. Probably the difficulty is that 
there are too few electrons in the upper lattice 
levels caused by the target bombardment. This 
result is considered important because others 
have observed some fine structure bands in the 
emission spectra of Cu and Ni and have specu. 
lated as to their origin. Cauchois* mentions a 
series of soft bands extending to the short wave. 
length side of Cu Kf: and a possible conduction 
electron origin of these bands. She gives no 
curves or data to indicate the exact position or 
extent of this structure. Hayasi* gives curves and 
data for emission structure for copper extending 
out 12 x.u. to the short wave-length side of 
Cu KBs. He ascribes this structure to double elec- 
tron jumps: one electron jumping from an optical 
level to a K level and another between two 
allowed lattice levels in the crystal. From the 
present work it seems probable, however, that 
self-absorption accounts for the structure (except 
perhaps that within about 6 x.u. of Cu Kf) 
observed by both of these authors. 


5 The region within 4 x.u. of Cu Kf: has been examined 
at high resolution with a double crystal spectrometer 
(Beeman and Friedman, Phys. Rev. 56, 392 (1939)) and 
showed emission structure caused by electrons in the 
lower permitted zones in the Cu lattice which are normally 
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PHYSICAL REVIEW 


Isotope Shift in Boron 


Joun P. 
Worcester Polytechnic Institute, Worcester, Massachusetts 


(Received September 7, 1940) 


Mrozowski’s microphotometer curves for the lines \2497, 2498 of B I due to 2°Py, 14y3—3*S; and 
43451 of B II due to (2s)(2p) 'P; — (2p)* 1D. showed asymmetries which he resolved into isotope 
shifts arising from the isotopes 10 and 11. Subtraction of the corresponding normal shift from 
each of his observed values leads to observed specific shifts of —0.373 cm™ for \2497, —0.366 
cm! for \2498, and +0.734 cm™ for 43451. The present note gives a calculation of these 
specific shifts with the aid of Morse, Young, Haurwitz wave functions. The results are —0.3665 
cm for either member of \2497, 2498 of B I and +0.566 cm™ for \3451 of B II. Reasons are 
suggested why the agreement with experiment is better in the case of the B I lines than in the 
case of the B II line, and also why the present calculation for the B I lines agrees better than 
that due to Opechowski and DeVries, who calculated —0.33 cm for 42497, 2498. The possible 
spread due to nuclear spin is estimated with the use of triplet separations and Morse, Young, 
Haurwitz functions and found to be of the order 0.03 cm™ in each case, thus explaining on the 
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basis of smallness the absence of nuclear spin h.f.s. in Mrozowski’s plates. 


ROZOWSKI' has recently investigated 


with high resolving power apparatus the 
boron spectrum excited in a hollow cathode 
discharge. Photometric investigation of the 
doublet lines 42497 and 2498 of BI, due to 
2P,—3*S,; and 2?P,,—3*S;, respectively, and 
of the line 43451 of BII due to (2s)(2p)'P; 
—(2p)? ‘D2 showed in each case an asymmetry. 
This asymmetry Mrozowski interpreted as indi- 
cation of an isotope shift. Resolution of the 
photometric curve of each line into two sym- 
metrical components, each of the typical shape 
due to Doppler effect, led to values for the 
abundance ratio B"/B"® in good agreement with 
mass-spectrograph values. Letting Av=v11— 
we may summarize his results as follows: 


BI, 42497: Av=—0.175+0.012 
BI, 2498: Av=—0.168+0.010 cm—, 
BII, 43451! Av= 0.877+0.004 cm-. 


There was no indication of h.f.s. splitting due to 
nuclear spin. 

It is the purpose of this note to give a theo- 
retical calculation of the expected isotope shifts 
due to nuclear motion and an estimate of the 
possible spreads due to nuclear spin. The nota- 
tion of the author’s paper on isotope shift in 
Mg I? will be followed. Thus for the line AB, 


Av= (Av)w+ (Ar), where the normal shift (Av)y 


1S. Mrozowski, Zeits. f. Physik 112, 223 (1939). 
*J. P. Vinti, Phys. Rev. 56, 1120 (1939). 


where, e.g., 74 is 
the term value in cm~ of level A, and the specific 
shift (Av), = —(1/11—1/10)1838-! 2 Ry(ka—ks) 
= 1.0855(k4—ks) cm. From the observed term 
values we have, to sufficient accuracy, (Ayv)w 
=0.198 cm for either \2497 or 42498 of BI, 
and (Av)y=0.143 cm~' for 43451 of B II. These 
values lead to specific shifts of —0.373 cm for 
2497, —0.366 for \2498, and 0.734 
for \3451. 

For the doublet \2497, 2498 of BI, the level 
A is (1s)*(2s)?3s and B is (1s)?(2s)*2p 
From reference 2, Eqs. (1), (4), (5), and the 
relation ¢= —2(m/M) Ry k, we obtain: 


ka=0 
1s) +J?(2p, 2s) (1) 


In deducing (1) the fact has been used that it is 
valid to calculate ¢ in the SLM representa- 
tion, in which the wave functions in the above 
simple cases are single determinants corre- 
sponding, respectively, to the configurations 
and 1s+1s—2s+2s-2p*m. In eval- 
uating the J’s the variational radial functions 
of Morse; Young, and Haurwitz* will be used. 
Use of the parameters for (1s)*(2s)?2p *P as given 
by them leads, after some simple integrations, 


to J*(2p, 1s)=0.7590 and J*(2p, 2s) =0.2539, 


*P. M. Morse, L. A. Young, and E. S. Haurwitz, Phys. 
Rev. 48, 948 (1935). This paper is eeeany referred 
to as MYH;; the radial factors of the MYH functions are 
written explicitly in reference 2. . 
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with the result that kg =0.3376, so that ka—kp 
= —0.3376 and (Av),= —0.3665 Thus for 
each member of the doublet line (1s)?(2s)?3s 2S 
—(1s)?(2s)?2p ?P there results the calculated spe- 
cific shift —0.3665 cm-, as against Mrozowski’s 
observed values of —0.373 cm~ for \2497 and 
—0.366 cm— for \2498. At this point it should 
be mentioned that a calculation by Opechowski 
and DeVries‘ of this shift gave a value —0.33 
cm~. Their calculation was carried out with the 
aid of Hartree radial functions obtained by 
Brown, Bartlett, and Dunn.* Thus their radial 
functions correspond to the solution of a 
variational problem in which the trial function 
has as much flexibility as is consistent with the 
requirement that it be a product of single-electron 
wave functions. The MYH functions, on the 
other hand, although not so flexible as the 
Hartree functions, take proper symmetry into 
account at the very start. This fact may account 
for the better agreement with observation 
obtained with the use of the MYH functions. 
For \3451 of B II the transition is (1s)?(2p)? Dz 
—(2s)(2p)'P;. For (1s)?(2p)? ‘Dz we use the fact 
that is the same in the representa- 
tion as in the SLJM, representation, and that 
its value in the SLM ,Ms representation is 
independent of Ms or Mz. We may then 
compute its value for Ms=0, M,=2, for which 
the wave function is the single determinant cor- 
responding to the configuration 1s+1s~2p*,2p~,. 
Then from Eqs. (1), (4), and (5) of reference 2 
and the relation connecting ¢ and k, we have: 


ka=§J*(2p, 15), (2) 


where the subscript A refers to the initial level 
(2p)? For the final level B= (1s)?(2s)(2p)'Pi 
we make use of the diagonal sum rule, in the 
manner of reference 2. In the use of this 
method we have to assume that a given single 
orbit function, e.g. 2p, is the same for (2s)(2p)'P 
and (2s)(2p)*P; in the final formulas for the 
singlet and triplet, however, we may use the 
proper (different) values of the parameters as 
found by MYH. This procedure is justified by 
the fact that if we actually wrote the wave 


‘w. nee and D. A. DeVries, Physica [6] 9, 
913 (1939). 

5F. W. Brown, J. H. Bartlett, Jr., and C. G. Dunn, 
Phys. Rev. 44, 296 (1933). 


function in each case as the proper linear 
combination of determinants and then calcy-. 
lated ¢, we should obtain essentially the same 
formulas as would be given by the method of 
sums, the only difference being that the correct 
(different) parameters would appear in the final 
formulas. The method of sums here goes very 
much as it does for (3s)(3p) of Mg I, for which 
the process was given in reference 2, p. 1124, 
thus for the triplet: 


—3(m/M) Ry [J?(2p, 1s) +J*(2p, 2s)] (3) 
and for the sum for singlet and triplet: 
6:+6,=—(4/3)(m/M) Ry J?(2p, 1s). (4) 
By difference, we obtain: 
—§(m/M) Ry [J*(2p, 1s) —J?(2p, 2s)], (5) 
so that 
1s) —J?(2p, 2s) ]. (6) 


Using in Eq. (2) the MYH parameters for 
(2p)?D and in Eq. (6) the MYH parameters 
for (2s)(2p)'P, we find after some simple inte- 
grations that ka =0.8019 and kg =0.2804. These 
values give (Av),=+0.566 cm, as compared 
with Mrozowski’s observed value of +0.734 
cm-', a discrepancy of 23 percent. The fact that 
the agreement is less satisfactory for the B II 
line than for the BI lines may possibly be due 
to the MYH functions. Examination of the term 
values given by the BI and BII functions of 
MYH shows that the wave functions for BI 
are probably much better than those for B II, 
MYH finding indeed that their term value 
discrepancies are always much less when there 
is only one electron in the outermost shell than 
when there are two or more. 

It still remains to be explained why no h.f.s. 
due to nuclear spin is observed in either case, 
since B® and B" both have nuclear magnetic 
moments. The obvious answer would be that the 
spread due to nuclear spin is too small to be 
resolved; this indeed turns out to be the case. 
The spread due to nuclear spin can be estimated 
for the line 3s 2S,;—2p *P; of B 1; for B" it is of 
the order 0.034 for j=}, 0.022 for 
j=%; for B® it would be still smaller. For 
(2p)? 1D.—>(2s)(2p)'P; of B II it comes out 0.031 
cm- for and smaller for Since shifts 
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even of 0.7 cm were not actually resolved by 
Mrozowski, but gave rise only to asymmetries 
in the microphotometer curves, it is clear that 
the splitting due to nuclear spin would be 
entirely unobservable in Mrozowski’s plates. 


Since Opechowski and DeVries‘ have indicated the 
method of estimating the nuclear spin spread for the doub- 
let line of BI, we shall outline only the somewhat less 
obvious estimate for the B II line. We use the notation and 
methods of a paper by Goudsmit,* in which a or b denotes, 
respectively, the interaction parameter in the sense of 
Goudsmit’s Eq. (3) for a single p electron or a single s 
electron with the nuclear magnetic moment, and A (J) is 
the coupling coefficient occurring in the expression A I-J 
for the interaction energy of the electrons with the nucleus, 
where I is the nuclear spin and J the total electronic angular 
momentum. The single electron parameters a or 6 are given 
by the expression . 

Ry = (7) 


where Ry is the Rydberg energy, « the fine structure con- 
stant, g(I) the nuclear gyromagnetic ratio, and r is ex- 

in units of the Bohr radius. For a hydrogenic non- 
penetrating orbit (1/r*)a has the value Z**/n*/(/+ 3)(/+1), 
Z* being the effective nuclear charge. 

For (2s)(2p) we denote the singlet A by A’(1) and the 
triplet A’s for J=1 and 2 by A(1) and A(2). Then the 
method of sums gives: 

A (2) =4a+}), (8) 
A(1)+A‘(1) =2a+}d. 


Applying now the method of sums to the interaction energy 
of the orbital magnetic moments only with the nucleus 
and then to the interaction energy of the spin magnetic 
moments only with the nucleus, we find: 


Aorp(2) = 3a, (0) 
Aom(1) +A? =1}a, 
Agp(2) = —(1/10)a+46, 
Agp(1) +4! (1) = 40, 


where the meaning of the symbols is evident, so A(2) 
=Apr(2)+Asp(2), etc. If we now apply Goudsmit’s 
Eq. (33) for the A(J) of a multiplet, we find: 


A(2) = 48, 
A(1)=}A+o+ 4B, (11) 
A’(1)=,. 


Here 4, ¢, and 8 denote, respectively, parameters for 
interaction with the nucleus of the total orbital angular 
momentum L, of the total spin S, and of the spin of an 
“unbalanced” s electron such as occurs in sp. It is worth 
mention that Goudsmit’s Eq. (33) can be derived very 


*S. Goudsmit, Phys. Rev. 37, 663 (1931). 


easily for all the special cases that are here considered. 
Equations (11) can now be decomposed as follows: 


Aor(2)= Asp(2) = — 48, 
Aor(1)= 4A, Asp(1) =o+ 48, (12) 
Al Aj (1) =0. 


Comparison of Eqs. (9), (10), and (12) leads to the rela- 
tions \=a, o= ja, B= 4b; thus the desired A’(1)=a. For 
(2s)(2p)*P;, therefore, we need only the parameter a. 
For (2p)? we denote the A of the 'D, by A’(2) and the 
A's of the *P by A(1) and A(2). Then, from the method 
of sums: 
A(2)+A'(2) =(8/5)a, 
A(1)=0, 


and from the general multiplet formula: 


A'(2)=), 
A(1)=}A+<e, (14) 
A(2)=}\—}e. 


Comparison of Eqs. (13) and (14) gives \=a and o = —a/2, 
so that the desired A’(2) =a. 

For (2p)? Dz the hydrogenic expression for (1/r*)y gives 
a=5.82(Z*/24)g(I)/1838 cm™. From the work of Mill- 
man, Kusch, and Rabi,’ g(J) is known to have the value 
1.788 for B"; also Z* may be estimated from the triplet 
separations for (2p)**P to be 2.80.* These values lead to 
a=0.0052 cm™ for (2p)? D2. For (2s)(2p) \P: we take Z* 
from the triplet separation for (2s)(2p) *P, viz. 2.82,* and 
obtain a=0.0053 cm. The (1/r*)y values may also be 
estimated from the MYH functions: for (2p)?!D, one 
finds ((1/r*)ay)2p= $(1.385)*=0.883, leading to a=0.0050; 
for (2s)(2p) one finds ((1/r*)y)op= 4(1.335)?=0.793, 
leading to a=0.0045. The two methods check about as 
well as could be expected ;* to obtain an upper limit for the 
spread in each case we take the slightly larger values of a 
obtained from the Z*. Then since A =a in each case, we 
have A =0.0052 cm™ for and A=0.0053 
for (2s)(2p) \P;. Each hyperfine level is then given by 


(15) 


where -F is the hyperfine quantum nimber and. J is the 
nuclear spin, for which the value 1} has been suggested’ 
as the most reasonable value for B". Use of Eq. (15) for 
each level plus the selection rule leads to the hyperfine 
components; the spread is then given by the sum of the 
absolute values of the largest positive and largest negative 


év’s. The result is (0.0078+0.0232) =0.031 cm™. 


7S. Millman, P. Kusch, and I. I. Rabi, Phys. Rev. 56, 
165 (1939). 

SL. Pauling and S. Goudsmit The Structure of Line 
Spectra (McGraw-Hill, 1930), p. 101. 

*For the level 2p *P of B1 an estimate of ((1/r*),y) 
from MYH functions leads to a value of A which is 1 
percent larger than that estimated from doublet separa- 
tions (the method of Opechowski and DeVries). 
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A Theorem on Nuclear Motion in Atomic Spectra 


Joun P. Vintr 
Worcester Polytechnic Institute, Worcester, Massachusetts 


(Received September 7, 1940) 


A theorem is derived concerning the energy shift of an atomic level due to nuclear motion. 
Since the expression obtained (Eq. (12) below) does not involve derivatives, it may be useful 
for calculations of isotope shifts due to nuclear motion in those cases where the radial wave 
functions are known not as analytical expressions but only as tables of numerical values. 


UPPOSE we consider an atom of atomic number Z with m electrons, the energy operator for a 
fixed nucleus being Ho. Then, as shown by Hughes and Eckart,' the effect of nuclear motion can 
be expressed exactly by the addition to Hp of a term 


G= (1) 


where M is the nuclear mass and p; the momentum of the ith electron relative to the nucleus. The 
calculation of isotope shifts due to nuclear motion involves the evaluation of the quantum-mechanical 
average of G for each level considered. It is customary to split G into two parts N and o,? where the 
average N of the “normal” term N can be obtained by use of the virial theorem and the energy of the 
given level, but where the average ¢ of the ‘‘specific’’ term o must be specifically evaluated for each 
level of each atom by the use of wave functions. The present note is the result of an attempt to find 
some general method which would put the calculation of the specific effect on a par with that of the 
normal effect. This attempt, although not successful, led to a simple expression for G which does not 
involve derivatives. It may, therefore, prove useful in calculations of isotope shifts due to nuclear 
motion in cases where the radial wave functions are known not as analytical expressions but simply 
as tables of numerical values. The derivation is very much like one which can be given for the general 


virial theorem. 
a) With the omission of G and spin orbit terms, the Schrédinger equation for a stationary state of an 


q atom can be written: 


~ (h?/8x*m) (2) 


where 
. V=—Ze 1/r;+e2/2 EE (3) 


We first derive a preliminary theorem, Eq. (7) below; m denotes the electronic mass, r, the radius 
vector of the uth electron, and the subscript ,, quaantum-mechanical averages over stationary states. 
To do so, we first apply the operator V, to Eq. (2), form the scalar product with ¥*r, (where the 
asterisk denotes the complex conjugate), and then integrate over configuration space. We obtain: 


— (h*/8x*m) f f f 


Some simple transformations now give: 


f va f (5) 


1D. S. Hughes and C. Eckart, Phys. Rev. 36, 694 (1930). 
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2 J. P. Vinti, Phys. Rev. 56, 1120 (1939). 
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We now sum Eq. (5) on k from 1 to n, transform the resulting first term on the right side by Eq. (2), 
and integrate by parts the resulting second term on the right side. (The last step here, as in the 
derivation of Eq. (5), involves the rejection of a surface integral at infinity and thus restricts our 
proof to discrete states.) We obtain 


= — (8x?m/h?) f V)y*dr—2 f (6) 


k=1 


We now insert Eq. (6) into Eq. (4), cancel a term common to both sides of the resulting equation, 
and use (h/27ri)V,=p,. There results: 


From Eqs. (1) and (7) we find: 
Equation (3) gives: 4 
V,V =Ze*r,/r,?—e? /T in; (9) 
so that 3 
VV=Ze? > iv. (10) 
itr 


The second term on the right side of Eq. (10) vanishes because of cancellation in pairs; physically 
this cancellation arises from the fact that this term is simply the sum of the interaction forces of the 
electrons. We now readily obtain from Eqs. (8) and (10): 


G=Z(e’m/2M)\> (11) 


With the Bohr radius a» and the Rydberg energy e*/2a» as units of length and energy, respectively, we 
may write: 
G=Z(m/M)Q, where (12) 


By use of the virial theorem an expression can now be written down for N as an average of a function 
of the coordinates, which together with Eq. (12) will lead to an expression for ¢ as a similar average. 
Since, however, the resulting expression for ¢ is somewhat more complicated than Eq. (12) for G, 
the usual procedure of finding N by the virial theorem and ¢ by calculation with wave functions would 
not be so simple as direct evaluation of the total energy shift G from Eq. (12). Then for a line A>B 
the shift Av=v(u,) — v(u2), where uw; and we are mass numbers of isotopes, is given by: 


The above derivation shows that Eq. (12) will give the same value for G as will Eq. (1), provided 
that exact wave functions are used. In any actual calculation there may of course be a discrepancy 
between the values obtained by the two methods, due to inaccuracy in the wave function. One 
cannot say on general grounds which value will be the better when such an approximate wave func- 
tion is used; agreement between the values obtained by the two methods might well serve as an 
additional criterion for the accuracy of a wave function. 

Finally it may be remarked that in a paper published some years ago by the writer,* there is 
contained implicitly an entirely different proof of the present theorem. A comparison of Eq. (5) of 
that paper with the equation about one-fourth of the way down the page on p. 436 leads directly to 
Eq. (8) of the present paper. 


*J. P. Vinti, Phys. Rev. 41, 432 (1932). 


58 
ora 
can 
(1) 

The 
nical 
e the 
f the 
each 
» find 
Ss not 
clear 
mply 
neral 
of an 

(3) 
radius 
states. 
re the 
btain: 

(4) 


NOVEMBER 15, 1940 


PHYSICAL REVIEW 


VOLUME 58 


Pressure Effects of Homogeneous Rubidium Vapor on Its Resonance Lines 


Cnr’EN SHANG-YI* 
Norman Bridge Laboratory of Physics, California Institute of Technology, Pasadena, California 
(Received September 23, 1940) 


The broadening of the resonance lines of rubidium in absorption under pressures up to 
150 mm of its own homogeneous vapor was studied by means of a 21-foot grating. Under 
pressures below 1 mm the broadening of the lines was very symmetrical and the line contours 
could be described by the dispersion formula, but when the pressure was high the lines exhibited 
asymmetrical broadening of van der Waals type. The ?P,2 component showed red, while the 
2P3;2 component showed violet asymmetry. Both lines showed the proportionality of the width 
with N, the concentration of Rb atoms. The experimental width at half-maximum is greater 
than that predicted by Houston’s theory by a factor of 1.5 and the width of the transition from 
3P 3/2 is greater than that from *P;2 by a factor of 1.6. A narrow diffuse band was observed near 
the shorter wave-length side of the ?P\;2 component and a similar one near the longer wave- 


length side of the *P3;2. component. 


HE development of the new corrosion re- 
sistant MgO windows! has made feasible 
many experimental researches. Using the new 
windows Hughes and Lloyd? first made an 
experimental study of the resonance broadening 
of potassium resonance lines under pressures up 
to 20 mm. In the present research the corre- 
sponding problem for rubidium was studied and 
observations were made for Rb pressures up to 
150 mm under which the two broadened com- 
ponent lines (with a separation of 147A) began 
to overlap. 


EXPERIMENTAL 


Absorption tubes of 1.5 cm inside diameter 
with MgO windows on both ends were used. 
Three tube lengths were employed, viz. 15, 7.5 
and 0.2 cm. The construction of the 2-mm tube 
is shown in Fig. 1. The tube was made of steel. 
A small capsule containing rubidium was placed 
inside the hole A, and the end of the hole was 
closed by a copper tap screw whose coefficient of 
thermal expansion was higher than that of steel. 
Thus the hole was always tightly closed through- 
out the experiment. B was an annealed copper 
- ring which served as a gasket when the tube C 
was compressed onto the tube D at high tem- 
peratures. The copper ring was too hard, but 


* Research Fellow of the China Foundation for the 
Promotion of Education and Culture. Now at Yenching 
University, China. 

1J. Strong and R. T. Brice, J. Opt. Soc. Am. 25, 207 


(1935). 
2D. S. Hughes and P. E. Lloyd, Phys. Rev. 52, 1215 


(1937). 


had the advantage over aluminum of not alloying 
with rubidium.* The length of the absorption 
tube could be easily determined from the thick- 
ness of the section S. 

The absorption tubes were mounted in a 
vacuum furnace. By means of a sylphon it was 
possible to move the tube C and thereby seal 
off the absorption tube while a vacuum was 
maintained in the furnace. The entire tube was 
heated uniformly and its temperature was kept 
constant within 1°C as indicated by a Chromel- 
Alumel thermocouple in contact with the ab- 
sorption tube. 

Light from a tungsten filament lamp passed 
through the absorption tube and was photo- 
graphed in the first order of a 21-foot grating ina 
Rowland mounting. Eastman type 1 R plates 
hypersensitized with ammonia were used, the 
second-order violet being eliminated with a 
Wratten F filter. The plates were calibrated by 
means of two step-weakeners mounted just in 
front of the plate on both sides of the absorp- 
tion lines. 


Fic. 1. The 2-mm tube. 


? An aluminum ring was tried. When the tube was heated 
to 300°C or higher a large amount of silvery liquid alloy of 
rubidium and aluminum appeared and the tube gave no 


more absorption lines. 
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TABLE I. The width at half-maximum of the two doublet components. 


VAPOR at 

LENGTH TEMPERATURE PRESSURE NUMBER OF ATOMS MAXIMUM (Api/2) X107 (Avij2) @P 3/2) 

INCM KELVIN IN MM PER CC DivipED By N _—_—___-— 

x T P N 3/2) (Avij2) @Pij2) 
15.06 409 1.87 x 10-3 4.42 X10" 1.84 1.41 1.3 
15.06 445 1.14 2.44 10" 1.06 0.62 1.7 
7.50 454 1.35 X10 2.88 X 1.04 0.63 1.7 
7.50 491 9.80 x 10-7 1.43 0.77 1.9 
0.20 563 9.83 x 10! 1.69 x 101° 0.72 0.42 1.7 
7.50 575 1,183 1.99 x 106 1.51 0.84 1.8 
7.50 601 2.361 3.89 x 1016 1.09 0.60 1.8 
0.20 600 2.32 3.75 X 10'* 1.06 0.61 1.7 
0.20 621 3.84 5.67 x 10'* 1.19 0.79 1.5 
0.20 651 7.57 1.1410" 0.76 0.54 1.4 
0.20 660 9.18 1.3410" 1.80 0.63 1.7 
0.20 664 9.98 1.45 x 10!" 0.52 0.36 1.4 
0.20 699 19.9 2.75 X10!" 0.29 0.19 1.5 
0.20 721 29.6 3.97 X 1017 0.16 0.12 1.3 
0.20 766 62.2 7.87 X 10"7 0.11 0.09 1.2 
0.20 797 98 ? 1.2 10" 0.09 ? 0.07 ? 1.3 
0.20 829 152? 1.8 X10" 0.07 ? 0.06 ? 1.2 


The line contours were traced by a Kruss 
microphotometer and were compared with those 
given by the dispersion formula: 


_ PY 
2m { 


in which p is the area under the absorption 
curve, y is the width, in wave-length units, of 
the absorption line measured between points of 
half-maximum absorption coefficient, a is the 
absorption coefficient at wave-length 6A from the 
central maximum, and x is the length of the 
absorption tube. If I/Jo is the fraction of the light 
of wave-length \+ 6d transmitted by the absorp- 
tion tube, this yields 


PY logio 
6A)? = 2 


Therefore a plot with (5A)? as ordinate against 
(logio Jo/Z)~! as abscissa is a straight line with 
slope (py logio e)/27 and intercept (y/2)?. The 
intercept (7/2)? was practically nil. However, 
since 


(1) 


+2 
p= f axd (5d) = (3) 


the width at half-maximum, y, could be calcu- 
lated from the experimentally determined slope, 
the number of Rb atoms per cc, N, the tube 
length, x and the oscillator strength fe; of the line 
components. N was found from the vapor 


pressure data 
logio P= —52.23(76/T) +6.976, (4) 


where T is the temperature and p is the vapor 
pressure in mm. The width at half-maximum was 
then expressed in wave number units (Avy) by 
means of the relationship Avij2=cy/d*. 


RESULTS AND DISCUSSION 


The lines were very symmetrically broadened 
when the vapor pressure was below about 1 mm 
Hg. But for plates taken when the temperature 
of the tube was raised to about 290°C or more, 
the lines exhibited more and more asymmetrical 
broadening of the van der Waals type, and a 
shadow of absorption was observed near the 
broadened line components. 


A. Observations of the line broadening under 
low pressures 


(a) Theoretical formulas for resonance broaden- 
ing.—It is possible to compare the symmetric 
broadening at the lower pressures with that 
predicted by the theory. The Lorentz collision 
theory has been quantum theoretically modified 
by various persons, including Weisskopf* and 
Margenau and Watson.® The result given by 
Margenau and Watson is 


(5) 


‘V. Weisskopf, Zeits. f. Physik 75, 287 (1932). 
20s and W. W. Watson, Rev. Mod. Phys. 8, 
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Fic. 2. The | broadened lines. T=766°K; 
p=62.2 mm Hg; x=0.196 cm. 


where fi2 is the oscillator strength of the line in 
question. The numerical factor 27 is the result 
of a crude approximation and might be in error 
by as much as a factor 2. 

The features of this expression that are subject 
to test are the dependence on N, the approximate 
correctness of the numerical factor, and the de- 
pendence on fiz. In Weisskopf’s work no dis- 
tinction is made between fiz, which is the same 
for both members of the doublet, and fz: which 
is twice as great for one member as for the other. 
The oscillator strength fiz is approximately } in 
this case, and Margenau and Watson indicate 
clearly that fi2 is to be used. 

Another treatment of the problem, with special 
attention to the relative broadening of the 
members of a multiplet, has been given by 
Houston.® His result is 


64S(1, 2) 3x6! eA 
2(2S,-+1) 


The quantity S(1,2) is the line strength as 
defined by Condon and Shortley,’ while J; and J2 
are the inner quantum numbers of the ground 
and the excited states, respectively. This ex- 
pression has the same dependence on WN as has 
Eq. (5), it has a somewhat different numerical 
factor, and it definitely predicts that all members 
of a Russell-Saunders multiplet should have the 
same width, except for a factor equal to the ratio 
of the wave-lengths. 


*W. V. Houston, Phys. Rev. 54, 884 (1938). 
7E. U. Condon and G. Shortley, The Theory of Atomic 


Spectra (Cambridge, 1935), p. 98. 


Nfiz- (6) 


CH’EN SHANG-YI 


(b) Comparison of theory with experiment.— 
In Table I are listed the experimental values of 
the width at half-maximum divided by N the 
number of absorbing atoms per cc. The ratios are 
listed separately for the two components of the 
resonance doublet corresponding to transitions 
from the *P3,2 and the ?P1;2 states. In the range 
below 1 mm of vapor pressure, in which the lines 
are observed to be symmetrical and in which one 
may expect the theory to apply, the width at 
half maximum is found to be proportional to NV 
as predicted by theory. In this range the average 
values of the ratio are 


1.22 4 10-7, 
4 10-”. 


In comparison Houston’s formula (6) predicts 
the ratio 
(Avy;2) (?P3/;2)/N= 0.65 4 10-”, 


(Avy) (?P1;2)/N= 0.66 10-7. 


The differences in the absolute values of the 
experimental and theoretical ratios are probably 
not greater than can be explained by the approxi- 
mations used in the theoretical derivation. On 
the other hand, the great difference between the 
observed line widths of the two components at 
every pressure is definitely not predicted by the 
present theories. 


B. Observational line broadening under higher 
pressures 

When the temperature of the tube was raised 
to 290°C, corresponding to a Rb vapor pressure 
of 1 mm, the absorption lines began to show 
traces of asymmetry which became more pro- 
nounced as still higher pressures were reached. 
The *P1;2 component was asymmetrical towards 
the red and the *P3;2 component towards the 
violet. Figure 2 gives sample curves showing this 
asymmetry. Obviously the whole of these asym- 
metrical lines cannot be represented by the 
dispersion formula (1). If (6A)? is plotted against 
(log Io/I)~ for each side of these lines separately 
it is found that a fairly straight line is obtained 
for the shorter wave-length wing of the *Pyy: 
component and the longer wave-length wing of 
the ?P1;2 component while the opposite wings of 
these components yield strongly curved lines. 
The values of the width at half-maximum divided 


poy 


| 
| 0.3 
| a0 60 BO 
| 
| | 
| 
| 


dicts 


FIELD CURRENTS 887 


by N listed in Table I for the pressures above 
1 mm were calculated from the slopes of these 
straight lines. The values for these higher 

ures do not therefore represent true widths 
at half-maximum. They are merely convenient 
parameters for indicating the form of one wing 
of each component in terms of formula (1). 

The fact that the shorter wave-length wing of 
the 2P1/2 component and the longer wave-length 
wing of the ?P3;2 component deviate strongly 
from the dispersion formula (1) is not only due 
to the asymmetry but is also due to the appear- 
ance at these higher pressures of a faint structure- 
less band near each of the component lines. 
Because the bands are diffuse and faint and are 


located near to the doublet lines the position of 
the bands could not be determined with great 
accuracy. The first of these bands extends for 
approximately 6A to the red of a fairly sharp 
edge at 7964.5A while the second extends about 
10A to the violet of an edge at 7817.6A. While 
these bands have apparently not been observed 
before, similar bands have appeared in the 
absorption spectra of other alkalis. 

In conclusion the author wishes to express his 
gratitude to Professors I. S. Bowen and W. V. 
Houston for their supervision and encourage- 
ment. He is also indebted to Dr. P. E. Lloyd for 
his generosity in discussing his experiences during 
research on K resonance lines. 
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The Energy Losses Attending Field Current and Thermionic Emission of 
Electrons from Metals* 


G. M. FLEmInGf AND JosEPH E. HENDERSON 
University of Washington, Seattle, Washington 
(Received September 10, 1940) 


The energy losses accompanying field emission and thermionic emission have been measured 
in the same experimental tube. This was accomplished by observing the thermal electromotive 
force in a junction from which the emission currents were drawn. No measurable temperature 
change was observed for field emission up to thermionic temperatures. For thermionic emission, 
temperature changes were observed which, when correlated with the power losses from the 
filament as a function of temperature, yield a calorimetric value for the work function of 
tungsten of 4.46+0.09 volts. A theoretical expression is given for the average net energy loss 
per electron emitted in thermionic and field current emission. The experiment gives strong 
additional evidence that in field emission the electrons escape by penetrating rather by sur- 
mounting the surface potential barrier as in thermionic emission. : 


INTRODUCTION 


T has frequently been observed that during the 
emission of large thermionic currents of elec- 
trons there is a very noticeable decrease in the 
temperature of the emitting body. This phe- 
nomenon provides direct evidence that the 
electrons furnishing the emission current partake 
of the thermal energy of the solid. Evidently the 
electrons leaving a metal at the high tempera- 
tures occurring in thermionic emission have more 


* Abstracts: J. E. Henderson and G. M. Fleming, Phys. 


Rev. 48, 486 (1935); 54, 241 (1938); 56, 853 (1939). 
t Now at Russell Sage College. 


temperature energy than those which are sup- 
plied to replace them. The difference in the 
average energy of the electrons emitted by the 
cathode and of those supplied to the cathode 
results in a loss of power and in a corresponding 
decrease of temperature. The power loss. pro- 
ducing the temperature change, when translated 
in terms of the average net energy loss per 
electron emitted, yields a value of the work 
function of the emitting metal in very satisfactory 
agreement with that obtained by independent 
methods.! 


1Saul Dushman, Rev. Mod. Phys. 2, 394 (1930). 
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Fic. 1. Section through the experimental tube. 


Since the knowledge of the energy losses by a 
metal emitting electrons provides information 
concerning the mechanism of emission it is 
important to make measurements in the case of 
field current emission where electrons are sup- 
posed to escape from the metal by penetrating 
rather than by surmounting the barrier as in 
thermionic emission. The only at all satisfactory 
theoretical explanations? of the current-voltage 
relationships for field emission, as well as the 
form of the energy distribution curves* for field 
emission support this mechanism. The almost 
complete independence of these energy distri- 
butions as a function of temperature‘ as well as 
the remarkably small temperature dependence of 
field emission characteristics as obtained by 
Ahearn® and others® makes the region extending 
to thermionic emission of particular interest. The 
present experiment deals with the problem by 
measuring the cooling of a field current emitter 
in this region. The amount of cooling depends, as 


in thermionic emission, on the amount of temper- 
ature energy possessed by the emitted electrons. 


METHOD 


The essential feature of the method is to 
incorporate a thermocouple junction in the 
emitting cathode and to use the changes in the 
electromotive force of this junction to measure 
changes in temperature of the cathode during 
emission. This information coupled with a knowl- 
edge of the physical properties of the cathode is 


— and L. Nordheim, Proc. Roy. Soc. A119, 

3 J. E. Henderson and R. K. Dahlstrom, Phys. Rev. 55, 
473 (1939). 

4 J. E. Henderson and R. K. Dahlstrom, Phys. Rev. 45, 
764 (1934). 

5 A. J. Ahearn, Phys. Rev. 44, 277 (1933); 50, 238 (1936). 
(1926) A. Millikan and C. F. Eyring, Phys. Rev. 27, 51 


G. M. FLEMING AND J. E. HENDERSON 


sufficient to determine the rate at which energy 
is carried away by the emission current. The 
cathode can be designed to permit measurements 
at temperatures extending from room tempera. 
ture into the thermionic range. 


APPARATUS 


Since large power losses must be dissipated at 
the target in this experiment, it was necessary to 
provide a water-cooled anode. Consequently, the 
thermocouple cathode assembly was installed in 
place of the regular cathode of a General Electric 
XP-DF x-ray tube the anode of which is water 
cooled. Figure 1 shows a section through the 
experimental tube. A seven-wire press was in- 
stalled to permit the testing of various cathode. 
thermocouple arrangements. Of the four different 
arrangements that were tried, only the one which 
yielded the best experimental results will be 
described here. Since very good vacuum condi- 
tions are necessary for stable field current 
emission, the tube was evacuated with a two- 
stage oil diffusion pump used with charcoal.’ 
Consistent space vacuums of the order of 10° 
mm of Hg were maintained throughout the 
experiments. 

The potentials necessary for field emission 
were provided by a filtered, valve rectified x-ray 
supply, of approximately 30 kv. To insure steady 
currents, a radiator type x-ray tube was used in 
series with the experimental tube. In the portion 
of the experiment dealing with thermionic emis- 
sion the potential was provided by batteries. 

Of the various cathode-thermocouple combi- 
nations tried, the one illustrated in Fig. 2 proved 
most satisfactory. The thermocouple was formed 
by wires of 0.004-inch tungsten and 0.005-inch 


Fic. 2. Schematic 


diagram of the point 
assembly and meas- 
uring circuit. 
DETAIL A = 
7 J. E. Henderson, Rev. Sci. Inst. 6, 66 (1935). 
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tantalum spot-welded together near A. These 
were bent and shaped to the form of points close 
to the weld and placed opposite the anode to 
serve as the field emitter. Legs 1 and 3 are one 
continuous wire of tungsten and legs 2 and 4 are 
continuous tantalum. The thermal electromotive 
force of the junction formed by legs 2 and 3 was 
measured by a type K potentiometer capable of 
reading to 10~* volt. Either leg, 1 or 4, furnished 
the path from ground for the emitted electrons. 

For temperatures above room temperature, the 
cathode was heated either by a current from a 
battery or from a 500-cycle generator. The 500- 
cycle current was found convenient, since it 
permitted isolation of the measuring circuits 
across various portions of the cathode by means 
of a choke. Temperatures were determined by 
resistance measurements up to where a Leeds and 
Northrup optical pyrometer could be used. Since 
the temperatures of the various parts of the 
cathode were not the same, the reference temper- 
ature is taken as‘ that of the region of maximum 
temperature, tm, along the heated legs. The 
temperature of the junction may be as much as 
200°C from this for the highest temperatures 
quoted. Corrections were made for absorption of 
light by the glass walls by measurements upon 
the temperature of a similar filament observed 
through both walls of the tube. The experimental 
plots in Figs. 3 and 4 show, respectively, the values 
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Fic. 3. The thermal electromotive force, K, of the 
tungsten-tantalum junction of the cathode shown in Fig. 2 
as a function of the maximum temperature, ¢,, along the 
heated wires. The dotted curve shows dK/dt, or the 
thermal electromotive power as a function of tm. 


1500 
TEMPERATURE 


2000 
(°C) 


Fic. 4. The power, P, supplied by the heating current to 
the point assembly as a function of the maximum tempera- 
- tm. The dotted curve shows dP/dt, as a function 
Of ba. 


of the thermal electromotive force, K, as a 
function of tm, and the power, P, supplied by the 
heating current as a function of ¢t». In the same 
figures are shown the derivative curves giving the 
rate at which the thermal electromotive force and 
power supplied by the heating current change 
with the maximum temperature /,,. These data 
were taken when the heating was produced by 
the current supply. In the experiment itself the 
changes in power corresponding to changes in ¢, 
as well as to changes in the e.m.f. of the junction 
will depend upon the distribution of this supply 
to the cathode. For example, during emission of 
electrons the power associated with this current 
is drained from the region of the cathode from 
which emission occurs, thus changing the distri- 
bution of the net power supplied to the cathode, 
and the temperature, ¢,, may not be that 
corresponding to the same total power sup- 
plied solely by the heating current. That is, 
(dP /dtm)(W), where the power change is pro- 
duced through the heating current, is in general 
different from (dP/dt,)(i), where the change in 
power occurs with the change of emission current. 
The relation between these two quantities de- 
pends upon the region from which emission 
current is drawn. 

Use was made of the values of the ir drops in 
the cathode in estimating the distribution of the 
areas active in field current emission. It was 
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Fic. 5. The change in electromotive force produced by 
emission. The full circles represent the emission current and 
“ open circles the corresponding changes in electromotive 
orce. 


found that at room temperatures emission oc- 
curred from the tips of the cathode, very close to 
the junction and nearest the anode; thus, in this 
case, the effectiveness of power changes due to 
emission to produce a change in temperature of 
the junction was twice that of power changes due 
to the heating current. The effect of higher 
temperatures of the cathode was, in general, the 
redistribution of the emission, frequently to 
regions of the surface at which the field strength 
was less than the maximum, and the predominant 
emission was not correlated with the position of 
highest temperature. This shifting with tempera- 
ture was perhaps due to the migration of an 
impurity patch of lower work function or to a 
change of the crystal structure along the surface. 
The estimation of the actual distribution of 
current could be roughly made, but the effective- 
ness of power losses in producing a change in 
temperature of the junction was uncertain be- 
cause of the ambiguity of the position along the 
cathode at which the emission occurred. Corre- 
spondingly the temperature of the portion of the 
metal from which the emission took place was 
uncertain within limits. At the highest tempera- 
tures for which measurements were made with 
field currents the effect of power losses was 
uncertain by a factor of about two, and the 
temperature of the emitting region by about 
200°C. Because of the nature of the results 


obtained, this uncertainty is unimportant except 
for the highest value of the temperature given for 
field emission. 

In the case of thermionic emission both the 
position of emission and the temperature can be 
specified because the emission current occurs 
predominantly at the region of highest tempera. 
ture, thus allowing a calculation of the effect of 
the power losses due to emission. Approximate 
expressions for (dP/dtm)(W) and (dP/dtm)(i) can 
be calculated from an equation obtained giving 
the temperature variations along the wire. The 
relation between them is: 


dP\ 1 
The highest temperatures along legs 1 and 4 
occurred at positions about one-third of the 
length of the leg away from the junction, and for 
thermionic emission the maximum along leg 1, of 
tungsten, was higher than that along the tantalum 
leg by a small amount varying from 30° to 50°C. 
The value of the ir drops in the cathode during 
emission showed that only a small part of the 
total emission current was furnished by the 
tantalum leg, this fraction varying from practi- 
cally zero at lower temperatures to one-fifth at 
the highest temperature used. The way in which 
the temperature of the junction, or the equivalent 
electromotive force, K, varies with t,, for the two 
types of power changes likewise must be known. 
With sufficient accuracy, calculations showed 
that 
(dK /dtm)(t) = (dK /dtm)(W). 


RESULTS 


Data were taken by measuring the difference 
in thermal e.m.f. of the junction with the 
emission current alternately on and off. Figure 5 
shows the results for both field current and 
thermionic emission. The remainder of the curve 
for field currents, omitted in the figure, extends 
for the range from 800°C down to room tempera- 
ture along the axis. The temperatures used in 
plotting are again those of the position of maxi- 
mum temperature along the cathode. The emis- 
sion currents are represented by full circles and 
the corresponding changes in e.m.f. by open 
circles. The field currents were produced by the 
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high voltage supply and were usually adjusted 
to be about 5 milliamperes. The thermionic 
currents, except for the one at the highest 
temperature, were produced with 140 volts 
accelerating potential. The last value, for which 
about 90 volts was used, shows lack of saturation 
by departure from the curve rising exponentially 
at lower temperatures. The changes in electro- 
motive force, K, for emission due to the high 
electric fields were less than the experimental 
error at all temperatures used. In the thermionic 
range the effects were easily measurable. Figure 6 
shows the trend of the points, giving the values of 
AK per unit emission current in microvolts per 
milliampere. The decrease in the ordinate for the 
thermionic range with increasing temperature is 
associated with two characteristics of the cathode. 
First, as the temperature of the cathode in- 
creases, the amount of power change necessary to 
produce a given temperature change increases (see 
Fig. 4), and, second, the thermoelectric power of 
the junction decreases in this range (seeFig. 3). 

The power loss per unit current for field current 
emission is less than that measurable by the 
apparatus. Therefore, only an upper limit can be 
given for the average net energy loss per electron 
emitted. At room temperature the sensitivity of 
the apparatus sets this limit at 0.002 ev/electron. 
At higher temperatures the sensitivity of the 
apparatus is less, and for the highest temperature 
at which field emission can be said to have been 
produced in this case, 1250°C, the energy loss per 
electron enfitted was not greater than”}0.08 
ev/electron. 


FIELD 
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1000 1500 2000 
TEMPERATURE (°C) 


_ Fic. 6. The change in electromotive force per unit emis- 
sion current as a function of temperature for field current 
and thermionic emission. 


FIELD CURRENTS 
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(ELECTRON VOLTS) 


4.52}---- ° 


1600 1800 2000 
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Fic. 7. Values of the work function of tungsten calculated 
from the experimentally measured energy losses found at 
different temperatures. 


The values of AK/i found for thermionic 
emission are used to calculate the associated 
power loss per unit emission current or the 
equivalent energy loss per electron emitted. The 
power change per unit current, AP/i, corre- 
sponding to AK/i is 


AP AK /i 
dtm (dK /dtm)(i) 
1 sdP AK/i 
-—(—)m 
1.16\dtn (dK /dtm)(W) 
where the quantities on the right can be obtained 
at each temperature from the curves in Figs. 3, 4, 
and 6. The result, expressed in watts per 


milliampere, can be alternatively given in terms 
of electron volts per electron by the relation : 


w=10°-AP/i, 


where w is the loss of energy per electron. 

The values of w as a function of temperature 
can be compared with those given by theory in 
the relation: 


w=¢o+3.2kT, 


¢ being the work function of the emitting metal, 
and 7, the absolute temperature. Conversely, by 
use of this relation a value of ¢ may be obtained 
which can be compared with the known value for 
the emitting metal. The values for ¢ obtained in 
this way are shown in Fig. 7. The accepted value 
for the work function of tungsten, 4.52 ev,® is 
indicated. The average of the experimental 
results displayed is 4.46+0.09 ev. 


THEORETICAL DISCUSSION 
It is interesting to compare the values obtained 
by theory with the experimental values of the net 


® See A. L. Hughes and L. A. DuBridge, Photoelectric Phe- 
nomena (McGraw-Hill, 1932), p. 76. 
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energy loss per electron emitted. The energy 
losses occurring during electron emission from 
metals can be computed for a simple triangular 
type of potential barrier such as assumed by 
Fowler and Nordheim.? The calculations are 
made for electrons within the metal having 
Fermi-Dirac statistics. 

The number of electrons per unit volume 
having velocity v in the range dv is 


2m* dv,dv,dv, v*dv 
f(v)dv =—— - 


given by the Fermi-Dirac distribution function. 
The electron parameter yu is determinable from 
the number of electrons per unit volume. Assume 
that the emitting metal has a temperature T and 
is connected to a metal at 0°K from which the 
supply of electrons to the cathode is taken. 
Electrons passing between the regions of different 
temperature undergo a redistribution in energy 
to correspond to the temperature of the region 
occupied. The group of electrons passing from a 
region of lower to one of higher temperature gains 
heat energy from the metal, while the group 
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passing in the opposite direction gives up energy 
to the cooler metal. When no emission occurs 
the net current within the metal is zero and the 
energy lost by the hot portion is gained by 
the cool portion of the metal. During emission the 
net current in the metal is not zero but is equa] 
to the emission current. Those of the emitted 
electrons which possess temperature energy, that 
is, for which E>uy, are replaced by electrons 
originally from the cool supply which have 
acquired heat energy on arriving at the emitting 
region at temperature 7. A net loss of energy 
occurs as a result of this replacement equal to the 
difference between the energy of the emitted and 
the added electrons. 

In order to calculate the heat energy which 
must be supplied to the electrons emitted it js 
necessary to find the average energy of the levels 
in the current supply out of which these electrons 
are raised. The number of electrons per unit time 
of each energy which gains temperature energy is 
proportional to the rate at which electrons of this 
energy enter the warmer metal and to the 
fraction of the number of levels of this energy 
which are vacant at temperature 7. 


The current of electrons of energy E coming from the 0°K supply is computed as follows. The 
number of electrons per unit volume having velocity v in dv which also have velocity v, (in the 


direction of increasing temperature) in dz, is 


f(v)dv-dv,/2v = /h*)vdvdv, 


and the current per unit area in this direction is 


(4am/h®)dEdW 
with W=}mz,. The total current density of electrons of energy E=4}mv? in dE is 


Ww=0 


The fraction of the total number of levels which should be vacant at T is 
/ /kT 4.1) 


and, therefore, the rate at which electrons of energy E increase their energy on entering the region of 


temperature T is 


4am 


d 


E dE. 


The average energy of the levels from which the electrons in the current are removed is 


f dE = kT =p 
0 0 
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In the emitted current the energy distribution is modified by the transmission coefficient of the 
barrier. The current density of emitted electrons at E is given by 


4am 
4 


Diw)aw, 


where D(W) is the transmission coefficient of the barrier for electrons of normal energy W. The net 
loss of energy associated with this current is, then, 


(E—ut1.2kT)dE 
4am (E-ut+ ) f Diw)aw 
0 


J, BAT 4-4 


and the average energy loss per electron emitted in the total current is 


The above expression is to be evaluated for two cases: for field current emission at temperatures up 
to the thermionic range, and for thermionic emission. For field emission Fowler and Nordheim? give 
the transmission coefficient as 


D(W) = (4/C)W*(C— W)! exp [—4«(C— W)!/3F], 


x = 2r(2m)!/h. 


The temperature of the metal from which field currents were produced was always below the range 
for which appreciable thermionic emission occurred, and therefore the upper limit of the integral on 
Ecan be put at C, the height of the barrier. Substituting for D(W), setting .»—- W=x and dropping 
terms of higher power than the first in x in the expansion about x=0, we find 


«x z x 


w= 


4 F 


These approximations are valid for the values of the field used® where the greatest contribution to the 
emission current comes from the range around W=uyz. The net energy loss per electron emitted is 
upon integration : 


1 1 1 
i—akT (1—akT)? 


1 
(= akT 2-—akT 


a= 2xo'/F=1.03 X 10°¢!/F. 
For T=300°K and F=108 volts/cm, akT<1, and 
w= 1.6a(kT)?=0.0025 ev/electron. 


where 


w= 


with 


*J. E. Henderson and R. K. Dahlstrom quote values of the field as high as 10* volts/cm; other estimates range from 
this value down to 107 volts/cm. 
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At 1550°K with F=108 v/cm, akT=0.29, and 


E. HENDERSON 


w=0.07 ev/electron. 


The magnitude of the computed energy loss per electron at 300°K is much too small to be detecteq 
by this experiment, but that calculated for 1550°K is very close to the limit detectable. 
In the thermionic range the accelerating potentials ae were relatively small, and the transmission 


coefficient is chosen as: 


D(W)=0, 


D(W) =D, a constant, 


W<C, 
W2C. 


The energy loss per electron is then readily seen to be 
w=o—3.2kT; ¢=C—u. 


The use of this relation between w and ¢ in obtaining a satisfactory value of the thermionic work 
function of tungsten from the experimental data establishes that the agreement between theory and 


experiment is good. 


DIscUSSION OF RESULTS 


This experiment demonstrates clearly that 
there is no net measurable loss in energy from the 
emitter in field emission over the temperature 
range from room temperature to that tempera- 
ture at which thermionic emission begins to 
furnish more than a negligible fraction of the 
total emission. Corresponding to this temperature 
and above there is a definite measurable energy 
loss. These conclusions show that there is a 
fundamental difference involved in the mecha- 
nism responsible for the two types of emission. 

The calculations show that the Fowler- 
Nordheim theory of metals provides a picture 
essentially in agreement with these results. The 
theory describes field current emission as taking 
place by penetration of the surface barrier and, 
consistent with experiment, predicts a cooling 
effect at room temperatures much smaller than 
that measurable and only slightly within the 
limits of error of the method at 1250°C. In the 
thermionic range, where emission occurs without 
the application of high electric fields, only the 
electrons of normal energy greater than that 
corresponding to the top of the barrier are 
emitted and thus carry away heat energy. On 
this basis from the net energy loss per electron 
emitted, the thermionic work function of the 
emitting metal has been computed, agreeing well 
with the accepted value. 

The relation between the temperature of the 
emitter and field emission has been attacked by 
two other types of experiments. In the first of 


these the striking independence of temperature 
of the total current in field emission was estab- 
lished first by Millikan and Eyring® and later for 
pure metals by Ahearn.’ The other type of 
experiment deals with the energy distribution of 
the emitted electrons. Henderson and Dahlstrom‘ 
show that the energy distribution for the electrons 
emitted under the action of intense electrical 
fields changes only slightly with temperature. 
Their work shows a small but detectable increase 
in the high energy groups with temperature. 
These changes in the high energy groups doubt- 
less represent temperature energy carried away 
from the metal but the small number of electrons 
involved make the total energy they represent 
undetectable in this experiment. 

Both of these experiments differ from this 
investigation in their fundamental approach in 
that they deal with the properties of the emitted 
electrons whereas this experiment deals with the 
effect of the emission upon the metal itself. 
Evidently the majority of the electrons involved 
in field emission are practically independent of 
the temperature of the emitter, at least up to 
that temperature at which thermionic emission 
becomes the predominant form of emission, and 
carry away an inappreciable amount of tempera- 
ture energy above that possessed by the electrons 
that replace them. 

The authors are much indebted to the General 
Electric X-Ray Corporation for the gift of the 
water-cooled x-ray anode which formed the 
basis of the experimental tube used in this 
research. 
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Interactions in the Tungsten Atom, W I, in a Magnetic Field 
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General theoretical expressions are given for the magnetic field dependence of the energies 
and intensities belonging to neighboring levels (regardless of coupling) of an atom. The energy 
dependence is tested experimentally for the levels 21,4482 and 21,453, of W I in fields up to 
4.5 cm= (=9.7 X 10‘ c.g.s.), where the pattern of levels extends more than 30 cm™. The values 
found for the constants are Avop=5.154 cm™, g(2142) = 1.4787, g(214,) =2.5113, k*(214;, 214s) 
=0.0447. The deviations vobs— calc, Obtained from individual determinations at different 
field strengths, are in some instances as great as 0.05 cm™, which is somewhat greater than the 
estimated error in the observations; they may be attributed in part to slight repulsive effects 
from other levels and similar disturbances of the reference lines, which have not yet been 


calculated in detail. 


INTRODUCTION 


N a magnetic field, each level of an atom splits 

into 2J+1 nondegenerate sublevels with mag- 
netic quantum numbers ranging from +J to 
—J. When the Larmor precession frequency 
cn(=el1/4xrmc) is small compared with the no- 
field frequency separations cAvp between the 
level under consideration and neighboring levels 
with the same parity that differ by not more than 
1 in J value, the sublevels are equally spaced, 
symmetric about the field-free level, and linear in 
the field (linear or ordinary ‘“‘anomalous” Zeeman 
effect), if we neglect a certain small effect! 
(diamagnetic quadratic Zeeman effect). When 7 
is not small compared with the Av’s, the de- 
pendence upon the field is more complicated: 
the approximate equidistance, symmetry, and 
linearity are lost through an interaction that 
may be described as a mutual repulsion of sub- 
levels with the same M. Thorough studies have 
been made of the interactions among the levels 
of an LS-coupling multiplet (Paschen-Back effect, 
called “incipient” when the asymmetries are 
slight), but hardly any attention has been paid 
to the more’ general case, presumably because 
levels other than members of the same “‘good”’ 
multiplet rarely lie close enough together to give 


'Shown in references 12 and 23 to be negligible in the 
cases considered here. 


rise to noticeable interactions in the fields com- 
monly used.?-* 

The work reported here is a detailed study of 
the mutual influence of any two levels of an 
atom in a magnetic field, with the experimental 
behavior of the levels 214; and 214, of W I as an 


? For the interesting case of Hg I 6s6d *D,, ‘""’D», which 
lie only 3.2 cm apart, see P. Zeeman, Proc. Amst. Acad. 
Sci. 10, 351 (1907); Physik. Zeits. 10, 217 (1909); P. 
Gmelin, Physik. Zeits. 9, 212 (1908); Ann. d. Physik 28, 
1079 (1909); Physik. Zeits. 11, 1193 (1910); M. M. Risco, 
Physik. Zeits. 13, 137 (1912); J. B. Green and R. A. 
Loring, Phys. Rev. 46, 888 (1934) ; E. U. Condon and G. H. 
Shortley, Theory of Atomic Spectra Pp ey University 
Press, 1935), p. 294. The interaction of these levels is com- 
plicated by their relatively great hyperfine structure; cf. S. 
Goudsmit and R. F. Bacher, Phys. Rev. 43, 894 (1933). 

* Asymmetric Zeeman effect patterns were reported for 
several Mo I and W I lines by R. Jack, diss. Géttingen 
(1908); Ann. d. Physik 28, 1032 (1909); Proc. Roy. Soc. 
Edinburgh 29, 75 (1908). 

‘ Incipient Paschen-Back effect in certain configurations 
2p*nl of Ne I, where the coupling is intermediate, has been 
studied by K. Murakawa and T. Iwama, Tokyo Inst. Phys. 
Chem. Research 13, 283 (1930); J. B. Green and J. A. 
ry Phys. Rev. 56, 54 (1939). 

5’ The qualitative explanation of certain of the asym- 
metries* in W I, and preliminary reports on the present 
quantitative study, were given by J. E. Mack and O. 
“oe. Phys. Rev. 51, 291 (1937); J. H. Roberson and 
J. E. Mack, Phys. Rev. 55, 1126 (1939); 57, 1074 (1940). 

* A surprising intensity anomaly in the Zeeman effect of 
Nil is attributed to interaction between levels 28 cm™'! 
apart labeled 3d*4s5s °F, and 3d%4d *D,, but unfortunatel 
without data on the positions of the sublevels, by H. 
Dijkstra, Physica 4, 81 (1937). 

‘J. B. Green and J. F. Eichelberger, whe Rev. 56, 51 
(1939) discuss the theory of the Paschen-Back effect for 
intermediate coupling. 

The general considerations in this apply to 
hyperfine structure levels if J is replaced by F throughout. 
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example. These levels are particularly suitable 
because, while their nearest neighbor with the 
same parity is almost 1400 cm~ away, they are 
only 5 cm apart, and, moreover, one of them 
has an exceptionally large g value. All the levels 
specifically referred to in this paper are listed in 
Table I;° of course, it is understood that except 
for parity and J, the designations of quantum 
numbers in this complex spectrum are only crude 
approximations. In a prospective paper’ a study 
will be made of several smaller interactions in 
W I and the dependence of the interaction con- 
stants k upon the structure of the atom. 


THEORY"! 


In considering the interaction of levels in a 
magnetic field, we shall assume that it is sufficient 
to take into account the mutual influence of only 
two levels at a time; and” that 


The energy of an atom in a magnetic field 
differs from the no-field energy by an amount 
proportional to 2j,+2s,., where j, and s, are the 
projections, respectively, of the total angular 
momentum and the spin angular momentum in 
the direction of the field, and the summation 
takes place over all electrons. In consequence of 
a general theorem regarding the matrix elements 
of a component of a vector," the matrix elements 
of the contribution of the magnetic field to the 
Hamiltonian of the atom can be written: 


(aJ,M, | Hag | bJ,Mz) =hcka fn 6(M.,M,), (1) 


M, if Je=J. 
where f=4(J?—M?)!, if |Ja—J| =1. 
0, if |Je—Jo| >1. 


Here J stands for the greater of the two values 
J., Js; and M, for the common value M,=M, 
where M is used. The coefficient k,, is the same 


®Q. Laporte and J. E. Mack, J. Research Nat. Bur. 
Stand., in preparation. 

10 J. H. Roberson, in preparation. 

11 We are indebted to Professor Breit for essential help in 
the ———— of this section. 
12 The term on the right is the value of 9 at which the 
diamagnetic quadratic effect would equal the normal effect 
(+5108 cm™ for the lowest, and 5X10? cm™ for the 
highest level listed in Table I; cf. reference 23). 

13Cf. M. Born and P. Jordan, Elementare Quantenme- 
chanik (Julius Springer, Berlin, 1930), Section 29. 


TABLE I. Levels of WI referred to in this paper. The desig. 
nation ‘PF’ for 132 indicates an approximately 


equal mixture of two L values. 

NAME DESIGNATION ENERGY gt 

Do 5d‘6s? 5Do 0.00 0 /0 
S3 5d56s 2951.27 1.975; 

13, 5d‘6s? 13777.71 

18, 5d‘6s? 18082.84 
214,° 5d‘6s6p F2 21448.65 1.4787 
214,° 5d*6s6p 7D, 21453.80 2.5113 
261;° 5d‘6s6p "Ds 26189.11 1.7403 

3532° 35311.46 

415.° 41583.16 
44, "Ds 44919.74 1.9479 


* The g’s given here were determined in this study. In the Cases 
indicated by dashes, it was unnecessary to know the g’s. 


for all the sublevels, characterized by different 
M values, of the two levels. The indices a and b, 
which are used simply to distinguish different 
states with the same parity, may be interchanged 
throughout any expression, here and in the 
discussion below. It will be observed that for a 
diagonal element the coefficient ka. is the ordi- 
nary weak-field g factor of the Zeeman effect. 

The simplicity of the general discussion in the 
previous paragraph arises from the fact that one 
may start with exact wave functions in no field; 
even though these functions are not explicitly 
known, they are known to have the above- 
mentioned properties. 

For the special case of LS coupling, all the k's 
are zero" except in the case of two states within 
the same multiplet, and in that case 


4 Here we assume n to be a good quantum number. Ina 
better approximation there is an interaction closely related 
to the famous sharing of intensity among members of a 
series; E. Fermi, Zeits. f. Physik 59, 680 (1929). Fora 
heavy atom like W the radial functions are not strictly the 
same for the orbits m/;,4 and nl;_;, and not strictly orthog- 
onal for nl; and n’‘l;(n’#n) to such a good approximation 
as in the lightest elements. For this reason, on the one hand 
there may be incipient ‘‘Paschen-Back’’ interactions 
between different series members (although a rough calcu- 
lation shows that these effects are negligible in practical 
cases), and on the other hand the exact validity of Eq. (2) 
is subject to doubt because the calculation of this formula 

resupposes that the radial functions are strictly the same 
or nl,,, and nl;_;. While this departure from simplicity, 
associated with large /s interaction energy, may be properly 
described as a breakdown of LS coupling, it is not always 
recognized as such because it cannot be read directly from 
the elementary vector model; for instance, it occurs even 
in one-electron spectra. ; 

16 P, Giittinger and W. Pauli, Zeits. f. Physik 67, 743 
(1931); cf. E. U. Condon and G. H. Shortley, Theory of 
Atomic Spectra, Eqs. 10°2a, 10°26. Giittinger and Pauli’s 
derivation does not cover the case of more than one 
multiplet with the same LZ and the same S in a configura- 
tion, but Professor Breit informs us that k=0 between 
different multiplets even in this case as a consequence of 
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RJ, J) =14+ +S(S+1) 
—L(L+1)]/2J(J+1), (2a) 


(2b) 


where it will be noticed that k(J, J) is the well- 
known LS coupling weak-field g value and 
k(J, J—1) is the coefficient leading to the well- 
known expression for the Paschen-Back effect. 
The sign of k depends, in general, upon a con- 
vention in the choice of phases for the wave 
functions, but we may, without loss of generality, 
use the positive sign throughout this paper. 

In any case where the values of the coefficients 
for the transformation between the actual case 
and the ideal LS coupling case are known, ka 
may be determined. The interpretation of ka» for 

icular cases, and the discussion of the in- 
verse problem of interpreting actual levels from 
experimental values of ka» is postponed." 

This paper is concerned especially with the 
nondiagonal coefficients ka, ab. Reference 
hereafter to k.» (or for brevity, k) implies a 
nondiagonal coefficient unless the contrary is 
specified. . 

The secular equation from Eq. (1) for the two 
sublevels with the same M, of the interacting 
levels a and 3, has the solutions 


(3) 


where vo’ and vo” are half the sum and half the 
difference, respectively, of the no-field wave 
numbers voa and yo, and g’ and g” are half the 
sum and half the difference, respectively, of the 
weak-field g factors, g, and go. 

Equation (3) is the equation of a hyperbola 
with its center at 


Neent = — vo Mg’’/(M*g"?+-k*f?), (4a) 
Veent = — vo" Mg"? +-k*f?) (4b) 
and asymptotic slopes 


(dv/dn). = Mg’ (4c) 
The wave numbers at cent, Which is the point of 


the orthogonality of the space components of the wave 
functions, following the general treatment of spectroscopic 


terms given in E. Wigner, Gruppentheorie (F. Vieweg, 
Braunschweig, 1931). 


minimum energy separation of the branches, are 
Vmin sep = Veent Rf( (4d) 

The intensity of a line component is 
Ta=(Ioa'x + I})?/(1+x*), (S) 


where Joa, Jo, are the intensity contributions of 
the two sublevels with the same M to the no- 
field lines, and . 


x=kfn/(v—voa— Mgan). 


It will be noticed that x=1 when 7 has the value 
—vo’/Mg"” at which the linear approximation 
would predict the crossing of the sublevels. At 
this field the intensities of two Zeeman effect 
sublines are equal in the important case where 
Ing or Ion is zero, for instance on account of a 
selection rule. When Jo, and Jo are both different 
from zero, the general expression for the in- 
tensity, Eq. (5), contains an ambiguity in sign 
that can be removed in any special case. 

Figure 1A shows the results of the calculations 
for a representative case. The values used for 
the seven constants Ja, J, Za; Zb, Yoa, Yoo, and k 
are those found expérimentally for the levels 
214, and 214, of W I. Although negative values 
of » have no physical meaning, it is convenient 
to draw the curves for the two values of M that 
differ only in sign, on opposite sides of the »=0 
axis, for then they form a single complete hyper- 
bola. This simple property arises from the fact 
that all the results depend, not upon the signs of 
M and » separately, but only upon the sign of 
the product Mn. Thus Fig. 1A shows the com- 
plete hyperbola representing the M=+1 sub- 
levels, and in addition the reflection in the 7»=0 
axis of the (dotted) portion representing M=-+1. 

An approximate solution, valid for incipient 
asymmetry, n<|Avo|, in agreement with second- 
order perturbation theory, is 


voat (7) 


where in any particular case only the “repulsive” 
sign is to be chosen in front of the last term. 
Neglecting this k*f*n? term leads, of course, to the 
linear Zeeman effect. 


APPARATUS AND EXPOSURES 


The study in fields »<2 cm was made at 
Wisconsin. The magnet, containing ? ton of iron, 
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was built in the shape of a 0. The current was 
supplied by means of 48 pancake coils of copper 
tubing in parallel. The tubes were cooled by 
distilled water, pumped through a closed circu- 
lating system and cooled by city water at an 
interchanger. The resistance of the coils was 
about 1 ohm. The field became sufficiently steady 
for Zeeman effect work after about 30 minutes 
of running. The current, which was about 10 
amperes for 7»=1 cm and 160 amperes for 
n=1.9 cm, was supplied by a 200-volt, 200- 
ampere motor generator. The pole pieces con- 
sisted mainly of Armco iron rods 12 cm in 
diameter; but each was furnished with a tip of 
Permendur'® carefully threaded into the pole 
piece with a minimum of clearance, and ending 
in a cone with a half-angle of 56°, truncated at a 
diameter of 6 mm. During the exposures the tips 
were from 3.0 to 4.5 mm apart. 

The source was an arc, enclosed in a Back 
chamber that has already been described'’ and 
was made over for this use. The pole tips were 


protected by quartz plates about 0.5 mm thick, 


each mounted upon a brass ring machined to fit 
the conical tip and to allow the quartz to rest 
almost in contact with the Permendur face; but 
for the strongest fields, requiring the greatest 
crowding, the quartz was omitted from the tip 
next to the grounded electrode. The arc was 
produced across a gap, longitudinal to the field, 
between electrodes made of pure tungsten rods'* 
or 50-mil tungsten wires placed in pairs, filed to 
form plates about 1 mm thick. One electrode was 
placed on a pivoted rod and agitated by a motor- 
driven eccentric shaft to make about 5 contacts 
per second, and an inductance was introduced 
into the arc circuit to help maintain a fairly 
steady current. 

The exposures at Wisconsin were made with the 
Physics Department’s 21-foot Rowland mounting 
and 5.6-inch, 15,000 lines/inch Rowland grating. 
The exposures were made in the second order, 
with a dispersion of about 1.33A/mm. Eastman 
33 and 40 plates were used ; the greatest exposure 
time required, with the strongest fields, was 
about 2 hours. 

The study in fields »>2 cm~ was made at 


16 Kindly furnished by the Bell veins Laboratories. 
17 A. E. Whitford, Phys. Rev. 39, 898 (1932). 
18 Kindly furnished by the General Electric Company. 


Massachusetts Institute of Technology. The 
apparatus and the technique of exposure have 
been described.'® Gratings F and G were used jn 
the first order. In exposures*® 2.4, 3.5, and 4.2, 
and in no others, the * and o components were 
separated with the aid of a Rochon prism. 


MEASUREMENTS AND REDUCTION OF Data 


The lines studied were 21,453 (=D )—214,), 
and vv23,465, 23,471 (=214;,2—442). The former 
was studied principally for »<2 cm and the 
latter, only for 7>2 cm~. The former was care. 
fully examined for self-reversal, with negative 
results. 

Each of the readings recorded on our data 
sheets, of the position of a spectrum line on a 
plate, is the result of several settings to 1 micron 
(0.004 to 0.006 cm) made on the large Société 
Génevoise comparator at Wisconsin. After each 
group of settings the plate was moved, so that 
the next group would be shifted to a different 
point along the length of the line, i.e., to different 
grains of the emulsion. 

A serious difficulty arose in our attempt to 
provide reference lines. If we had attempted to 
use field-free lines of any element as standards, 
it would have been necessary to try to avoid the 
possibilities of error inherent in any attempt to 
superimpose the radiation diffracted by the 
grating from two sources separated in time or in 
space. Instead, we used as actual reference 
points for each line the centers of the patterns 
of two neighboring lines (or only one, as dis- 
cussed in the next paragraph). For Dy— 214, the 
reference lines were v21,533 (=132—353.) and 
v21,359 (= D;— 2613), and for 214), 2—442, »23,500 
(18, — 4152) and v23,237 (=.S;— 2613). Our method 
has the serious defect that any departure of the 
center of the pattern of any reference line from 
the no-field wave number introduces an error 
into the result. For that reason, these pairs of 
reference line patterns were selected in each case 
only after a careful study had shown both 
patterns to be symmetric and their centers to 
have the same separation as the no-field lines, 
within 0.02 cm~. Yet, even for lines thus chosen, 
we have been forced to tolerate the undesirable 


1” F.C. Bitter and G. R. Harrison, Phys. Rev. 56, 15 


(1940). 
2 Each exposure is designated by the first two digits ing. 
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possibility of bodily shifts away from their no- 
field positions as a possible source of error in 
our final result. 

Since every sufficiently intense and well-re- 
solved line between 23,471 and v23,237 was 
found to wander with respect to v23,500 and 
y23,237, and experience showed that vv23,471, 
73,465 could be measured as accurately with 
reference to v23,500 alone as by interpolation 
between v23,500 and v23,237,some of the measure- 
ments upon vv23,471, 23,465 were made with 
reference to v23,500 alone. These measurements 
required the independent determination of the 


Fic. 1A. Typical behavior of two otherwise isolated 
atomic energy levels in a magnetic field, predicted by Eq. 
(3) and found for the W I levels 21,453.80;, g=2.511, and 
21,448.652, g=1.479. B, Sharing of properties between the 
levels shown in Fig. 14, as predicted by Eq. (5) and 
qualitatively found. The thickness of each line (above that 
of the thinnest) is proportional to the intensity of a 
ay ary toa J=0 level, permitted for J=1 but forbidden 
or J=2., 


dispersion curve of each plate involved, for 
different plates from the same grating were 
found to differ as much as 0.1 percent in dis- 
persion. 

In order to eliminate false relative shifts of 
the differently polarized components arising from 
possible imperfections of adjustment in the cases 
of the exposures taken with the Rochon prism, 
all measurements of any line component in these 
exposures were made with respect to reference 
line components similarly polarized. 

Although the wave-length dispersion dd\/dx is 
sensibly constant over any one line pattern, the 
approximation involved in considering the wave 
number dispersion dv/dx to be constant over 
a pattern would have introduced errors (al- 
most 0.002 cm for our widest patterns, 17 
cm-! wide) in our calculations. Of course the 
words ‘“‘symmetry,” “center,”’ etc., in our dis- 
cussion of the patterns, refer to the patterns 
reduced to a wave number scale, rather than to 
the photographic plates themselves. The equa- 
tions used in our routine reductions were as 
follows, where Xo refers to the wave-length in air 
of the no-field line (or of the no-field line in- 
volving the level 214, in the case of patterns 
arising from 214, and 2142, which were always 
considered as single patterns); 5A, 6v, and 6x to 
deviations from the no-field line; and 6»’ to an 
approximation to the vacuum wave number 
deviation, 5v, obtained in an intermediate step: 


= dxdd/dx 
bv’ = (8) 
bv = dv’ / Mair t+do(dv’)?. 


FIELD DETERMINATIONS 


Because no quantitative determinations have 
been made heretofore of magnetic fields as great 
as those used in this work, it was necessary to 
give careful consideration to the question of field 
determination. The procedure finally selected 
was: first, to determine as our principal field 
standard" the field = (4.093+0.006) cm found 
for an auxiliary exposure, calculated as discussed 
below from certain patterns chosen for theo- 


retical reasons from among the lines of several 


*t + here refers to probable error. 
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TABLE II. Individual determinations of sublevels of 2141.2 from exposure 3.5 a 552 a. scripts in the body 
the table give M for the even sublevel of the transition. A = shift of with ay symmetric 
44,2, determined to be —0. 036 cm by adjustment to best agreement with line (1444s). ean indicate lines 
too faint to measure, or missing. 

G +9.089% +0.40071 —8.158° —5.572* — 11.186° — 15.662 
.092*1 408° 171-1 —0.082* 
371% 072° 545-1 173? 
2142-44. F 42871 149° 
115 407° 155-1 351% .654-2 
.1392 .368*! .077° 5507 
Do—214, F 139°+A .456°+A .0799°+A .046°+A .526°+A -119°+A 
mean +9.105 | +0.413 —8.141 +5.360  —0.079 —5.557 — 11.165 — 15.656 


TABLE III. Observed sublevels of Sats: 214, and 214, referred to vo(214,). All quantities are expressed in cm™; uncertainty 


~ 0.010 B—Impossible to 


determine value on account of line 


blends. b— Value with unusually large uncertainty 


on account of partly resolved blends. Exposures listed in brackets are plotted together in Fig. 2. 
f 1.586 +4.068 +0.117 — 3.835 
1.586 +4.080 +0.098 —3.835 
1.588 +-4.093 +0.098 —3.878 
1.646 +4.222 +0.123 —4.013 
1.664 +4.258 +0.123 —4.062 
1.666 +4.277 +0.117 —4.050 
4 1.667 +4.277 +0.129 —4.050 
1.679 +4.295 +0.123 —4.074 
1.693 +4.338 +0.110 —4.099 
| 1.698 +4.350 +0.110 —4.154 
1.915 +4.897 +0.153 —4.627 
2.252 +5.727 +0.179 —5.420 B — 1.885 —5.330 — 8.696 B 
2.426 +6.208 +0.234 —5.821 +2.035 — 1.649 —5.326 —9.015 — 12.276 
3.552 +9.105 +0.413 —8.141 +5.360 —0.079 —5.557 — 11.165 — 15.656 
4.285°  -+11.016 +0.602 B +7.540 +0.966 —5.688 —12.721 — 17.802 
4.533 +11.677 +0.674 —9.947b +8.272 +1.320 —5.795 — 13.307 — 18.5845 


elements introduced as impurities; and second, blackening for measurement. The two directions 


to determine the field for each exposure from 
the separations, assumed to be linear in the field, 
within the W I line v23,237 (=S;— 2613), which 
was especially selected as the reference line from 
theoretical and experimental considerations dis- 
cussed below. The necessity of considering the 
value of a transformation factor between field 
and energy separation is automatically elimi- 
nated by the use of the Lorentz unit (7=1 cm™, 
or about 2.142X10* c.g.s.) instead of the c.g.s. 
unit of magnetic field, because the unit is dimen- 
sionally the same as that in which we measure 
the separation of the levels. 

The impurity lines used to determine our 
principal standard of field were the resonance 
doublets of Ca II, Sr II, and Ba II. The concen- 
tration of each impurity was sufficiently low to 
produce sharp lines without self-reversal. The 
exposure was gauged to bring out the impurity 
lines and the W I reference lines with optimum 


of polarization were not separated. The field was 
calculated by making use only of the five in- 
tervals in each doublet that would have been, 
for weak fields, symmetric about no-field lines. 
This method compensates for all nonlinear effects 
that do not influence the interval within any 
pair of sublevels of a level that differ only in the 
sign of M, i.e., to a good approximation, for the 
diamagnetic quadratic Zeeman effect and for 
incipient mutual repulsion effects such as the 
Paschen-Back effect (cf. the last subsection of 
this paper). The error introduced in the calcula- 
tions by the approximate treatment of the 


- Paschen-Back effect can be seen from Eq. (7) to 


be, to the next approximation, 
2MAdg- 


This error, which we have neglected, amounts to 
no more than 1.3X10-* cm™ in line separation, 
or 2.0X10-* cm~ in 9, in any case. Earlier tenta- 
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Change in k produced by on error of 0.0! cm” in D 


A 


25 

~ 


Fic. 2A, Values of k(214;, 2142) observed for each sublevel in each of a number of fields. Full lines, J=1; 


dashes, J =2; filled 
should be labeled “0.211.” B, Discre 


circles, relatively unreliable values obtained from line blends. The ordinate labeled “0.21” 
NCies, Yobs— ¥eale, between the energy levels observed and those calcu- 


lated from Eq. (1) with the best values of the constants. The scale of ordinates (discrepancy) is spread out 


tenfold compared with that of abscissas (field). 


tive field values calculated from Al, Ca, and Sr 
impurities in two of the exposures listed in 
Table III, were »=3.543 and 4.267 cm~, lower 
by 0.3 and 0.5 percent, respectively, than the 
accepted values ; but these values were abandoned 
on account of the relatively poor observational 
quality of the impurity lines, some of which were 
self-reversed. 

The W I reference line »23,237 was chosen 
primarily for experimental reasons: its great 
width (=5n), the relative slightness of its ob- 


served asymmetry, its stable distance from cer- 
tain other relatively symmetric patterns, and its 
general favorable appearance on the plates. 
A secondary reason for the choice was the fact 
that the theoretically estimated’® nonlinearities 
for both levels are small (although, as in all 
other levels in the spectrum, there are enough 
levels near to each level to make one expect 
some asymmetry). The observed pattern in the 
largest field (expressed as Av/n, where »=4.533 
was 
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The greatest discrepancy between this pattern 
and the weak-field one calculated from the 
accepted g values, g(S3)=1.9751 and g(2613) 
= 1.7408, is 0.022 (~5 microns). 


RESULTS AND COMPARISON WITH THEORY 


Although only the energies and not the in- 
tensities were measured, it may be remarked that 
the intensities were found, by visual estimation, 
to be qualitatively in agreement with Eq. (5) 
and Fig. 1B. 

By means of an auxiliary study, the level 44, 
was found to be very nearly symmetric, with 
g(442) =1.9470. In the calculation, levels Dy and 
44. were assumed to follow the weak-field for- 
mula. 

Table II shows the results of all the readings 
recorded for levels 214;,2 from exposure 3.5. It is 
typical except that the discrepancies (average 
deviation from accepted means, 0.011 cm~) are 
rather larger than in most exposures and that it 
is derived from both the line pair vv23,471, 23,465 
and the line v21,453, the former having been 
photographed simultaneously with two gratings, 
while usually the data for one field were deter- 
mined from one line or pair photographed on one 
plate. In Table II and throughout the paper, a 
subscript (except 0, a, 6) stands for J and a 
superscript stands for M. Table III is a sum- 
mary, showing the final values of the sublevels 
observed in each field. The no-field interval 
between the levels was found to be 


In the calculations discussed below, all line 
components were given equal weight, except for 
partly resolved blends (marked 65), which were 
recorded but given no weight when other read- 
ings were available for the same quantity, and 
unresolved blends (marked B), which were not 
recorded. Each interval was weighted by a 
factor proportional to its magnitude. In averag- 
ing the k’s, each value was weighted by a factor 
inversely proportional to the change that would 
result in k from an arbitrary small alteration of a 
single reading (Fig. 2A, bottom). 

The calculations and the comparison of the 
observed with the calculated sublevels proceeded 
as follows: g2 and g; were determined in that 


order, and then k (called ky), independently for 
each sublevel, for all the exposures with 7>2 
cm; and finally the wave number difference 
Vobs— Veale between the observed and the calcy. 
lated value was determined for each sublevel jn 
each field. Each calculation for one constant used 
the best (the weighted mean) value for each one 
previously listed ; e.g., all the calculations for g, 
used g2= 1.4787. From the linear theory 


82 = (vot? — /4n. (9) 


From the linear theory, modified by an inter- 
action with spur invariance: 


81= (vit! + — — vo) ge. (10) 


The results for the g’s are as follows: 


n(cm~!) 
2.252 B 1.9936 
2.426 1.4754 1.9986 
3.552 1.4792 1.9940 
4.285 1.4785 B 
4.533 1.48115 1.99935 
weighted mean 1.4787 1.9950 


An independent determination of k;” was 
made for each sublevel in each field, where k,;” 


was available, by solving Eq. (3) for k: 

ky" =((vs" — vo’ — Mg’n)? 

” ” M=-+1,0, —1 

— (v0! + Magn)? 2. (iy) 


In addition to these six ky” there were calcu- 
lated for each field the following weighted means: 
three k” from and ko”; two ky, from 
k,;®, and ky; and one & from all six ky”. 
Finally the grand average k was computed from 
all the fields with complete (six) ky” data, i.e., 
from fields 2.2, 2.4, 3.5, and 4.5. The values are 
shown in Fig. 2A. The grand average k is 


k(214;, 2142) =0.211. 


The best value of k squared (formerly® called 
K and given a tentative value of 0.06) is 0.0447. 

The ky,” show small systematic trends, and it 
should be noted that the values of k,™” are con- 
sistently larger than those of k2™. 

The values of v1 were obtained from Eq. (3) 
with the best values of the constants. The re- 
sulting curves are shown in Fig. 1A. The con- 
stants of the hyperbolas, calculated from Eq. (4), 


d 
| 
re 
er 
ex 
ha 

in 
; en 

in 
| | 
we 
ant 
de 
est 
line 
(~ 

be 
hac 
| ( 
tior 
acti 
ass 
lated 

with 
(Tab 
| fields, 
cause 


are: 

M n(center) v(center) (dv/dn) 

0 0.000 — 2.577 cm™ +0.424 


The values of vous— Veale are shown in Fig. 2B. 
The discrepancies are somewhat greater than our 
estimated experimental uncertainty, and are 
almost all positive. 

In interpreting Fig. 2B, two distinct types of 
discrepancy ought to be considered: (1) Experi- 
mental error. Of all the factors possibly contribut- 
ing to the accumulated experimental errors in our 
results, two deserve special mention here: First, 
the wandering of the reference levels. Each 
energy determination for 214;,2 depends upon the 
assumption that each of five reference levels 
exhibits in certain respects the weak-field be- 
havior described in the second sentence of this 
paper.” The fact is that most of the discussion 
in the next paragraph is applicable to the refer- 
ence levels as well as to levels 214; ». Second, the 
repeated application of the same measurements 
in the determination of several constants: For 
instance, the same readings upon »;*! and »,~! 
were used successively in determining g(442), g: 
and k, each of which depends upon its pre- 
decessors, as shown above. Thus, although we 
estimate our probable error in setting upon the 
line components to be not more than 0.010 cm 
(~2.5 microns), the total error in yobs — Veale May 
be considerably larger than if independent data 
had been used for each constant with a conse- 
quent statistical canceling of errors. 

(2) Inadequacy of the calculation. The calcula- 
tion took explicitly into account only the inter- 
action between the two levels under special 
consideration. Each sublevel is subject to other 


"Since the energies were measured with respect to 442 
(assuming the field-independence of the zero point interpo- 
lated from 18, —415, and S;—261;) in fields 7 >2 cm™, and 
with respect to 5D» (assuming the field-independence of the 
zero point interpolated from 132.—3532 and D;—261;3) in 
fields <2 cm™, a false discontinuity might be expected to 
a at the junction of the two regions in Fig. 2B (as 

las in Table III). On the assumption that the effect is 
peeertsonal to 7, the value found for the quantity A 
(Table II, caption) at »=3.5 cm™ leads one to expect the 
magnitude of the discontinuity to be about 0.01 cm=, the 
ordinate for strong fields being less than that for weak 
fields. Indeed such a discontinuity, probably from this 
cause, may be recognized in Fig. 2B. 
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disturbances that we have not tried to evaluate, 
such as the diamagnetic quadratic Zeeman effect 
and incipient interactions with many other levels 
in the neighborhood. The disturbances specified 
can be expressed as a series in ascending powers 
of ». The first and most important member of 
the series is the 7? one, which includes the 
incipient interactions with more distant levels, 
insofar as they can be described by Eq. (7), 
as well as the diamagnetic quadratic Zeeman 
effect?* term: 


Av diamagnetic quadratic = 


Thus the 7? member is, to a close approxima- 
tion, an even function of M. Advantage was 
taken of this property in the determination of the 
constants, as follows: Eqs. (9) and (10) depend 
not upon individual v’s, but only upon differences 
vyt™—y,-™“; and in the determination of the 
grand average k, Eq. (11) was used in effect only 
to determine the values of k,, from pairs vy*™” 
and v,;~”. Consequently the constants gi, ge, and 
k used in the calculations are independent of all 
disturbances that are even functions of M and ». 

No such immunity can be claimed, however, 
for the individual determinations shown in 
Fig. 2. Parabolic terms are not evidently pre- 
dominant in the discrepancies shown in Fig. 2B. 
Nevertheless, it is hoped that the outstanding 
discrepancies between the observations and the 
calculations can be improved, especially with 
respect to the general upward trend in-Fig. 2B 
and the associated fact that the values of k,” 
are on the whole higher than those of ko”, after 
a more detailed study" of the influence of more 
distant levels. 


% J. H. Van Vleck, The Theory of Electric and Magnetic 
Susceptibilities (Oxford, Clarendon Press, 1932), Chapter 6, 
Eqs. 105, 106. The summation, which is to be carried out 
over all electrons, is simplified by the fact that the contri- 
bution of any closed subshell is just 22(r*),,/3. From 
certain equations and susceptibility data (especially for 
Au*) in Chapter 8 the value of Av (diamagnetic quadratic) 
for WI may be roughly estimated to range from about 
+0.003 cm to about +0.005 cm for the levels listed in 
our Table I, leading to absolute changes not greater than 
about 0.002 cm™ for any line with which we are concerned. 
Thus the shifts of the lines under study relative to the 
reference lines are too small to detect. 
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A simple method is suggested which allows ionic mobilities in insulating liquids to be com- 
puted from conductivity-time curves obtained after application of a d.c. potential. The ionic 
radii and the concentration of the corresponding particles can be computed from the mobilities 
obtained. The method is applied to some data of J. B. Whitehead, and ionic mobilities of the 


order 10-7 cm/sec. per volt/cm are obtained. 


1. INTRODUCTION 


HIS paper is concerned with the possibility 
of determining the mobility of ions present 
in insulating oils. This is one of the fundamental 
problems regarding insulating materials in gen- 
eral. Measurements of the conductivity supply 
only the product of mobility and ionic concen- 
tration, and a separation of the two factors is not 
possible without some additional information. In 
aqueous or alcoholic solutions such additional 
information is obtained from osmotic data 
or electrode potential measurements, methods 
which, with very few exceptions,' have not up to 
now been extended to insulating liquids. 

It has long been recognized that the decay of 
the current passing through an insulating liquid 
after the application of a d.c. potential is, partly 
at least, caused by the removal or accumulation 
near the electrodes of the ions present. If this 
assumption holds, it should be possible, by an 
analysis of the current-time curve, to draw 
certain conclusons with regard to the ionic 
mobilities, and from them to the size and concen- 
tration of the free ions. 

K. F. Herzfeld? has applied the theory of ionic 
distribution to a potential distribution curve 
obtained by J. B. Whitehead,* and computed 
values for the ionic mobilities on the basis of 
rather involved equations. Continued interest in 
this matter is shown also by some recent publica- 
tions of a more mathematical nature by B. Gross.‘ 


1 Andrew Gemant, “Electrochemical potentials at the 
boundary copper-insulating oil,’ pees at the National 
Research Council meeting, Cambridge, Massachusetts, 
November, 1939. 

2K. F. Herzfeld, Phys. Rev. 37, 287 (1931). 

See J. B. Whitehead, Impregnated Paper Insulation 
lew York, 1935). 

. Gross , Phys. ag 57, 57 (1 —-. Co Compare also a 
i. paper by J. L . Oncley and ollibaugh, 
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A simple theory allowing ionic mobilities to be 
computed from the initial current-time curves 
will be presented. The greatest difficulty jn 
interpreting such curves is the fact that several 
mechanisms are overlapping. Because of the 
widely differing time constants of these single 
processes it is possible, however, approximately 
at least, to separate several regions of the 
current-time curves, corresponding to different 
mechanisms. 

Orientation of dipoles usually has a very small 
time constant, and the separation of this process 
is generally easy. An identification of charging 
currents, corresponding to internal surfaces is, on 
the other hand, a hard task, and a definite 
solution of this problem does not yet exist. 

Electrolytes dissolved in homogeneous media, 
with which this paper is concerned, can be of two 
kinds. First there are ions which can be con- 
sidered as pertaining to an electrolyte in a state 
of “‘total dissociation.’’ These are generally large 
ions of a more or less colloidal nature. In re- 
stricting our consideration to liquids having a 
total charge of zero, these large ions must have 
their oppositely charged equivalent counterpart 
either in equally large, or smaller ions. Second, 
there are electrolytes in a state of partial 
dissociation of a type investigated by Kraus' and 
co-workers. Ions of this type are of a molecular 
size. The first type is relatively quickly removed 
from the liquid and is expected to cause a sharp 
drop in the conductivity after application of a 
d.c. potential in a comparatively short time. The 
second type, however, is continuously regenerated, 
and its effect upon the conductivity is, therefore, 
presented at the Summer Convention AIEE, 1940. They 


mention that ionic mobilities can be computed from current: 


time curves, but do not give any ape. 
. A. Kraus, Science 90, 281 (1939) 
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IONIC MOBILITIES IN LIQUIDS 


either negligibly small, or just a very slow 
decrease. A separation of the two types appears 
possible on this basis. 

The dissociation process will not be discussed 
here. It is the fast drop in conductivity with 
which this paper is concerned, and it should be 
shown how the mobilities of the respective 
charge carriers can be computed from that part 
of the current-time curve. | 

The ions accumulating near the electrodes 
might either be discharged or persist as charged 

icles. Which of the two cases actually prevails 
will essentially depend on the chemical nature of 
the ions. In case they persist, it is to be expected 
that a space charge develops, and, after removing 
the voltage, a discharge current flows. Whereas 
the space charge was detected by J. B. White- 
head, as also by others, the discharge current 
to be more elusive. The fact that the 
original value of the conductivity is sometimes 
restored after removal of the d.c. potential® 
clearly indicates that the ions are not discharged, 
but persist, at least to a certain extent, in a 
number of cases. 

In deducing the equations controlling the cur- 
rent-time curve, this accumulation of ions near 
the electrodes has to be taken into account. One 
can see that there will be two reasons for the 
decrease of the conductivity. One is the decrease 
of the number of carriers in the liquid phase ; the 
other is a decrease of the voltage drop across the 
liquid, since part of the voltage is used up in the 
polarization voltage at the electrodes. This latter 
effect is analogous to that encountered in the 
electrolysis of aqueous solutions. 


2. GENERAL EQUATIONS OF CURRENT FLOW 


Let us consider a condenser of 1 cm? cross 
section and separation a. The total initial con- 
ductance A, of the cell is the sum of Ao, the 
totally dissociated part referred to above, and A,, 
the final conductance, the change of which 
should be negligible. Thus 


A,=Aot+Ay 
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and, similarly, the conductance at any time ¢ 
after application of a d.c. voltage Ep is given by a 


sum: 
A'+Ay. 


Considering the term A’, let us take the 
general case that the ions removed accumulate 
near the electrodes forming a double layer of 
capacitance C and voltage drop E£,. The total 
charge removed is then CE., corresponding to a 
loss of CE./a in charge concentration. This again 
corresponds to a loss of mCE./a in specific 
conductivity, if m=mobility of ions. The loss in 
conductance is then mCE,/a*. We have then the 
relation 

A'=A,—bCE, 
with 
b=m/a’. 

It is assumed here that the remainder of the ions 
is always redistributed uniformly by means of 
diffusion, an assumption which holds only ap- 
proximately. The charging of the capacitance C 
will then be controlled by the differential 
equation 


CdE./dt = (Ay—bCE,)(Eo—E.). 
In replacing E. by E according to 
E=E,—E., 
where E is the voltage drop across the liquid, the 
differential equation can be solved. With the 


boundary condition E = Ey for t=0, the equation 


gives: 
(Ao—bCEo) 


Ao—bCEce* 


where 
B=Ao/C—bE. 


The total current is now given by 
E(A’'+Ay,) 
and the apparent conductance A by: 
A=E(A'+A;)/Eo. 


Putting the above value of E into the equation 
for A, we obtain 


*J. Sticher and D. E. F. Thomas, Elect. Eng. 58, 709 (1939). 
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CURRENT 


DIAGRAM OF CURRENT-TIME CURVE OF 
INSULATING OILS FROM DATA OF J.B. 
WHITEHEAD 3 


LOGARITHM OF TIME 
Fic. 1. 


Equation (2) shows the following peculiarities: 
For t=0 
A =A, 
as required. 

For t=, two cases must be considered, 
according as @ is larger or smaller than zero. If 8 
is positive, 

A)/bC> Eo. 

The expression Ao/bC is the maximum possible 
polarization voltage corresponding to the case 
that all ions contributing to Ao had left the liquid. 
A positive 8 then means that the liquid will 
finally be deprived of any finite electric field. 
Equation (2) for ‘= « gives in this case, A =0, as 
expected. 

On the other hand, for a negative 8, Eq. (2) 
gives, for t= ~, 


A= (Eo—Ao/bC)As/Eo. 


This again is according to expectations, since the 
bracket contains the final voltage drop across the 
liquid. 
Another useful equation to be derived from 
Eq. (2) is the following: 
(dA —Ao(bEo+A./C), 


which is the slope of the A (¢) curve at the origin. 


3. SPECIAL EQUATIONS 


Equation (2) shows that it has to be used in its 
complete form only if the term Ao/dC is of the 
order of Eo. If the former is, at least, say, ten 
times larger than Eo, then the change of the 
voltage drop is the dominating factor, the process 
then amounts to charging a condenser across the 
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liquid as resistance.’ This special case will not be 
considered here any further. 

If, on the other hand, Ep is, at least, say, ten 
times larger than Ao/dbC, or in other words, if 


bCE)/A o> 1, (3) 


then the decrease of the ionic concentration is the 
chief factor. Equation (2) then becomes con. 
siderably simplified. In order to show that Eq, 
(3) usually holds, it shall first be given the form: 


meF/4mdoo> 1, (3a) 


where e=dielectric constant of liquid, F= applied 
field, \= thickness of double layer accumulating 
at the electrodes, oo = specific conductivity, corre- 
sponding to the large ions, being equal to Aw, 
Taking likely orders of magnitude, we now have 
m~10-7 cm/sec. per volt/cm, e~3, 
volts/em, \~10-* cm, abs. e.s.u., and 
the left-hand side of Eq. (3a) becomes 25. This 
means that the second approximation will gener- 
ally hold, and we shall now consider the simplified 
equation referring to this case. 

Equation (2), with Eq. (3) valid, takes the 


form 
A=Ace*'+A; (4) 


or, writing Eq. (4) in terms of specific con- 
ductivities ¢, and replacing b by m/a’*: 


c= oe + (4a) 


Equation (4), as does Eq. (2), operates with 
one mobility m, which seems to be an unlikely 
specialization. The equation can be generalized 
easily by taking account of several ions with 
different mobilities. On the other hand, it is 
always desirable to operate with equations witha 
minimum of constants. Thus Eq. (4) should be 
extended only to two different mobilities, those 
of the positive and negative ions. In this case the 
relative factors of the two terms, generally 
unknown, also can be computed. If m and m’ are 
the mobilities, the respective transference num- 
bers are m/(m-+m’) and m’/(m+m’). Conse- 
quently we have, instead of Eq. (4a): 


m m 
~e~mEot/ a? Eot/a?. (5) 


oo m+m’ m+m’ 


7 Compare S. F. Acree, E. Bennett, G. H. Gray and H. 
Goldberg, J. Phys. Chem. 42, 871 (1938). 
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IONIC MOBILITIES IN LIQUIDS 


Equation (5) is the final answer to our question, 
for with it ionic mobilities may be computed 
from curves. 

It should be pointed out here that according to 
Eq. (5) the o(¢) curve depends upon the voltage. 
If conditions were such as to make the expression 
bCE,/Ao smaller than unity (say for small field 
intensities or large initial conductivities), then 
the o(¢) curve should become independent of the 
potential applied. This circumstance supplies an 
experimental method to differentiate between the 
two cases. 

In order to estimate the size of the ions, it first 
has to be decided whether the Stokes equation or 
the Helmholtz cataphoresis equation holds. This 
depends upon whether the thickness of the 
double layer around the ions is larger or smaller 
than the particle radius. The thickness of the 
double layer (ionic cloud) around the particles is 
given by the expression :* 


where R=gas constant, T7=absolute tempera- 
ture, F=Faraday constant, c=molar ionic con- 
centration per cc, and if c is estimated from the 
residual conductivity to about 10-4, then the 
thickness becomes of the order 10-* cm. Since the 
particle size is much smaller than that, it follows 
that Stokes equation will control the motion. 
Applying this equation, we have for the ionic 
radii : ‘ 


r=en/6xnm, (6) 


where e=electronic charge, »=valency of ion, 
n=viscosity of liquid. 

Finally the number of particles can be com- 
puted from oo and the mobilities obtained, as 
outlined in the introduction. Evidently this 
number is given by: 


(7) 


oo/en(m+m’). 


4. NUMERICAL APPLICATION OF THE THEORY 


The last equations now shall be applied to 
experimental data on mineral oils obtained by 
J. B. Whitehead? and co-workers. Figure 1 shows 
diagrammatically a typical current-time curve. 

The part 1—2 is a decay of very short duration 
(around 10-* sec.) which is possibly due to the 


*A. Gemant, J. Phys. Chem. 43, 743 (1939). 
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orientation of dipolar constituents. At 2 the 
current becomes constant until about 1 sec. and 
then starts to drop again fairly rapidly until the 
region indicated by 3 (around 10° sec.) is reached. 
The region 3—4 is one of a very slow continuous 
decrease where the curve frequently runs nearly 
horizontally with a scarcely perceptible slope. 
Thus we might assume that the portion 2-3 has a 
mechanism of its own, and corresponds to the 
above discussed removal of large ions. If 3—4 has 
a slope different from zero, then this decrease is 
due to a mechanism different from that between 
2 and 3, and very likely corresponds to electrolytes 
in a state of partial dissociation. 

Figure 2 shows a plot, carried out on White- 
head’s data, referring to his oil B-4 at 60°C. 
The abscissae are seconds, the ordinates are 
logio ¢o/(¢—«,), oo, and oy being taken from the 
experimental data. It can be seen that the curve 
gives the impression of being the superposition of 
two different slopes, one effective between 0 and 
0.6 of the ordinate scale, the other beyond about 
1.2, the intermediate region showing a smooth 
transition. It should be mentioned that other 
curves of Whitehead show a similar behavior. 
This is in agreement with Eq. (5). 

Applying Eq. (5) to Fig. 2, in which case Eo 
was 1500 volts, a=0.1 cm, op=114X10-" 
mho/cm, and o;=36X10-* mho/cm, the follow- 
ing values were obtained: m=3.1X10~’ and 
m' =0.37X10-7 cm/sec. per volt/cm. Table I 
shows the degree of agreement between experi- 
mental and calculated o¢ values. 

Table II shows the results referring to another 
oil, B-6, the figures being somewhat less certain 
in that o; had to be estimated (Point 3 in Fig. 


LOGARITHM OF CONDUCTIVITY VS 
TIME FOR OIL B4 AT 60 C 
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1) in contradistinction to the previous case 
for which the branch 3-4 was horizontal. Here 
a and E, have the same values as before, 
oo=3.24X10-", and o;=0.81 X The values 
1.0X10-7 and 0.31X10-? were obtained for m 
and m’, respectively. The agreement is somewhat 
less satisfactory than in the previous case. 
Figure 3 shows the theoretical curves and experi- 
mental points for both oils considered. The 
larger mobilities, m, probably correspond to the 
positive ion (solvated H-ion?), the smaller ones, 
m’, to the negative ion. 

In computing ionic radii according to Eq. (6), 
n for the cation will be taken as one. The vis- 
cosities are 1.1 poise for the B-4 oil at 60°C, and 
9.0 poise for B-6 oil at 30°C. The corresponding 
cation radius for the first oil is 2.5 10-7’ cm, and 
9.4X10-* cm for the second. 


TasLeE I. J. B. Whitehead’s B-4 Oil at 60°C. 


TIME IN SEC. EXPERIMENTAL ‘ALC, WITH Eq. (5) 
10 108 112 x10-" 
40 oF 61.5 
100 46 44 
200 40.5 40.0 
400 37.5 37.3 
1000 36.0 36.0 


TABLE II. J. B. Whitehead’s B-6 Oil at 30°C. 


TIME IN SEC. witH Eq. (5) 
10 3.67 X 10-5 3.62 X 
20 3.19 3.24 
40 2.32 2.68 
100 1.45 1.73 
200 1.25 1.19 
400 0.96 0.93 


CURVES DRAWN ACCORDING 
TO EQUATION (5) 
\ © EXPERIMENTAL POINTS 
\ 
3 
8 
beg 2 
\ 
TAN 
10 100 1000 


TIME IN SECONDS 


Fic. 3. Comparison of conductivities computed from Eq. (5) 
with experimental values. 


The number of charges carried by the anion is 
unknown; because of association it probably 
increases with decreasing temperature, while at 
higher temperatures it approaches unity. Thus 
for the first oil at 60°C, and with »=1, the value 
is 2.1X10-* cm, whereas for the second oil at 
30°C we can only write 0.3X10-* m cm. 

In computing the concentration of ions from 
Eq. (7) that for the cations (n=1) is 2X10" per 
cc for oil B-4 and. 1.5 X 10" per cc for oil B-6. The 
number of anions is the same with the restriction 
that it contains 1/n as a factor. 

In spite of these uncertainties it is felt that 
some useful information is obtained by the 
method suggested, and that it should be worth 
extending Whitehead’s method to insulating 
liquids of known chemical constitution. 
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Ferromagnetic Anisotropy and the Itinerant Electron Model 


HarRvEY Brooks 
Harvard University, Cambridge, Massachusetts 
(Received October 3, 1940) 


The purpose of the present investigation is to explore the 
possibilities of the “itinerant” or “collective” electron 
picture in the theory of metals to explain the quenching 
of orbital angular momentum in solids and, through the 
introduction of /,s coupling, the phenomenon of ferro- 
magnetic anisotropy in cubic crystals. The approach used 
js that of the Bloch “approximation of tight binding’’ in 
the theory of metals, the exchange energy being treated as 
a Weiss internal field as in the work of Stoner and Slater, 
and the spin-orbit coupling being introduced as a perturba- 
tion. The anisotropy is shown to appear in the fourth 


approximation, and to have the correct order of magnitude 
for iron and nickel. The model also predicts the correct 
sign of K, in nickel and iron, but this prediction is not 
entirely satisfactory because computational difficulties 
prevent the inclusion of all the d-wave functions in the cal- 
culation. A qualitative discussion of the behavior of iron- 
nickel alloys is given. The chief weakness of the model is 
its failure to take account adequately of Russell-Saunders 
coupling within the atom, and the dependence of many of 
its predictions on details of the model which are not very 
well established. 


INTRODUCTION 


CONSIDERABLE amount of experimental 

data has accumulated in recent years on 
the magnetic anisotropy of ferromagnetic cubic 
crystals." Most of the observed effects can be 
interpreted on the assumption that the crystal is 
composed of a number of magnetized domains,’* 
whose energy depends upon the orientation of the 
magnetic intensity in the domain relative to the 
crystalline axes. Thus the energy of a domain 
takes the form 


E=Eot+f(a1, a2, as), (1) 


where a, @2, a3 are the direction cosines of the 
magnetization. Fortunately the form of the func- 
tion f(a1, a2, a3) may be deduced from the sym- 
metry class* of the crystal without any assump- 
tions as to the mechanism responsible for the 
anisotropy. In a cubic crystal, for example, it 
may be shown that the energy must take the 
form 


E=Ey+ Ki 
(2) 


plus higher degree terms which have never been 
found necessary in practice. The magnitude and 
general properties of the coefficients K; and Ke 
are not, of course, determined by the symmetry, 
1R. M. Bozorth, J. App. Phys. 8, 575 (1937); E. C. 
Stoner, Magnetism and Matter (Methuen, 1933). 
*Cf., for example, W. Heisenberg, Zeits. f. Physik 69, 287 
(1931); R. M. Bozorth and E. C. Stoner, reference 1. 
Introduction to Ferromagnetism (McGraw- 


but may be deduced from experiment.! The 
elucidation of the magnitude, sign, and tempera- 
ture and alloy composition dependence of K, and 
K; is a task for the electron theory of metals. The 
Heisenberg exchange coupling, which accounts so 
nicely for the existence and general properties of 
the ferromagnetic state,‘ is of no use in explaining 
the anisotropy since it depends only on the 
relative orientations of neighboring spins and not 
on their absolute orientation with respect to the 
crystal. The earliest such mechanism proposed 
was the magnetic dipole interaction between 
spins. This interaction contains a term depending 
on the orientation of the spins with respect to the 
line joining them, and such terms are naturally 
structure dependent since they depend on the 
actual positions of the spins in the lattice. In 
cubic crystals, however, the magnetic interaction 
vanishes in first approximation because of the 
symmetrical arrangement of nearest neighbors, 
and the second approximation of perturbation 
theory gives much too small a magnitude for the 
anisotropy constants. Bloch and Gentile® sug- 
gested instead that the spin-orbit coupling, 
acting in conjunction with the electrostatic 
coupling of the orbit to the crystal, would provide 
an anisotropy of the proper magnitude, and at 
present this is the generally accepted view. The 
two models are not dissimilar from the mathe- 
matical viewpoint, since, as Van Vleck*® has 


‘ W. Heisenberg, Zeits. f. Physik 49, 619 (1928). 


5 R. Becker, Zeits. f. Physik 62, 253 (1930) ; F. Bloch and 
G. Gentile, Zeits. f. Physik 70, 395 (1931). 
6 J. H. Van Vleck, Phys. Rev. 52, 1178 (1937). 
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indicated, the mechanism of Bloch and Gentile 
gives rise to an apparent dipole-dipole as well as 
quadrupole-quadrupole coupling which is of much 
greater magnitude than the purely magnetic 
interaction. Calculations using this somewhat 
phenomenological substitute for the true mecha- 
nism were undertaken from different viewpoints 
by Van Vleck and van Peype, a pupil of Kramers.’ 

In spite of their success these methods 
are not entirely satisfactory because of their 
semiphenomenological approach, and it therefore 
seems worth while to attempt a calculation with 
a more explicit model, which gives a more 
detailed insight into the origin of the coupling 
terms. Van Vleck’s theory is based on the many- 
electron (Heisenberg) model, which is too compli- 
cated for rigorous solution. At Professor Van 
Vleck’s suggestion, accordingly, we have treated 
ferromagnetic anisotropy from the standpoint of 
the ‘‘itinerant” or ‘‘collective’’ electron model of 
Bloch, Stoner, and Slater,*® in which each electron 
is treated as belonging to the metal as a whole 
and as moving independently of the other elec- 
trons in a self-consistent field determined by the 
potential of the ions and the ‘‘time exposure” 
field of the other electrons. Although this model 
is a far poorer approximation for the tightly 
bound d electrons involved in ferromagnetism 
than the Heisenberg model, it leads to a secular 
equation which is reasonably manageable, and 
which, with the introduction of spin-orbit inter- 
action, gives rise to ferromagnetic anisotropy in 
a straightforward manner. 


QUENCHING OF ORBITAL MOMENTUM 


The whole mechanism of ferromagnetic ani- 
sotropy bears an intimate relation to the 
quenching of the orbital angular momentum in 
solids. Determinations of the gyromagnetic ratio 
for the ferromagnetic metals and alloys give a g 
factor only slightly less than 2, showing that even 
though the electrons involved in the magnetic 
properties are d electrons, their orbital angular 
momentum is not free to orient and precess in an 
external field. This means that there is a powerful 
force within the crystal tending to orient the 


7 W. van Peype, Physica 5, 465 (1938). 


8 F, Bloch, Zeits. f. Physik 57, 545 (1929); J. C. Slater, 
Phys. Rev. 49, 537 (1936); E. C. Stoner, Proc. Roy. Soc. 
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charge cloud around each atom in a definite 
manner with respect to the crystal axes. The 
weak magnetic coupling of the spin to the orbit 
insures that this force will be ‘‘felt” slightly by 
the spin, whose energy is therefore not entirely 
independent of orientation. The most obvioys 
orienting or ‘‘quenching’’ mechanism would be 
the crystalline Stark effect, lifting as it does the 
degeneracy between the various magnetic sub- 
states of a given total angular momentum, but 
unfortunately this is not an adequate mechanism 
for cubic crystals, where the asymmetry of the 
electrostatic field is not sufficient completely to 
destroy the degeneracy with respect to the 
magnetic quantum number. In a cubic field, as 
shown by Bethe,® an atomic d level subdivides 
into two levels, forming the representations I, 
and I; of the cubic group. The I’; level is doubly 
degenerate but magnetically inactive, but the I, 
level is triply degenerate and possesses matrix 
elements of the magnetic moment between its 
substates. If the level is not completely filled in 
the solid state, therefore, we should expect 
contributions to the magnetic moment from the 
orbit as well as the spin, and a g factor con- 
siderably less than 2. In order to remove this 
discrepancy Sommerfeld and Bethe’ suggested 
an additional quenching arising from an antifer- 
romagnetic coupling between orbital angular 
momenta on different atoms, and this is the sort 
of quenching postulated by Van Vleck and van 
Peype. In the itinerant electron picture the 
quenching arises naturally from the hopping of 
electrons between adjacent atoms. Even in the 
tightly bound d band the hopping takes place 
with sufficient frequency to destroy any ori- 
entation in an external magnetic field. In more 
quantum-mechanical language, the situation may 
be described as follows. In the Heisenberg model 
each electron is confined to its own atom, and its 
energy levels are determined by a field having the 
symmetry of the point group about a fixed atom. 
In this field the representations and energy levels 
are degenerate as described above. In the 

®H. Bethe, Ann. d. Physik 3, 137 (1929) ; Zeits. f. Physik 
60, 218 (1930); H. Kramers, Proc. Amst. Acad. 32, 1176 
(1929); 33, 959 (1930); W. G. Penney and R. — 
Phys. Rev. 41, 194 (1932); 42, 166 (1932); }. H. Van 


Vleck, Electric and Magnetic Susceptibilities (Oxford, 1932), 


Section 73. 
1” A. Sommerfeld and H. Bethe, Handbuch der Physik 


(second edition, Berlin, 1933), Vol. 24/2, p. 613. 
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itinerant model, on the other hand, each electron point group alone. The reality of the resultant 
moves in a. field whose symmetry is determined nondegenerate wave functions then precludes the 
by the space group of the whole crystal, including existence of diagonal matrix elements of the 
translations as well as rotations. For the majority magnetic moment, and such nondiagonal ele- 
of wave vectors k the representations are non- ments as there may be join states whose energy 
degenerate, the space-group insuring the removal differences are large compared to the magnetic 
of the degeneracy which may exist under the _ perturbation. 


APPROXIMATE STRUCTURE OF D BANDS 


In order to put these ideas on a quantitative basis it is necessary to say something about the form 
of the d bands in cubic crystals. This is a case where the Bloch ‘‘approximation of tight binding’’" is 
reasonably reliable, since the overlapping of the d shells is small. The energy levels are just those of 
the single atom subject to the crystalline Stark effect, except that each level is broadened out into a 
band, arising from the possibility of leakage of electrons between atoms. When the original atomic 
levels are degenerate, however, the substates give rise to a number of overlapping bands, and con- 
siderable care is necessary in choosing the correct linear combinations of wave functions which 
diagonalize the perturbing energy. 

In setting up the secular equation we confine ourselves to the I; level, although this is a poor 
assumption, in that there are matrix elements of the magnetic moment between the I’; and I; levels 
as well as between the substates of T;. For the purposes of this paper we shall assume that the 
separation between I’; and I’; is large compared to the broadening of the individual levels. The atomic 
wave functions of the I’; level are: 


oi=yaf(r), d2=sxf(r), o3=xyf(r). (3) 


In the approximation of tight binding we start with the linear combinations :"' 


¥n(k) = N-! exp (ial-k)¢,(r—al), (4) 
1 


where a is the lattice constant, N the total number of atoms in a domain, 1 a vector representing the 
lattice points. With these wave functions and a perturbing potential equal to the difference between 
the crystalline self-consistent field and the potential due to the free ion, the secular equation factors 
according to k, but for each value of k there are in general matrix elements of the ‘‘hopping”’ energy 
(perturbing potential) between states of different ». Thus for each value of k, the wave vector, the 
energy is obtained as a root of a third-degree secular equation: 


H,-E Fy. Hy; Hy= S Hy (k)dr, 
Ay, H,—E | =0, (5) 
FT o3 H;-—E Hy = S (k) Ho(k)dr. 


This equation was first obtained by Jones and Mott” for the body-centered cubic lattice. In this 
case the matrix elements are: w76 


H,=H.,=H;=—8A cosétcosyncost, Hi:.=8B sin é sin 7 cos 
H3=8B cos sin n sin ¢, H,3=8B sin cos sin (6) 


"F, Bloch, Zeits. f. Physik 52, 55 (1928); N. F. Mott and H. Jones, Properties of Metals and Alloys (Oxford, 1936), 


Chapter II. 
'®H. Jones and N. F. Mott, Proc. Roy. Soc. Al62, 49 (1937). 


inite 

The 
orbit 
y by 
irely 
vious 
d be 
s the 
sub- 
but 
nism 
f the 
ly to 

the 
d, as 
vides 
is T; 
he r; 
atrix 
n its 
ed in 
Kpect 

> sort 
1 van 
» the 
ng of 
n the 
place 
more 
may 
nodel | 
ig the 
atom. 
| the 
+hlapp, 


912 HARVEY BROOKS 


where &, n, ¢ are defined by §=}k,a, »=}k,a, {=4$k,a and where 
A=— y+4a, s+4a)dV, 
B= — y+}, 
Since the perturbing potential is negative, A and B are positive." 
For the face-centered cubic the matrix elements are readily obtained. The diagonal elements are: 
H,=4A;(cos & cos 7+cos cos {) cos 9 cos 
H,=4A;(cos £ cos n+cos cos cos cos 
H;=4A;,(cos cos ¢+cos cos ¢) —4A>2 cos cos n, (8) 
— Sxyf(r)H' (x+ 43a) (y+4a)f(| (x+4$a)?+ 
A simple argument shows that A; and A; are positive as defined and that A2>A,." 
The off-diagonal elements are 
His=—4B sin ésin (9) 


and B as defined is positive. 
SOLUTIONS OF SECULAR EQUATION 


1. Body-centered cubic 
The solution has been obtained in approximate form by Jones and Mott." The result is expressed 
in terms of a variable x= (E/8B)+(A/B) cos = cos n cos ¢. For most values of the wave vector the 


approximate solutions are: 
x=0, (10) 


where 
a=sin ésinyncos{, c=cos ésin ¢. 

With these energies the approximate form of the diagonalizing transformation is: 
$1 


0; be ae" 
+) 4v2re* (11) 
—|—43v2re" = 


where 


This gives the energy correctly except for terms of order abc/e*, which happens to be small for most 
values of the wave vector k. Thus the continuum generated by the three I’; functions really consists 
of three overlapping bands. Each band is referred to an atomic function which is a linear combination 
of the functions (3), and which is different for each choice of propagation vector. 


2. Face-centered cubic 
We divide the secular equation through by 4B and define 


«= (E/4B) —(A,/B)(cos cos »+cos cos §+cos cos), Q=(Ai+A:2)/B. (12) 


Here it must be remembered that B <(Ai+ Az), and that therefore Q is large compared to unity. This 
being the case, the secular equation is largely dominated by its diagonal elements, except for the 
special values of k for which they become nearly equal. This means that in the face-centered cubic 
the original ¢1, ¢2, ¢3 are nearly the true atomic wave functions for the three bands of the I’; level. 

Similar results apply in the simple cubic lattice, but here it may be shown that the off-diagonal 


matrix elements vanish altogether so that the functions (3) are just correct. 


18 This follows from the fact that xyf(r), yaf(r), and sxf(r) are even with respect to reflections along the body dia , 
\ A, is positive since the perturbing potential is negative while x and x+ $a have opposite signs on nearest neig bors. 
A; is positive because the integrand is positive in the region between nearest neighbors. A; >A, because y~4, i.e., the 


wave function has a node in the region of maximum perturbing potential. 
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Sprn-OrBIT INTERACTION 


In order to carry on the discussion it is necessary to obtain the matrix elements of the spin-orbit 
jnteraction in the system of representation in which the band or “hopping”’ energy is diagonal. This 
interaction may be expressed in the form 


(13) 
1 


where L(r—al) is the angular momentum operator about the atom at lattice position 1, and ¢(|r—al|) 
js a function falling off as 1/r* at large distances and depending only on the absolute magnitude of the 
radius vector.'® O’ is an operator which is periodic in the fundamental lattice translations, and hence 
is rigorously diagonal in k.'* Its matrix elements between Bloch functions of the same k and different 
are furthermore independent of k to a first approximation. It is easily shown that 


(nk| §(|r—al|)S-L(r—al) | n’k) = exp (—tal-k) (14) 


If ¢(r) falls off sufficiently rapidly the integral in (14) will be inappreciable unless 1=0.!7 Thus we see 
that the matrix elements of the spin-orbit interaction with respect to Bloch functions are practically 
the same as between the corresponding atomic functions, that is 


S V =A (n'm,'|S-L|nm,). (15) 


In this equation A is the ordinary spin-orbit interaction parameter for free atoms. It is the interaction 
for an individual electron, corresponding to the small a; in Goudsmit’s'* notation, and must not be 
confused with thé spin-orbit interaction parameter for the atom as a whole. The point is that we are 
here assuming the “hopping” energy sufficient completely to destroy the Russell-Saunders coupling 
within the atom, so that there is no question of the individual electronic angular momenta combining 
to form a resultant L and S. This is a poor approximation since, owing to the tight binding of the d 
shell, there must be a rather high degree of internal coupling, but so far no way has been developed 
of taking it into account in the itinerant electron model, where the electrons are treated as essentially 
independent. We assume that the secular equation for the I’; bands is diagonalized by a set of real 
transformation coefficients a(r, m). In terms of this diagonal representation (labeled by the index r), 
we find for the matrix elements of the spin-orbit interaction 


(rm,|S-L| r’m,’) => a(r, n)(nm,|S-L|n'm,’)a(r’, n’) 
= {a(r, n)a(r’, n’) —a(r’, n)a(r, n’)}(nm,|S-L|n'm,’), (16) 


n<n’ 


where the index n refers to the original set of atomic wave functions [Eq. (3) ]. In this representation 
the spin is quantized along the z axis. In practice, however, the spin must be quantized parallel to the 
external field. In this new system of representation the matrix elements are as follows: 


(r1, $|S-L| ro, = —i cos 6’, (71, $|S-L| 73, $) sin sin ¢’, 

—}) =exp (—i¢’)[—cos ¢’—i cos 6’ sin ¢’], 
(ro, $|S-L| 73, —$) =exp (—i¢’)[—sin $’+ cos 6’ cos ¢’]. 


The angles 6’, ¢’ are measured not with respect to the x, y, z axes but with respect to a new set of axes 
determined by the unit vectors i’, j’, k’, which are defined in terms of the transformation coefficients 
as follows: 


i= [a(r:, 1), a(n, 2), a(n, 3)], = [a(re, 1), a(r2, 2), a(r2, 3)], k’= [a(rs, 1), a(rs, 2), a(rs, 3)]. (18) 


4% E. U. Condon and G. Shortley, Theory of Atomic Spectra (Cambridge, 1936), p. 120. 

*N. F. Mott and H. Jones, reference 11, p. 59. ° 

"A more detailed analysis of the spine: effects would necessitate inclusion of the spin-other-orbit interaction. 
For atoms distant from the one in question this would tend to cancel out almost exactly the 1/r? term arising from the 
ions, and so in effect ¢(r) would oecend fall off exponentially rather than as 1/r*. 

*S. Goudsmit, Phys. Rev. 31, 946 (1928). 
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914 HARVEY BROOKS 
TABLE I. Secular equation for Pomme anisotropy, including spin-orbit interaction. Angles measured with respect to 
special skew axes determined by transformation coefficients, as shown in Eq. (18). 

72,4 13,4 m1, 72, 73, 
Tl, —iA cos0’ iA sin @’ sin ¢’ 0 —iA sing’ —Acosp’—iA cosé’ sin 
72, iA cos 6’ E:+5—E —iA sin cos ¢’ iA sin 0’ 0 —Asing’+i 
—iA sin 6’ sin ¢’ iA sin cos ¢’ ‘Acos¢’-+iA cosé’sing’ Asing’—iAcosé’cos¢’ 
—iA sin’ Acos¢’—iAcosé’sing’ iA cosé’ —iA sin sin ¢’ 
72,—$| tA siné’ Asin ¢’+iA cos6’ |—iA E:—3—E iA sin 
73, —} |—A cos¢’+iA cos@’ sing’ —A sing’—iA cos¢’ 0 +i. ‘sin 6” sin ¢’ sin 0’ cos ¢’ 


Because of the unitary property of the transformation matrix the vectors i’, j’, k’ form an orthogonal 
unit set, and 6’, 9’ are the polar angles of the magnetic field vector measured with respect to this set. 

The complete secular problem, including the spin-orbit interaction, is of the sixth order, corre- 
sponding to the three overlapping bands of the I's level, and the two possible orientations of the spin. 
The secular determinant is shown in Table I. E,, E2, and E; are the energies of the three bands, and 
are, of course, functions of the wave vector k. In this problem we have made the usual approximation 
of replacing the exchange energy by an internal ‘“‘molecular” field 6=kr’d. Here X is the fractional 
magnetization I/Jo, where Jo is the saturation intensity, and kr’ is mgJ, where mz, is the number of free 
spins per atom (electrons or holes) corresponding to the saturation magnetization; also J is the 
atomic exchange integral,’ that is, the exchange integral between electronic wave functions normalized 
per unit cell rather than for the whole crystal. In the system of representation we have chosen, the 
energy matrix is diagonal in everything except the spin-orbit interaction, which may be considered 
as small compared to the difference between diagonal elements in Table I. Thus the solutions of the 
secular problem may readily be expanded in powers of the spin-orbit parameter A, the successive 
powers being obtained in the successive orders of perturbation theory. 


DETAILED DERIVATION OF G FACTOR where ¢1, $2, ¢3 are the wave functions for the 
three bands and the indices +, — stand for the 
spin directions. We must find the mean of the 
magnetic energy with respect to these wave 
functions. Unfortunately, since it is linear in 
the coordinates, the vector potential cannot be 
treated as a small perturbation throughout the 
crystal. This difficulty may be obviated by a 
method of gauge transformations due to Peierls 
and F. London,” in such a way that the magnetic 
energy for Bloch functions is the same as for 
the corresponding atomic functions.” The cor- 
rection to the energy arising in this way is 


We have already seen that the removal of the 
degeneracy of the cubic levels under the space 
group accounts for the almost complete quench- 
ing of the orbital angular momentum. There is, 
however, a secondary effect which allows the g 
factor to differ slightly from 2. Owing to the spin- 
orbit interaction there is in effect a small internal 
field tending to orient the orbital angular mo- 
mentum, in addition to the large field acting on - 
the spins and arising from exchange. The result 
of this field is always to make the g factor 
slightly less than 2. From the secular equation 
we obtain the first-order correction to the wave- Baie ee as 
functions by the usual perturbation method : 
‘Here a’, a2’ are direction cosines with respect to 


iA sin 6’ sin ¢’ iA sin 6’ cos ¢’ 

or + the skew axes defined by (18). Equation (20) 
a i, gives the energy correction for one electron. To 
cos ¢’+iA cos sing’ obtain the magnetic energy for the crystal as a 
+28 (19) whole we must sum over the occupied electronic 
nag states of the Fermi distribution, taking proper 
Asin ¢'—14 cos 6 cos ¢ a account of the symmetry properties of the energy 
ok. Peierls, Zeits. 763 (1933); F. London, 

S One’ neturelly neglects diamagnetic effects arising 


19 J. C. Slater, E. C. Stoner, reference 8. from the migration of electrons among many atoms. 
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surfaces.” Finally dividing by the total angular 
momentum leads to a g factor 


g=24+-9A(e (21) 


The average is taken over all the occupied states 
for both directions of spin. Since the second 
term of (21) averages to zero when all the d 
levels are fully occupied, we may replace the 
average over the occupied states by the negative 
of the average over unoccupied states or ‘‘holes’’ 
of the distribution. In nickel and iron the band 
is nearly full, and most of the holes are in the 
highest of the three energy bands so that the 
energy differences ¢:(=E;—E;) and e.(=E;— 
are positive for the overwhelming majority of the 
states involved in the averaging processes. Thus 
the second term of (21) is always negative, and 
the g factor is slightly less than 2, since A is 
positive for single electrons. This is in agreement 
with observation, and means simply that there 
is a slight tendency for the orbit to aline itself 
parallel to the .orienting fields (external and 
internal). Assuming a spin-orbit parameter of 
about 500 cm, and a mean band width of 
about 20,000 cm, we obtain a g factor in the 
neighborhood of 1.97. Experimentally the de- 
partures from 2 range between 0.1 given by 


A‘tj1 1 1 1 2 1 1 1 1 1 
+ +44 ———— + \\—- 
eot+25) Ler? (e1+25)? (€2+26)? 


€2 +25 €2+26 


€2 


1 
Au=A—-— 


1 1 : .% 1 1 1 

€:+26 26 € lee 
where now 6=ngJI/Ip. 


In Eq. (22) we must expand the primed direction 
cosines in terms of the unprimed ones, using the 
definition of the primed axes: 2) an, 
where the a, are defined so that a,;= az, a2= ay, 
a3=a,. Substituting into Eq. (22), and using the 
symmetry argument once again, as well as the 
unitary property of the transformation coeff- 
cients, we obtain a relatively simple formula for 


* Brouckaert, Wigner and Smoluchowski, Phys. Rev. 
50, 58 (1936) ; Conyers Herring, Phys. Rev. 52, 365 (1937). 
Cf. particularly Manning and Chodorow, Phys. Rev. 56, 
787 (1939), discussion in fine print at end of paper. 


Barnett”* to 0.004 obtained by Sucksmith and 
Bates for Ni.** Our calculated departure of 0.03 
is thus of the correct order of magnitude. 


CALCULATION OF ANISOTROPY CONSTANT 


The calculation of the anisotropy constant K, 
is made by applying fourth-order perturbation 
theory to the secular equation shown in Table I. 
The energy corrections arising from lower orders 
of perturbation theory are easily shown to 
possess no angular dependence when the summa- 
tion is carried out over all the occupied states. 
This result follows from the symmetry properties 
of the energy surfaces and is analogous to the 
argument usually employed to fix the form of 
the angular dependence in cubic crystals.*° In 
carrying out the fourth-order calculation, we 
need include only terms whose angular de- 
pendence is of the form a;’a;’?.2* In this way we 
obtain for the fourth-order energy an expression 


A + A +A (22) 


in which a’, a2’, a3’ are the direction cosines 
of the magnetic field with respect to the skew 
axes defined by (18), and Ax, Ais, Aos are given 
by the following expressions: 


1 1 1 


1 1 1 1 1 
€2 €2+26 26 €1 +26 €1 €: +26 


€2 +26 


2 
(23) 


the anisotropy constant K,, namely 
Ki=}3 A,,[1 —5J,»] 


us 
X (24) 


23S. J. Barnett, Rev. Mod. Phys. 7, 129 (1935). 
** W. Sucksmith and L. F. Bates, Proc. Roy. Soc. A104, 
L. F. Bates, Modern Magnetism (Cambridge, 


). 
25 Details of this argument, which is really quite obvious 
pee” are included in the author's thesis, Harvard 


26 It is to be understood, however, that all terms in the 
fourth power of the direction cosines such as a,‘ are to 
be expressed as: a3"). 
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where the J,, are defined by a biquadratic ex- 
pression in the transformation coefficients: 


(25) 


On account of the terms in J,, the anisotropy 
bears a very close relation to the structure of 
the energy bands. If the original wave functions 
xyf(r), y2f(r), 2xf(r) are a good set of atomic 
functions for setting up the Bloch waves which 
diagonalize the band energy, then the J,, will be 
small and the anisotropy, as will be shown, is 
negative. But if the bands are not well described 
by the functions xyf(r), etc., but rather by linear 
combinations such as ¥,= > ,a(r, )¢n, then the 
Jy» May become important, and the anisotropy 
may even change sign and become positive as we 
shall show to be the case for iron. This effect has 
no counterpart in the theories of Van Vleck and 
van Peype, but seems to be characteristic of the 
particular model we have adopted. 


NICKEL 


The above discussion is well illustrated by the 
two common ferromagnetics nickel and iron. 
As can be seen by reference to the section on the 
structure of the d bands, nickel is a good example 
of the case where the functions xyf(r), etc., are a 
good approximation and the J,, are small. To 
simplify the discussion we shall consider only the 
case of magnetic saturation, i.e., very low tem- 
peratures. The expressions for the A,, in Eq. (23) 
must be summed over all the occupied states in 
the three bands. The summation is carried out 
most readily by expressing each of the A,, as a 
sum f:(5)+/2(5) of an even function of 6 plus an 
odd function of 6, respectively (f1(5)=f:(—§4), 
f2(5) = —f2(—6)). Here 5, it must be remembered, 
is a quantity proportional to the exchange 
energy, whose sign is negative for spins parallel 
to the external field and positive for spins anti- 
parallel to this field. The terms in the summation 
of the odd function cancel out in pairs for all 
states which are occupied by both directions of 
spin, so that we are left with a summation over 
only those states which contribute directly to the 
magnetic moment, whose spins are parallel to 
the external field. On the other hand, the summa- 
tion for the even function must be carried out 


BROOKS 


over all occupied states since the function has 
the same sign for both spin directions. However, 
>/f:(6) vanishes when taken over all states, un- 
occupied and occupied alike, as follows by 
application of the diagonal sum rule to the 
matrix of Table I. Thus the sum over occupied 
states may be replaced by the negative of the 
sum over unoccupied states, and since the sign 
of 6 is immaterial for the function f1(5), the 
summation may be carried out over exactly the 
same states as fz, namely, those contributing 
directly to the magnetic moment. Thus we have 
the result: K,;=sum over all occupied states of 
(A12+A23+A13)=sum over states of unpaired 
spin of —f:(5)+f2(—6)= The 
last is just the negative of Eqs. (23) summed 
over all the magnetic states, but with § 
treated as a positive quantity. Since we are 
dealing with a band which is well over half full, 
most of the unpaired spins lie in the uppermost 
of the three I’; bands, so that both the quantities 
€, and ¢2 are positive. With these facts it is 
readily shown that A,, in Eq. (23) are all posi- 
tive, irrespective of the numerical magnitudes of 
€:, €2, and 6. Thus the sign of K, in nickel is 
proved to be negative under quite general con- 
ditions. 

The most recent results on the anisotropy of 
iron and nickel are due, for nickel, to Bozorth 
and Williams?’ who find K,=8X 10° ergs/cc by 
extrapolation to absolute zero; for iron to L. P. 
Tarasov"* who gives K,=5.3 X 10° ergs/cc at room 
temperature. The latter value would probably 
not be increased by a factor of more than 2 on 
extrapolation to absolute zero. It is of course 
impossible to reproduce these values to within 
better than a factor of 10 by means of the theory, 
but with this accuracy the theory predicts the 
correct order of magnitude. In his paper on 
ferromagnetism Slater®® has given curves showing 
the distribution of states in the d band for 
copper, which he assumes to be similar to that of 
nickel. The six highest energy electrons, which 
presumably correspond to those occupying the 
states of the I’; level in our model, occupy an 
energy range of about 50,000 cm~. Thus we 


asst J. Williams and R. M. Bozorth, Phys. Rev. 56, 837 
28. P. Tarasov, Phys. Rev. 56, 1231 (1939). 
29 J. C. Slater, Phys. Rev. 49, 537 (1936). 
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should not be far off in taking as a mean value 
for the separation of levels of the same k, (i.e., 
the quantities ¢, and ¢, in Eq. (23)), a value of 
10,000 cm. Slater has also estimated the value 
of the exchange integral 5, assuming that it 
arises from intra-atomic exchange. On the face 
of it, this seems highly unlikely, but the work 
of Stoner,” on the other hand, certainly indicates 
that the exchange energy cannot be much less 
than the band spread if ferromagnetism is to 
occur at all. We should probably not be far off 
in taking 25 equal to 8000 cm™ for nickel. 
Estimates of the spin-orbit parameter A have 
been made for the transition elements by 
Goudsmit.'® Again these estimates cannot be 
taken too seriously in the solid state, not only 
because of the uncertainty arising from the 
incomplete destruction of Russell-Saunders coup- 
ling, but also because of the entirely altered 
screening effect of the conduction electrons in the 
solid phase.** With Goudsmit’s estimate of 
A=630 cm for nickel, we obtain a value of K, 
between 10° and 10° ergs/cc, which is the correct 
order of magnitude. The screening effect of the 
conduction electrons would tend to reduce the 
estimate of the spin-orbit parameter somewhat, 
and provide worse agreement with experiment. 
For iron Goudsmit estimates A = 390 cm-, which 
would give an anisotropy for iron less than } that 
of nickel, other factors being equal. However, the 
order of magnitude is still correct, and most of 
the discrepancy can be accounted for by the 
greater exchange energy in iron, and the greater 
number of spins per cc. 


oF K,—IRON 


Iron is a good example of the case where the 
Jy of Eqs. (24) and (25) are not negligible. 
This arises from the fact in the body-centered 
cubic the diagonal elements of the secular equa- 
tion for the “hopping” energy are all the same. 
For the highest of the three I’; bands, namely the 
band corresponding to x= + in Eq. (10), the J,, 


may be obtained from the transformation matrix 


*E. C. Stoner, Proc. Roy. Soc. A165, 372 (1938). 
* This effect is discussed in detail by M. F. Manning and 
M. I. Chodorov, Phys. Rev. 56, 787 (1939). 


(11) by means of (25), and are as follows: 
Jo=R, Jis=}(1—R), 
R= /e4, (26) 
e= (a?+b?+c?)!. 


a, b, c are defined in terms of the components of 
the wave vector, as shown in Eq. (10). The 
maximum value of R is }. 

With these values of J,,, the anisotropy con- 
stant follows from Eq. (24). The exact expression 
is very complicated, but it may be very closely 
approximated by the formula: 


45° 1 
3« 


(27) 


This expression applies only to the uppermost of 
the three bands, in which most of the electrons 
will lie in practice. The expression is easily seen 
to be positive for all positive values of ¢ and 4, in 
agreement with the observed positive sign of K, 
in iron. 

The theory is very attractive in that it makes a 
satisfactory unambiguous prediction of the cor- 
rect sign of K, in nickel and iron. When we come 
to examine the data on alloys, however, some 
difficulties begin to appear. The theory suggests, 
for example, that the sign of the anisotropy is 
to be correlated with the crystal structure, the 
difference between iron and nickel arising from 
the difference in the form of the energy bands 
between the face-centered and the body-centered 
lattices. Yet the observations on the Fe-Ni alloy 
system® show no marked discontinuity in the 
anisotropy as we pass from the face-centered to 
the body-centered phase. In fact the change in 


sign of K, takes place well within the y-phase. 


(face-centered) at a composition corresponding 
to one free spin per atom. All the alloy data 
suggest a much closer correlation of the aniso- 
tropy with the number of free spins per atom 
than with the crystal structure. In spite of this 
it is possible to devise a crude qualitative ex- 
planation of the alloy data, based on the itinerant 
electron model, somewhat as follows. In our dis- 
cussion of the band structure of nickel we 


assumed that the off-diagonal matrix elements of 


® R. M. Bozorth, reference 1. 
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the ‘“hopping’’ energy could be neglected. This 
is true in the approximation of very tight 
binding, but at the lattice spacing which actually 
obtains in nickel it is probably a very poor 
assumption except in portions of the Brillouin 
Zone in which sin £, sin 7, and sin ¢ are very 
small, where, as a consequence, the off-diagonal 
matrix elements are negligible. The electrons or 
holes in the highest energy states occupy just 
such a portion of k-space. Now, nickel has only 
0.6 free spin per atom, and most of these lie in 
states very near the maximum. We should there- 
fore expect the anisotropy of nickel to be nega- 
tive, since small off-diagonal elements of the 
“hopping” energy imply small J,, in Eq. (24) 
and hence negative K;. As iron is added the 
number of free spins per atom increases rapidly. 
The additional holes are forced into parts of 
k-space where the off-diagonal ‘‘hopping”’ in- 
tegrals are large, and these holes will therefore 
tend to have positive K,, opposing the effect of 
holes higher in the band, until finally the interior 
holes predominate and cause the K, for the whole 
crystal to change sign. For the body-centered 
phase a somewhat similar explanation applies in 
reverse. Here the experiments indicate a gradual 
falling off of K, with increasing percentage of 
nickel. In our model of the body-centered 
structure we have neglected the effect of next 
nearest neighbors. This effect will indeed be 
small except when the off-diagonal elements of 
the hopping energy are small, namely, near the 
top of the band. In this case the secular equation 
for the hopping energy takes a form similar to 
that for the simple cubic, since the next nearest 
neighbors lie in the same relative positions as 
the nearest neighbors in the simple cubic. Thus 
holes very near the top of the band in iron may 
be expected to make a negative contribution to 
K,, since K, for a simple cubic may easily be 
shown to be negative in our model. Now the 
addition of nickel to iron decreases the number 


BROOKS 


of free holes, and so increases the relative number 
of states of negative anisotropy, as these lie 
closest to the top of the band and so tend to 
hang on to the bitter end, as it were. There will 
thus be a gradual decrease of the positive 
anisotropy of iron with increasing percentage of 
nickel, in agreement with experiment. 

This explanation of the alloy data is very 
rough, and an attempt at quantitative verifica- 
tion of it with the present model did not prove 
successful. Nevertheless, the effects calculated 
were in the right direction, if not of sufficient 
magnitude to substantiate the qualitative dis- 
cussion. This is hardly surprising in view of the 
very approximate nature of our treatment of 
the band structure. 


It would be desirable to obtain some idea of 


the temperature variation of K,. Unfortunately, 
this would involve integration over the distribu- 
tion with a Fermi distribution function, a task 
in which the labor involved is hardly com- 
mensurate with the results obtainable from so 
crude a picture of the band structure. It was at 
first believed that the temperature variation 
could be represented by a formula such as (27), 
the anisotropy being expressed as a function of 
magnetization as in Van Vleck’s theory. The 
variation obtained from this formula, however, 
is much too slow, being of the order of the 
square of the magnetization instead of the tenth 
power as observed experimentally for iron.** The 
agreement with experiment in this respect is 
considerably worse than in Van Vleck’s theory. 
Probably most of the discrepancy can be taken 
care of under the assumption of variations of the 
J,» in Eq. (24) arising from thermal excitation 
of the band electrons. 

In conclusion I should like to express my 
thanks to Professor Van Vleck for suggesting the 
problem and for many helpful discussions and 
suggestions during the course of the work. 


3 N. Akulov, Zeits. f. Physik 100, 197 (1936). 
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The considerations of the first paper of this series are modified so as to be consistent with the 
special theory of relativity. It is shown that the inertia of energy does not obviate the necessity 
for assuming the conservation of matter. Matter is to be interpreted as number of molecules, 
therefore, and not as inertia. Its velocity vector serves to define local proper-time axes, and 
the energy momentum tensor is resolved into proper-time and -space components. It is shown 
that the first law of thermodynamics is a scalar equation, and not the fourth component of the 
energy-momentum principle. Temperature and entropy also prove to be scalars. Simple 
relativistic generalizations of Fourier’s law of heat conduction, and of the laws of viscosity are 
obtained from the requirements of the second law. The same considerations lead directly to 


the accepted relativistic form of Ohm’s law. 


INTRODUCTION 


N the second paper of this series,' the theory 

of e-substitutions was outlined, and it was 
shown that this device can be used to simplify 
the derivation of some thermodynamic formulae. 
However, the author was reluctant to use it in 
the derivation of any fundamental formulae 
because the e-substitution depends on the fact 
that the internal energy contains an arbitrary 
additive constant, while the principle of the 
equivalence of mass and energy appears to give 
an absolute meaning to the energy content of 
any portion of matter. Thus it seemed possible 
that in a relativistic theory of thermodynamics, 
«substitutions might be inadmissible. The con- 
siderations Which follow show that this is not 
the case. To establish this result, it was necessary 
to develop a complete relativistically invariant 
theory of the simple fluid—a project which is of 
interest in itself, even though the departures of 
such a theory from a classical theory must be 
negligible in most cases. 

In I, four general principles were used: 
(a) The conservation of matter; (b) the con- 
servation of momentum; (c) the conservation of 
energy; and (d) Kelvin’s hypothesis concerning 
the thermodynamic temperature. It might seem 
as though the equivalence of mass and energy 
would make (a) superfluous in a relativistic 
theory. This is not the case, as becomes clear if 
matter matter is interpreted as number of molecules 


“10, Eck Eckart, Phys. Rev. 58, 267, 269 (1940); hereafter 
cited as I and. II, respectively. 


rather than inertia. The principles (b) and (c) 
combine into a single tensor equation, as is 
well known. This is somewhat disconcerting, for 
the first law of thermodynamics is a scalar 
equation ; its relation to the energy-momentum 
principle must be discovered. Moreover, the 
correct form of the energy-momentum tensor is 
still a matter of discussion, and some assumption 
must be made before the theory can be com- 
pleted. Finally, Kelvin’s hypothesis involves 
Fourier’s law of heat conduction, and the 
relativistic generalization of this presents a 
difficulty, arising from the fact that in the 
classical theory there is a three-vector that 
represents the flow of heat, but no density of 
heat to combine with it to form a four-vector. 
A somewhat similar difficulty arises in connection 
with the stress-tensor.? It will be shown that all 
of these problems can be solved in a systematic 
manner if one defines proper-time and proper- 
space in terms of the velocity of matter. 


MATTER, PROPER-TIME, AND PROPER-SPACE 


A Galilean coordinate system (x°=ct, x!, x*, x*) 
will be assumed, in which the metric tensor is 


010 | (1) 
001 


? R. C. Tolman has solved this problem for the stresses; 
see his Theory of the Relativity of Motion (University of 
California Press, 1918), Chap. X, and Renan 
dynamics and Cosmology (' ord, 1934), Cha 
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and this tensor will be used to raise and lower 
indices in the usual manner. Matter will be 
represented by a four-vector, m*, which has the 
units g/cm*, and depends only on the molecular 
weight and on the motion of the molecules. The 
inertia and internal energy of the molecules will 
be introduced later, so that the molecular weight 
is to be considered only as a convenient con- 
version factor. The conservation of matter is 
expressed by the equation 


Om*/dx*=0. (2) 

It is to be assumed that m°>0, and that the 
invariant density 

m =(—gagmem*)* (3) 


is real and positive. The (dimensionless) velocity 
of matter is 
u*=m*/m, (4) 
so that 
Ugu*= —1 (S) 
and 
= uq(du*/dx*) =0. (6) 
The vector u* will be a function of the co- 
_ ordinates, and will determine a direction at each 
point of four-space—the local axis of proper- 
time. Three directions orthogonal to this are the 
local axes of proper space. Any vector F*, then 
gives rise to a scalar and a vector by means of 
the equations 


f=—uaF*, (7) 
fax (8) 
= + (9) 


The scalar f is the projection of F* on the axis of 
proper-time, and f* is the projection of F* into 
proper-space. It is easily seen that 
Fe=fur+f« 
and that 
UeSp*=0, (10) 
Covariant vectors can be resolved into proper 
components in an analogous manner. 
The tensor sg* has several important proper- 
ties: it is easily seen that 
Sp%S (11) 


and 


sf = 
=g%+ (12) 
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is a symmetric tensor. A simple calculation shows 
that 


so that, because of Eq. (5), 
0, (14) 


the equality holding only when F,=fua. 

This resolution of vectors into time and space 
components is quite consistent with Einstein’s 
theory of simultaneity, since there is no reference 
to separated points of four-space. This becomes 
clearer on defining the proper-rate of change of 
a quantity: 

(15) 
This may be called the derivative of ¢ with 
respect to proper-time. However, proper-time 
itself has not yet been defined; one may try to 
define it by the equation 


Dr=1 


but this does not determine a unique function 7, 
To make 7 unique, it is necessary to specify one 
surface in four-space on which r=0, i.e., one 
locus of simultaneous events. The essence of 
Einstein’s theory is that this locus cannot be 
specified in a physically unique manner, so that 
Tt remains somewhat arbitrary. The differential 
operator D is unique, however, and corresponds 
closely to the classical operator D/Dt. 
The Eq. (2) is readily transformed into 


du*/dx*= —vDm=mD2, (16) 
where v=1/m is the invariant specific volume 
(cf. I, Eq. (2)). Equation (6) yields 

u*Du,=u,Du*=0. (17) 
RESOLUTION OF THE ENERGY-MOMENTUM 
TENSOR INTO PROPER COMPONENTS 


The energy-momentum principle will be as- 
sumed in the usual relativistic form, 


OW*/ax*=0, = Wee, (18) 
the units of W being erg/cm*. 


The equation obtained by setting 6=0 in 


Eq. (18) is variously called the equation of 
conservation of matter, of mass or inertia, and 
of energy. It seems to the author that, while 
each of these designations can be justified in 


cases, the first two are misleading in 
general. He would prefer to retain only the third, 
and refer to the equation W*°=W* as the 
principle of the inertia of energy. Even the third 
js misleading unless a distinction is made between 
the principle of conservation of energy and the 
first law of thermodynamics, the two being 
jdentical only when matter is at rest (see below). 
For the moment, no assumption will be made 
concerning the form of W, but it will be used 
in defining other quantities, such as the internal 
energy and heat flow. Assumptions concerning 
the form of these quantities will be made below. 
The tensor W may be resolved into proper- 
components: if 


w= (19) 
w* = (20) 
wt =s, 2s PW, (21) 


then 


It will appear that each of these components 
has an important physical interpretation. Thus, 
w is the invariant density of energy (erg/cm*). 
The internal energy, « (erg/g=cm?*/sec.’), is 
appropriately defined by 


m(e+a)=w, (23) 


a being an arbitrary constant. It will appear 
that a does not enter into the first law of thermo- 
dynamics. The vector 


g* =cw* (24) 
is the heat flow in erg/cm? sec. ; the equation 
=0, (25) 


which follows from Eq. (10) and (20), corre- 
sponds to the absence of a density of heat in the 
classical theory. The symmetric tensor —w* is 
the stress tensor (dynes/cm?) ; the equation 


=0 (26) 
reduces the number of its independent compo- 


_ nents to six, as in the classical theory. 


THE Finst Law OF THERMODYNAMICS 


The first law of thermodynamics can be derived 
from Eq. (18) and the definitions just adopted. 
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Multiplying Eq. (18) scalarly by —ws and 
rearranging, we have 


(upW**) + (27) 
because of Eqs. (19), (20), (22) 


— = wut+we 
=m(e+a)u*+q*/c. 


Since mu*=m-, Eq. (2) results in the disappear- 
ance of a from the equation :* 


(8/ax*)(upWe8) = 
=mDe-+ (1/c)(dq*/dx*). 


Equations (6) and (22) result in 
(dug/dx*) = 
so that finally Eq. (27) becomes 


mDe-+ (1/c)[(dq%/dx*) +q*Due 
+w*(dug/dx*)=0. (28) 


This is to be compared with the classical form 
of the first law (cf. I, Eq. (7)) 


m(De/Dt)+V-q—(p-V)-V=0, 


q being the heat flow, V the velocity and p the 
total stress. The terms of the two equations are 
in an obvious correspondence, except that the 
term (1/c)g*Du. does not appear in the classical 
case. This term is very small in all ordinary 
cases, and may be interpreted as the work done 
by a flow of heat through accelerated matter— 
a phenomenon not envisaged by the classical 
theory. 


* It will be noted that, up to this point, the introduction 
of the matter vector, m*, served no purpose t to 
define the density and velocity of matter, m oa ue, 
Tolman, reference 2, and others avoid introducing m*, but 
cannot avoid introducing u*. At this place in the considera- 
tions, however, it becomes important that matter be 
conserved (Eq. (2)). If this were not the case, the constant 
a would enter into the first law of thermodynamics, Eq. 
(28), in the form of a term a(4m*/dx*). In other words, the 
creation or annihilation of matter requires energy ab- 
sorption or release. This is as it should be, and indicates 
that the theory of e¢-substitutions may be used even in 
relativistic considerations, except when matter is being 
created or destroyed. The theory of such nuclear processes 
is beyond the scope of the present paper, ge | 
insofar as electric charges and radiation are involved. 
Insofar as only heavy particles are involved, the device of 
using mass numbers (integers) in the definition of m* (and 
thus including the packing-fraction energy in ¢) will insure 
the validity of Eq. (2). 
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THE SIMPLE FLUID AND THE SECOND LAW 
OF THERMODYNAMICS 


The theory of the simple fluid is of considerable 
importance, since in this case it is possible to 
prove the existence of temperature and entropy 
(cf. I, Eq. (9) et seg.). The appropriate definition 
of the hydrostatic pressure is 


p= jw." (29) 


and it is to be assumed that the viscous stress 
tensor, 

—w*+ (30) 
is a linear function of (du,/dx*) and proportional 
to the scalar coefficient of viscosity, \. This 
together with the equations 


UeP®=0, P.*=0, (31) 


which follow from Eqs. (29) and (30), restricts 
the form of P very materially. The equation 
(du,/dx*) + (us/dx) 

— (32) 
is very strongly indicated, if not uniquely 
determined. At any rate, this expression for P 
satisfies Eq. (31), and differs from the classical 
expression only by terms that are ordinarily 
negligible (see below). 

Since 
= (du*/dx*) =mDv 

by Eggs. (6), (12), and (16), the Eq. (28) becomes 


m(De+pDv) + (1/c)[(dq*/dx*) +9*Dua 
— P8(dug/dx*)=0. (33) 
For the simple fluid, ¢ is a function of p and 
v only, so that there are two functions 6 and » 


such that 
De+pDv=6Dn; (34) 


the proof is identical with that given in I. Since 
¢, p and v are scalars, 6 and 7 will also be scalars.‘ 
Substituting in Eq. (33) and rearranging, we 
obtain 
mDn + (d/dx*) (q*/cd) = — (1/c#*)g*[ (90/dx*) 

+ (1/0) P%*(dug/dx*). (35) 


‘The scalar character of » appears to be generally 


accepted (cf. Tolman, Relativity of Motion, Chap. XI). 
However, M. Planck introduced a nonscalar quantity, 7, 
which he identified with temperature, and which appears 
to be related to the present @ by the equation 7=06w 
Tolman introduces a quantity 
identical with 6. 


which appears to be 


‘The inequality 


(36) 


cannot be proved without some assumption 
concerning the relativistic form of Fourier’s law 
and Kelvin’s hypothesis. The simplest assump- 
tion leading to Eq. (36) is that 


= + 0Dus], (37) 


which satisfies Eq. (25) and reduces to the 
classical Fourier’s:law when ordinarily negligible 
terms are canceled (see below). The scalar 
k(20) is the thermal conductivity of the fluid. 
The term s%*(00/dx*) is obviously the relativistic 
temperature gradient. The term —@Dug has no 
classical analog, and implies an isothermal flow 
of heat in accelerated matter, in the direction 
opposite to the acceleration. It is ordinarily 
small, and may be explained as due to the 
inertia of energy. Equation (37) shows that 
q*=0 when the vector 00/dx*°+6@Duzg is parallel 
to the proper-time axis. This seems strange at 
first, but on reflection, it is seen to correspond 
to the fact that, in the classical theory, there 
will be no flow of heat if V@=0 even though 
00/dt~0—i.e., that there are adiabatic changes 
of temperature. 

As was discussed in I, Eq. (34) does not 
determine @ and 7» uniquely; however, if it -be 
assumed that Eq. (37) holds for two of the 
possible functions 6, @ and all possible values of 


s8(dv/dx*), s(dp/dx®), s®Dug, 
then it can be shown as in I that 
6/09=const. 


Thus Eq. (37) fixes the definition of the tempera- 
ture except for a choice of unit. 
From Eqs. (14) and (37), it follows at once that 


(38) 
It therefore remains to show that 
P°8(dug/dx*) 2 0, (39) 


the positiveness of \ and @ being assumed. For 
this purpose, the attention may be fixed on a 
single arbitrary point, A, in four-space. If the 
time-axis of the Galilean coordinate system is 
chosen parallel to the proper-time axis at A, 


and 
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then at this point, 


u*=1, 0, 0, 0, 
ua=—1, 0, 0, 0, (40) 
and 
000 0 
010 0 
— 
(41) 
0001 


Substitution into Eq. (32) shows that 


P%= 
while 


P™=2rc{ [(0u1/dx!) — (du2/ dx?) 
—(Au3/dx*)]}, (42) 


+ (du2/dx') ], etc., 


which are the classical expressions for the viscous 
stresses. It then follows, as in I, that Eq. (39) 
is true at the point A, and thus everywhere in 
four-space. 


THE MOMENTUM OF HEAT 


The general principle of the inertia of energy 
leads one to expect that a flow of heat g will be 
associated with a momentum g/c*. This result 
can be read out of the form of the energy- 
momentum tensor (Eqs. (22) and (24)), but it is 
interesting to consider the question in somewhat 
more detail. It will be supposed that the fluid is 
at rest in the Galilean system x*, so that Eqs. 
(40) and (41) are valid at all points of four-space 
and P*=0. Then Eq. (37) yields 


q’=0 
—kd6/dx*, x=1, 2, 3, 
while Eq. (2) becomes 
dm /dt=0 
so that m and v are functions of x', x*, x* only. 
The components of the energy-momentum 
tensor are 
W°=m(e+a), 
pb. 


The Eqs. (18) become 
m(d/dt) + (dg*/dx*) =0 (43) 


and 


(0/dt)(q*/c*) + (dp/dx*) =0. (44) 


Equation (43) is the equation of conservation of 
energy, and is (in this case only!) identical with 
the first law, Eq. (33). The Eqs. (44) are the 
momentum equations, and it is immediately 
seen that g*/c? is a momentum, since its time 
rate of increase is equal to the negative gradient 
of the pressure. 


Oum’s Law 


The general problem of thermoelectromagnetic 
phenomena is too complex for brief discussion, 
but it is of interest to indicate how these phe- 
nomena may be included in the present theory. 
Consider an ideal fluid whose internal energy is 
independent of the electromagnetic quantities, 
but which is an electric conductor. A molten 
metal approximates such a fluid, but shows 
thermoelectric effects, etc., which indicate a 
slight dependence of its internal energy on the 
electromagnetic variables. 

Let j* be the electric current density in 
e.m.u./cm? sec. =e.s.u./cm', and 


0 —E, -—E, -E, 

E 0 -B. B 

‘ap — =z z 
© 


E, —B, B, 0 


be the field tensor, the E’s being the electric 
field in e.s.u., the B’s the magnetic induction in 
e.m.u. The energy-momentum principle is now 
to be formulated as 


(46) 


The tensor W can be resolved into components 
as before, and the first law becomes 


— (ugW**) + We8(dus/dx*) 
= —jaugF*, (47) 


This may be transformed into an equation 
similar to Eq. (33), except that the right side is 
not zero, but the same as that of Eq. (47). 

The entropy and temperature may be intro- 
duced in exactly the same manner as above, 
and an equation, differing from Eq. (35) only 
by the additional term 


— (1/0) jaug Fe? 
can be deduced. The Eq. (36) requires this 
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term to be positive for all possible values of j. 
' and ug; this will certainly be the case if (with 
o>0 and pa scalar) 


oF (48) 


since Fu, is a vector in proper-space. However, 
Eq. (48) is precisely the well-known invariant 
form of Ohm’s law;’ in the case of stationary 
matter it reduces to the classical form: 


j*=p, cE,, cE,, cE;. 


Thus p is the charge density and o the electrical 
conductivity. 

The ease with which this law can be incorpo- 
rated into the equations lends strong support to 
the theory developed above. It will be noted 
that Maxwell’s equations have not entered this 
discussion. However, the Eqs. (47) and (48) 
are too few to determine all the variables, and 


My Tolman, Relativity, Thermodynamics and Cosmology, 
p.- 


the Maxwell equations are needed to complete 
the theory. 
CONCLUSION 

It has been shown that the classical considera. 
tions of I can readily be extended so that they 
are consistent with the special theory of rela- 
tivity. The further extension into the general 
theory should offer no special difficulties. 

One result has appeared clearly : It is necessary 
to introduce the current-density of matter sepa- 
rately from the energy-momentum tensor. This 
is also apparent from other treatments of energy 
and matter, but has rarely been emphasized, 
From the macroscopic standpoint, therefore, 
matter cannot be considered a form of energy, 
even though inertia is a property of energy 
rather than of matter. Matter has inertia because 
it has energy, but is not a form of energy. This 
same conclusion can be reached from a considera- 
tion of recent attempts to construct microscopic 
theories of the ultimate particles. 


Correction: The Thermodynamics of Irreversible Processes. II. 
(Phys. Rev. 58, 269 (1940)) 


Cart EcKART 


Ryerson Physical Laboratory, University of Chicago, Chicago, Illinois 


On page 274 of the above paper, the sentence beginning ‘‘It is not easy to 
see . . .”’ should be replaced by the following: 


It can be shown that 


20 


for all values of A;. Let 


Ae/RO, 
x" = Ar/ Ré, 


so that this sum may be written 


RY [exp (x,’) —exp (x.”) ]. 


Since the exponentials are monotonically increasing positive functions of 
their arguments, it follows that both brackets in the above expression always 


have the same sign. Since a,2 0, the required inequality follows at once. 
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LETTERS TO THE EDITOR 


Prompt publication of brief reports of important discoveries in physics may be secured by 
addressing them to this department. Closing dates for this department are, for the first issue of the 
month, the eighteenth of the preceding month, for the second issue, the third of the month. Because of 
the late closing dates for the section no proof can be shown to authors. The Board of Editors does 
not hold itself responsible for the opinions expressed by the correspondents. 


Communications should not in general exceed 600 words in length. 


Hyperfine Structure of Ionized Mercury Lines 


In connection with his investigations on the isotope shift 
in the band spectra of ionized mercury hydride, it has 
become necessary for Mrozowski to investigate the hyper- 
fine structure of certain ionized mercury lines. He has re- 
cently reported his results! regarding the analysis and the 
interpretation of the structure of some ionized mercury 
lines. Reference has been made by him to a short note in 
which the structure of \3984A had just been mentioned by 
Venkatesachar and myself;? but it is unfortunate that 
Mrozowski seems to be quite unaware of our subsequent 
detailed paper* on the “‘Hyperfine structure of some HglII 
lines,” in which the structures of \\3984, 2848 and 2262A 
had been investigated as early as in 1934. Mrozowski’s 
results in regard to the structure as well as the scheme of 
hyperfine levels of the line \3984A agree with our earlier 
findings. Mrozowski, just like ourselves, has been unable 
to confirm the satellite at +1.950 cm™ reported by 
Murakawa‘ earlier, and it is now certain that no such 
satellite exists. The total spin separations of the two levels 
and Sd* 6s* given by Mrozowski agree with 
those previously determined by us and are consistent with 
theoretical expectations based on Goudsmit’s formulae.’ 
The scheme of hyperfine levels worked out by us for ex- 
plaining the observed structure of the line is identical with 
the scheme proposed by Mrozowski; the agreement in the 
nuclear spin separations of the upper and lower levels is 
shown in Table I. 

In the case of \2848A, the analysis of the line has been 
rendered incomplete for want of sufficient resolving power. 
But it is certain that the line exhibits no even isotope shift 
and has, in agreement with theoretical expectations, a few 
close components arising from the odd isotopes of mercury. 
Large isotope shifts of the order of 0.680 cm™ between 
successive even isotopes have been observed by us in the 
level 5d* 6s 6 *Ds;2 as shown from the analysis of \2262A 


in mercury is 
Hg 199 Hg 201 
V. and S. 0.050 0.127 0.060 0.135 
Mrozowski 0.049 0.118 0.054 0.131 


It is abundantly clear that Mrozowski, who has analyzed 
3894A completely unaware of our earlier work, has lent 
valuable support to our conclusions. 


L. SrBparya 
of 


Bangalore, India, 
August 20, 1940. 
1S. Mrozowski, Phe. Bev. | 57, 207 (1940). 
2 B. Venkatesachar and L Sibaiya, Naturwiss. 52, 1041 (1931). 
3 B. Venkatesachar and L. Sibaiya, Proc. Ind. Acad. Sci. 1, 8 (1934). 
*K. Murakawa, Sci. of Inst. of Phys. and Chem. ‘Research, 


18, 299 ( 
*S. Goudsmit, Phys. Rev. 43, 636 (1933). 


On the Nuclear Gyromagnetic Ratios of the 
Chlorine Isotopes* 

The nuclear gyromagnetic ratios of the chlorine isotopes 
have already been determined' by the molecular-beam 
magnetic-resonance method. On the basis of the depths of 
the resonance minima the observed g value of 0.546 was 
assigned to the Cl® nucleus and that of 0.454 to CI*’. 
Although this assignment seemed reasonably certain it was 
thought worth while to eliminate any possible doubt on this 
point by experimenting with a molecular beam of an 
alkali chloride containing an enriched proportion of Cl*’. 

One of us (E.F.S.) has obtained,? by the use of the 
thermal diffusion method, a sample of HCI in which the 
abundance of Cl*’ was raised from the normal 24 percent 
to about 50 percent. A portion of this sample was used to 
prepare about 20 mg of LiCl. 

Two resonance minima, corresponding to the g values 
previously reported for chlorine, were again observed. The 
ratio of the depth of the 0.546 resonance minimum to that 
of the 0.454 minimum was found to be 1.2 for the enriched 
sample instead of 3.1 to 3.4 previously reported for normal 
samples. This constitutes conclusive evidence that the 
assignment has been correctly made. 

A recent determination* of the spin of Cl*? by the ob- 
servation of alternating intensities of lines in bands due to 
Cl*"Cl*" yields 5/2 as the most probable value of the spin 
of Cl*’. If we assume this value of the spin the moment of 
Cl" is 1.135 nuclear magnetons. 


E. F. SHRADER 
Yale University, New Haven, Connecticut 


Ss. M 
Columbia University and City College, New York, New York 
P. Kusca 
Columbia Usivenky, New York, New York 
October 25, 1 


* Publication ne Ernest Kempton Adams Fund for 
Research of Columbia University. 


Physical 
ve. Kusch and S. Millman, Phys. Rev. "56, 527 (1939). 
2 E. F. Shrader, Phys. Rev. 58, 475 (1940). 
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Identification of Two Radioactive Xenons from 
Uranium Fission 


The observation of two long-life radioactive xenons 
growing out of radioactive iodine from uranium fission has 
been reported by us! and confirmed by Dodson and 

_Fowler.? Recently some samples of cesium and barium were 
placed immediately behind a beryllium target bombarded 
by 16-Mev deuterons. After the bombardment the active 
cesium or barium was dissolved in water and the radioactive 
gas was then extracted into an ionization chamber. In the 
case of cesium only one radioactivity was observed which 
decayed with a half-life of 5 days, but in the case of barium, 
both 9.4-hr. and 5-day periods were present. The identity 
of the periods and of the absorption curves of these 9.4-hr. 
and 5-day radioactive gases from cesium and barium with 
those of radioactive xenon from uranium and thorium 
fission proves that they are the same substances. 

Cesium has only one stable isotope with a mass number 
of 133, and since its bombardment gives only one radio- 
active xenon, we assume that the 5-day xenon is produced 
by the following reaction 


ssCs'3(n, 


and assign to it the mass number 133. 

Although barium has seven known isotopes (130, 132, 
134, 135, 136, 137 and 138), only three of them would be 
able to produce radioactive xenon (except for isomers of 
stable nuclei) by a(m, a) reaction. These are 130, 136, and 
138. We have seen that Ba gives by a(n, a) reaction the 
5-day xenon and we interpret this as 


ssBa!**(n, a) sixes, 


Since the abundance of Ba™ is only 1/700 of that of Ba™, 
it is very likely that the barium of mass number 138 is 
responsible for the formation of the 9.4-hour xenon accord- 
ing to the reaction 

ssBa!™(n, 


Based on these results, the two chains found in uranium 
and thorium fission may be identified as 


60 min.— 531! 22 hr.—5Xe"™™ days—s5sCs' 


and 
sele'S~10 min.—s531 6.6 hr.—5;,Xe!™ 9.4 hr.—s5sCs™?. 


Iodine was also separated from the fast-neutron-activated 
cesium sample. The separated iodine decayed with a single 
period of 12.6 hr. This confirms the assignment? of 12.6-hr. 
iodine to 531 resulting from a reaction 


ssCs'3(n, ox) 


The author wishes to thank Prufessor E. O. Lawrence for 
his interest in this work and Dr. E. Segré for his helpful 
suggestions and criticisms. 

Cc. S. Wu 


Radiation Laborato 
University of California, 
Berkeley, California, 
October 17, 1940. 


C- Ss. 8. Raw. $83 ¢ 552 (1940). 


W. Dodson R. D. Fowler, aye, Rev . 57, 966 (1940). 
‘J. Livinecod and G. Seaborg, ys. Rev. $4, 775 (1938). 
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Vibration Analysis of the Chlorine Dioxide 
Absorption Spectrum 

The electronic band system appearing in the absorption 
spectrum of chlorine dioxide between 2600A and 5225A has 
been measured by several investigators'* and vibration 
analyses have been attempted. The author offers evidence 
in support of a new analysis. Urey and Johnston were the 
first to attempt a detailed analysis of this system. However, 
none of their frequencies for the ground state agreed with 
the two infra-red frequencies found later by Bailey and 
Cassie;* namely, 945 cm~! and 1109 cm™. For the third 
fundamental frequency of the ground state Bailey and 
Cassie accepted 528 cm which appeared as a constant 
difference in two pairs of strong progressions in Urey and 
Johnston’s analysis. Later Ku made more accurate and 
more complete measurements on the bands and attempted 
an analysis using the above frequencies but did not use the 
selection rules of Herzberg and Teller.® 

Using the data of Ku, one finds that the 528 frequency 
comes out to be either 528 cm™ or 524 cm™ depending on 
which pair of progressions is used in determining its value, 
This disagreement is greater than experimental error. The 
author proposes to replace the fundamental frequency of 
528 cm™ by 447 cm™ which appears as a constant differ- 
ence in three separate pairs of progressions. 

Photographs taken by the writer with the absorbing gas 
at dry ice temperature give clear evidence in support of this 
change. 

From electron diffraction experiments, Brockway’ finds 
that the apex angle of ClO; is 137°+15°. On this basis one 
must make the following assignment to the fundamental 
frequencies of the lower state: 


447 cm", 
ve’ = 945 
v3’’=1109 cm™, 


where »” designates the frequency of the bending vibra- 
tion, v2” that of the breathing vibration, and »;3"’ that of the 
antisymmetrical vibration. The apex angle of the lower 
state of ClO: calculated from these frequencies on the basis 
of the valence force model is approximately 122°. 

As for the symmetrical frequencies of the upper state, 
the author essentially agrees with Urey and Johnston as 
well as Ku. These are: 

cm", 
| 708 

The bands have been photographed by the writer in the 
second order of the 30-ft. 30,000-line per inch grating at 
the University of Chicago. Each band is made up of sub- 
bands, some of which show clearly defined rotational series. 
This rotational structure is now in process of analysis. The 
sub-bands, like the bands themselves, degrade to the red. 
This is sufficient to show that both dimensions of the mole- 
cule increase as it passes to the upper state. 

Further, the 708 frequency is excited up to 20 quanta or 
more, while no more than two quanta of the 290 frequency 


are excited with enough intensity to be observed. If 
qualitative use is made of the Franck-Condon principle, 


I 
t 
fc 
1! 
tl 
al 
0, 
4 gr 
Ca 
th 
th 

Th 
is 
as 
ord 
pic 
side 
7 
for 
adv 
Dr. 
‘2 
Lae 
(193: 
| 
A 
on ‘ 
tung 
Pott 
tion 
resol 
of A 


LETTERS TO 


the foregoing facts require the assignment 


= 290 
vo’ = 708 


In addition, if use is made of two series in which it appears 
that vs’ changes by two quanta, the present analysis gives 


v3’=790 


On the basis of these approximate frequencies, the angle 
for the upper state is found to be about 114° compared to 
122° in the lower state. This small change is consistent with 
the already mentioned fact that only v2’ progressions are 
appreciably excited. The numerical value of each frequency 
given above corresponds to a one-quantum jump from the 
0, 0, 0 state of vibration. 

All distinct bands excluding isotope bands are now 

grouped into twelve progressions. All other reported bands 
can aa classed as sub-bands, or are in the violet region of 
the spectrum where predissociation causes confusion. In 
the notation these progressions are: 


STRONG MEDIUM 
(0 0 0) © 1 0) 
a 0 (1 
2) WEaK 

( v2’ 2 0) 

MEDIUM (1 v2" 0) 
(O v2’ 0)—(1 1 0) 


The (0 0 - 0 0) band, contrary to previous analysis, 
is the 21,016 cm™ band. 

The fact that (0 v.’ 2)+-(0 0 0) and (0 v,’ 0)+-(0 1 0), 
as well as (0 v.’ 0)+-(0 2 0) are stronger than would 
ordinarily be expected is consistent with the qualitative 
picture of the Franck-Condon principle, taking into con- 
sideration the change in shape of the molecule. 

The author wishes to thank Professor R. S. Mulliken 
for suggesting the problem and for generous interest and 
advice throughout the work. The author is also indebted to 
Dr. H. G. Beutler for helpful experimental suggestions. 


r J. B. Coon 
Chicago University, 
Chicago, Illinois, 
October 25, 1940. 


1C. F. Goodeve and C. P. Stein, Trans. Faraday vy = 738 (1929) 

2H. C. Urey and H. Johnston, Phys. Rev. 38, 2131 ( 1931). 

3W. Finkelnburg and H. J. Schumacher, Zeits f. physik. Chemie 
Bodenstein-Festband, 704 (1931). 


‘C.R. Bailey and A. B. D. Cassie, Proc. Roy. Soc. A137, 622 (1932). 
*G. Herzberg and E. Teller, Zeits. f. physik. Chemie B21-—22, 410 
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' at O. Brockway, Rev. Mod. Phys. 8, 231 (1936). 


Remarks on the Temperature Dependence of the Work 
Function of Tungsten 


A few remarks in criticism of the paper by J. G. Potter! 
on “Temperature dependence of the work function of 
tungsten . . .” are called for, since in my opinion Dr. 
Potter has overestimated the reliability of the interpreta- 
tion of his work when he writes: ‘The coefficient obtained 
resolves the discrepancy between the experimental value 
of A in the Richardson equation and the theoretical factor 
of 120 without the introduction of a reflection coefficient.” 
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The so-called “experimental” value of A is taken as 60 
amp. per sq. cm per degree squared. Its determination 
depended on the assumption that the entire area of the 
filament studied had a constant work function and that this 
could be obtained by measuring the emission with a strong 
accelerating field and then extrapolating back according 
to the Schottky theory to zero in order to evaluate the 
“theoretical” emission which should obtain at zero field. 

Experiments by Johnson and Shockley,? Martin,*? and 
Nichols‘ demonstrate beyond any doubt that the emission 
from a tungsten wire is far from uniform and, in fact, prob- 
ably over 95 percent of all of the emission comes from less 
than 50 percent of the surface. Although no theoretical in- 
vestigation has as yet shown that an appreciable reflection 
effect is to be expected when very slow electrons pass over 
a potential barrier, al/ experiments designed to investigate 
the energy distribution of electrons emitted thermionically 
or photoelectrically show a very considerable deficiency in 
the number of slow electrons. The experimental fact is that 
the order of 50 percent or more of the electrons normally 
expected are apparently reflected at the surface and do not 
contribute to the observed emission. Either the effect of 
the non-uniform emission or the reflection effect alone 
would be capable of “resolving’’ the discrepancy between 
the “experimental” value of A =60 and the theoretical 
value of A = 120. A positive temperature coefficient of the 
work function as reported by Dr. Potter makes matters 
worse when all of the recognized experimental results are 
considered. 

The negative temperature coefficient of the work- 
function as determined by Langmuir® has been shown® to be 
capable of resolving the difficulty and yet it is necessary to 
stress the fact that this problem cannot be said to be 
“resolved” until every step in the argument has been 
checked and cross-checked by independent experiments 
since the experimental difficulties are known to be so great. 

Dr. Potter leads the reader to think that Langmuir 
attributed his observed negative temperature to “‘effects 
resulting from imperfect vacuum conditions.” A careful 
reading of Langmuir’s paper shows that the more likely 
source of error was to be attributed to possible alterations 
in the exact relative positions of his tube elements as the 
temperature was changed. As a result of a recently com- 
pleted investigation of field emission by J. H. Daniel,” it 
has been possible to demonstrate that even under the most 
favorable vacuum conditions, which can be obtained only 
with sealed-off and properly ‘‘gettered” tubes, a surface of 
tungsten cannot be maintained at temperatures below 
900°K more than a very few minutes without considerable 
alteration of the work function because of the adsorption 
of gas. Most experimenters in this field know that serious 
adsorption takes place many times more rapidly when the 
experimental tube is not sealed off from the pumping system 
even though the pressure as measured by an ionization 
gauge is of the order of 5X10-* mm. Therefore it is not 
surprising that Dr. Potter experienced so much difficulty 
in this respect. At this pressure a monolayer could be 
formed in about five minutes and only a fraction of a mono- 
layer is needed to make a very considerable alteration in 
the work function. This would confuse any measurements 
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made on the temperature coefficient of the work function 
of tungsten over the range of temperature below 900°K, 
where the effect sought amounts to only 0.04 volt. 


W. B. NotTInGHAM 
George tory of Physics, 
Massachusetts Institute of Technology, 
Cambridge, Massachusetts, 
October 18, 1940. 
:G Potter, Phys. Rev. 58, 623 (1940). 


: Wr Stoke 


H. Daniel, 


Reply to Remarks on the Temperature Dependence of the 
Work Function of Tungsten 


In the paper' criticized in the letter of W. B. Nottingham? 
the statement, that the positive temperature coefficient 
obtained for the work function of tungsten “resolves the 
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discrepancy between the experimental value of A in the 
Richardson equation and the theoretical factor of 129 
without the introduction of a reflection coefficient,” does 
not imply that the coefficient is uniquely capable of resoly- 
ing the discrepancy. As cited by Dr. Nottingham, interpre. 
tations of thermionic emission phenomena are possible on 
the basis of several other lines of experimental evidence, 
and more research will be necessary before a decision can 
be made in favor of any one interpretation. 

In the paper criticized, the effects of adsorbed gas were 
taken into account, whereas in the paper by Langmuir,? in 
which it was stated that the vacuum was not “‘unimpeach- 
able,” no criterion was established for distinguishing be- 
tween true temperature effects and contamination effects, 


James G. Porter 
it of 
th Debeta State School of Mines, 
Rapid City, South Dakota, 
ober 28, 1940. 


Fetter. Phys. Rev. 58, 
B. Langmuir, Phys. —, 49, 428 (1936). 
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